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The lattice constants of the binary alloys of silicon 
with silver, gold, lead, tin, zinc, cadmium, antimony and 
bismuth have been measured. It is shown that the lattice 
constants of the elements in these alloys are practically 
identical with those of the pure elements. This fact leads 


to the interpretation that mutual solubility of the elements 
in these alloys occurs to only a very slight extent. It is 
further observed that silicon forms no compounds with 
any of these elements. 


R each of the eight binary silicon alloy 
systems included in the present report, 
analogy with silicon alloys of closely related 
elements and, in at least one case, the earlier 
results of thermal and microscopic analysis would 
lead us to expect compound formation or ap- 
preciable ranges of solid solution or both. It was, 
therefore, surprising that not a single case of 
compound formation was observed in these 
alloys, and, further that in none of these systems 
was there a range of solid solution of sufficient 
magnitude to be detected by the method applied. 
The method used in the present investigation 
was the determination of the lattice constants 
by means of the x-ray crystal structure analysis. 
The usual procedure for determining the limit 
of a solid solution range is to measure the lattice 
constants at several different compositions within 
the homogeneous solid solution range and then 
for the same phase when in equilibrium with the 
neighboring solid phase. By projecting this last 
value to the curve for the homogeneous range, 
the equilibrium composition can be established, 
frequently with a high degree of accuracy. 
Aside from the precision attainable in the 
measurement of the lattice constants, the accu- 


racy of the method depends upon the difference 
in radii of the atoms concerned and the actual 
extent of the solid solution range. These two 
factors, however, are not entirely independent; 
general experience shows that if the atomic radii 
differ by more than about 10 percent, the solid 
solution range will be limited. The limit seems 
to be set by the amount of distortion caused by 
the foreign atoms which the lattice can with- 
stand. This qualitative rule is subject to further 
restrictions due to influences which would cause 
positive deviations from Raoult’s law or give 
rise to compound formation. 

In the present cases, the usual procedure could 
not be applied because of the entirely insignifi- 
cant changes in the lattice constants observed 
even when the neighboring solid phases were in 
equilibrium. The cameras used were of the 
Bohlin type as modified by Phragmén.'! They 
were calibrated empirically with NaCl as the 
standard. Only the most deviated reflections of 
radiation from an iron target were used for 
calculating the lattice constants. The error for 
substances having cubic symmetry is certainly 


1Westgren, Trans. Inst. Min. Met. Eng.; Inst. of 
Metals Div. 1931, p. 13. 
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not greater than +0.001A; in the tetragonal 
and hexagonal systems, +0.001—0.002A for the 
smaller constant and possibly +0.003—0.004 
for the length of the c axis when this is more 
than 10A. 

As a first approximation, which is certainly 
sufficient for the present cases, we may assume 
that the change of the lattice constant is a linear 
function of the atomic composition and further 
that it may be calculated from the difference in 
the radii of the atoms concerned. The differences 
between the radius of the silicon atoms and the 
radii of the other atoms varies about 15 to 50 
percent. Calculation shows that with the accu- 
racy attainable in determining the lattice con- 
stants, solid solubility of one element in another 
in these alloys, amounting to 0.2 atom percent, 
should have been indicated with certainty, while 
for Si in Pb, Zn or Cd, as well as for Pb or Bi 
in Si the detectable limit is about 0.1 atom 
percent. As already noted, no solubility of this 
order was observed. 


MATERIALS 


Silicon of exceptional purity was supplied by 
the Electro Metallurgical Company through the 
courtesy of Dr. F. M. Becket. It had the follow- 
ing analysis: Si, 99.79 percent; Fe, 0.027; Al, 
0.008; Ca, 0.014; C, 0.004; Os, 0.002; Ne, 0.1; 
Hg, trace. 

Bismuth: Merck’s C.P. with total maximum 
impurities of approximately 0.04 percent. 

Antimony and cadmium: Kahlbaum “Zur Ana- 
lyse.” 

Gold and silver: Specially refined material from 
the U. S. Metals Refining Company, in 
each case 99.99-+ pure. 

Lead: Merck’s C.P. with total maximum im- 
purities amounting to 0.1 percent. 

Tin: Baker’s ‘‘Analyzed’”’ with total metallic 
impurities of approximately 0.03 percent. 

Zinc: Eimer and Amend, “Tested Purity” with 
total metallic impurities amounting to ap- 
proximately 0.2 percent chiefly Pb and Cd. 


PREPARATION OF ALLOYS 


For each of the systems, alloys were made by 
at least two different methods. It was early 
realized that the miscibility of the two elements 
and the formation of sound and uniform ingots 
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might be seriously interfered with by the 
presence of thin oxide films, especially on the 
silicon. Accordingly it was attempted in each 
system to make at least one alloy at tempera- 
tures above the melting point of silicon, ice., 
above 1420°C. This was successful except in the 
cases of zinc and cadmium. These elements when 
added to liquid silicon in a hydrogen atmosphere, 
boiled off completely; the alloys were, therefore, 
prepared by heating mixtures of the two elements 
just below the boiling point of zinc and cadmium 
respectively for several days. 

In most cases, weighed amounts of the two 
elements were sealed off in an evacuated quartz 
tube, heated to the desired temperature and 
quenched; in the remaining cases melting was 
conducted in a hydrogen atmosphere. Micro- 
scopic examination of the ingots showed a 
relatively uniform distribution of the phases in 
each case. Filings of these ingots were then made 
and annealed in vacuo for several days or, in 
some cases, weeks to allow equilibrium to be 
reached. These final annealing temperatures are 
recorded in Table I. It is to be emphasized that 
every care was exercised, (1) to maintain the 
original purity of the materials, (2) to provide 
adequate opportunity for any reaction or solution 
phenomena to take place, and (3) to permit 
final equilibrium relations to be established. 


RESULTS AND DISCUSSION 


The experimental results are summarized in 
Table I. All results are included except those 
from photographs on which the lines did not 
permit accurate measurement or of alloys which 
had separated into two distinct layers before 
solidification. The main result has already been 
mentioned, namely, that no clear-cut evidence of 
solid solution was found and for none of the 
systems examined did the films contain lines 
other than these corresponding to the two 
elements of the system. 

Arrivant? found silver to dissolve in solid 
silicon to the extent of 10 percent. The maxi- 


2 Arrivant, Comptes Rendus 147, 859 (1908). Cf. L- 
Loskiewicz, Congr. Intern. Mines (Liege 22-28 June, 
1930). 

3 According to J. L. Haughton (personal communicatio®) 
the original thermal diagram by Arrivant has been it 
correctly reproduced in Int. Crit. Tab. 2, 422. 
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X-RAY STUDY OF BINARY ALLOYS OF SILICON 


TABLE I. 
Alloying %Si_ Final anneal Lattice Constants 
element Temp. °C Alloying element Silicon 
Pure Si 5.4170 
Ag Pure Ag — 4.0765 — 
5.0 825 4.0769 5.4174 
40.0 635 4.0765 5.4176 
40.0 635 4.0768 5.4186* 
40.0 635 4.0771 5.4173 
40.0 825 4.0769 5.4181* 
50.0 770 4.0764 5.4167 
70.0 635 4.0773 5.4175 
80.0 825 t 5.4176 
85.0 750 4.0761* 5.4174 
Pb Pure Pb 4.9385 
50.0 130° 4.9375 5.4156 
50.0 150° 4.9382 5.4160 
Au Pure Au 4.0693 
50.0 240° 4.0682 5.4176 
a c 
Sn Pure Sn _ 5.8167 3.1728 — 
12.0 595° 5.8156* 3.1744* 5.4165 
50.0 150° 5.8178 3.1728 5.4184 
80.0 595° 5.8151* 3.1748* 5.4174 
Zn Pure Zn — 2.6581 4.9341 — 
50.0 273° 2.6578 4.9353 5.4164 
Cd Pure Cd 2.9713 5.6046 
50.0 yi 2.9709 5.6058 5.4181 
Sb! Pure Sb 4.295 11.247 
40.0 595° 4.296* 11.247* t 
50.0 580° 4.297 11.243 5.4176 
80.0 595° 4.296 11.250 5.4182 
Bi! Pure Bi 4.535 11836 
50.0 240° 4.534 11.830 5.4161 


* Value uncertain due to poor lines. 
} Lines present but too poor to measure. 


‘Reported in terms of hexagonal axes instead of usual 
rhombohedral axes. 


mum solubility from the results of Table I would 
be of the order of 0.2 percent by weight. 

For the antimony alloys, R. S. Williams‘ found 
a solid solution of antimony in silicon up to 
1 percent by weight and of silicon in antimony 
up to 0.3 percent. He also found that bismuth 
dissolved up to 0.8 percent of silicon. The 
Present results indicate that the solubility of 
antimony in silicon is under 0.5 percent by 
Weight. Rough calculation shows that the pres- 


. ence of 0.3 percent Si in antimony or 0.8 percent 


in bismuth would change the lattice constants 


*R. S. Williams, Zeits. f. anorg. Chemie 55, 1 (1907). 
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of these two elements by approximately 0.25 
percent and 1.5 percent respectively and it is 
evident from Table I that no such change has 
occurred. A probable reason for the difference 
between the older results and the present is the 
greater purity of the silicon used in this investi- 
gation. The silicon used by Williams, for ex- 
ample, contained 0.95 percent Fe and 0.4 percent 
Al. 

The failure to find any solubility in the Au-Si 
system agrees with the work of Di Capua.® 
Tamaru’s® similarly negative results in Pb-Si 
and Sn-Si are also confirmed. 

Note added in proof: In view of the facts (1) 
that the preparation of silicon of a high degree 
of purity has been developed orly during the 
last few years and (2) that many of the binary 
alloy systems of metallic elements with silicon 
have not yet been systematically investigated, it 
is of some interest to inquire what metallic 
elements may be expected to form compounds or 
extensive solid solution with silicon. In this 
connection it is interesting to note that silicon 
has been observed to form compounds with all 
of the first twenty elements of the periodic 
system except its immediate neighbors Al and P, 
the rare gases, and the two alkali metals Na 
and K. In the case of Al a small range of solid 
solution of silicon in aluminum is observed at 
high temperatures. The alloys of the alkali and 
alkaline earth metals with silicon, excepting Li 
and Ca, appear not to have been studied. Of the 
metallic elements remaining in the periodic 
system, compound formation or extensive solid 
solution may be expected only in the alloys with 
transition elements, i.e., elements having in- 
complete inner electron shells (atomic numbers 
21-28, 39-46, 57-78). The only thoroughly 
established exception to this statement known to 
the writers is copper and it has been observed by 
other workers that this element in alloys some- 
times exhibits the behavior of a transition 
element. The basis for these generalizations was 
found chiefly in the results of x-ray structure 
analysis. 


5 Di Capua; Atti della reale accademia nazionale dei 
Lincei 291, 111 (1920). Cf. L. Loskiewicz, reference 2, who 
finds that Si dissolves to a limited extent in Au. 

6S. Tamaru, Zeits. f. anorg. Chemie 61, 40 (1909). 
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Oil Lenses on Water and the Nature of Monomolecular Expanded Films 


IRvinG LANGMurIR, General Electric Company, Schenectady, New York 
(Received September 8, 1933) 


A higher hydrocarbon, such as tetradecane, on the 
surface of water forms a circular lerts. The condition for 
the formation of such stable lenses is that the spreading 
coefficient Fgs=yi1—y2—vyi2 shall be negative. As the 

‘volume of the oil lens is increased, the thickness ¢ at the 
center slowly approaches a limiting value given by ¢,,? 
= —2F,;/gp2(o1—p2). Equations are given by which Fs 
can be accurately determined from measurements of the 
radii of large lenses of known volume. The magnitude of 
the linear tension at the lens boundary is also calculated. 
Experiments with tetradecane on water at 25°C give 
Fs=—6.2 dynes cm™ and y=26.9. When hydrophilic 
molecules are introduced into the interface between the 
lens and the water, they give a spreading force F,_. which 
causes a decrease in —Fg and in ¢t. When ¢<0.1 mm 
gravitational effects are negligible and the lens degenerates 
into a duplex film. If such a film is confined by a movable 
barrier (piston), the force per unit length F exerted on 


the barrier is equal to Fs. If the adsorbed molecules at 
the interface give a force Fi. in accord with the gas laws, 
then (F—F )(a—ao)=kT gives the equation of state of 
the duplex film, where Fy is the value of Fs in the absence 
of hydrophilic substances. It is shown that this equation 
applies to expanded oil films on water and that these are 
properly regarded as duplex films although they are 
monomolecular. The presence of a kink in the observed 
F—a curves-is attributed to the sudden appearance of 
micelles, each containing §-molecules when rises to a 
critical value. Fer the fatty acids on dilute HCl, 6=13. 
This small value is due to the crowding apart of the heads 
because of hydration. The esters show little or no hydration 
and @=60. The heat of formation of the micelles in the 
fatty acids is found to be 5.5+0.24n kg cal. per mole 
where » is the number of carbon atoms. The values of 
Fy and ao for various substances are determined from data 
given by Adam. 


ATTY acids and similar organic substances 
spread over the surface of water to form 
monomolecular films in which the active groups 
or heads are in contact with the water while the 
hydrocarbon chains or tails form a layer above 
the plane containing the heads. When such films 
are compressed by a force F (dynes per cm), they 
form condensed films having an area 20.5A? for 
each molecule. This area is identical with the 
cross-sectional area of the molecules in stearic 
acid crystals in a plane parallel to that occupied 
by the heads of the molecules.' Since this plane 
makes an angle of 63° 38’ with the axis of the 
hydrocarbon chains, there is strong evidence that 
the hydrocarbon chains in these condensed films 
of normal fatty acids on water are inclined about 
27° from the vertical. 

With the lower fatty acids, such as lauric and 
myristic acids, when the films are only slightly 
compressed, they occupy much larger areas than 
20.5A*. There are many different substances 


1 Alex. Miiller, Proc. Roy. Soc. (London) A114, 542 
(1927). 


756 


which give, under these conditions, a rather 
definite area of about 50A? per molecule. Adam 
has called these films liquid expanded films. The 
tendency to form these expanded films increases 
as the hydrocarbon chain is shortened, or as the 
temperature is raised, the effect of each additional 
CH, group in the chain being equivalent to a 
lowering of temperature of about 8°C. 

In two papers?:? in 1925 I offered an expla- 
nation of the existence of these expanded films. | 
pointed out that in these films there are really 
two surfaces which form the upper and lower 
boundaries of the hydrocarbon part of the film. 
Since the thicknesses of these films are several 
times that of the hydrocarbon chain, the region 
between these boundaries must be regarded as 
essentially a typical hydrocarbon liquid in which 
the molecules possess all the natural freedom of 
motion of such liquids with the single restriction 


21. Langmuir, Colloid Symposium Monograph 3, # 
(1925). 

3]. Langmuir, an article, The Effects of Molecular 
Dissymmetry on the Properties of Matter, in a book entitled 
Colloid Chemistry, Jerome Alexander, 1, 525 (1926). 
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OIL LENSES AND MONOMOLECULAR EXPANDED FILMS 


that the molecules are constrained to keep one 
end in contact with the lower surface. A higher 
normal hydrocarbon liquid on a water surface 
does not spread; it must therefore require a 
definite spreading force F to pull a hydrocarbon 
liquid out over a water surface to form an 
expanded film. Such a force must have its origin 
in some property imparted by the heads. At the 
interface between the hydrocarbon liquid and the 
underlying water there are, in the case of a fatty 
acid film, carboxyl groups. Even if these do not 
exert repulsive forces on one another, they should 
at least tend to spread within the available 
surface as a two-dimensional gas, and thus exert a 
two-dimensional gas pressure F which will 
counteract the normal tendency of the hydro- 
carbon layer to draw together. 

It was shown that this theory leads to an 
equation of state of the form 


(F—Fy)(a—ay) =kT. 


By taking the values Fy = —13 and ao= 18, this 
gave very satisfactory agreement with Adam’s 
early data** for myristic acid at 32.5°. This theory 
as given in the papers referred to was presented 
briefly in about four pages and there was no 
detailed analysis of the experimental curves. 

Adam has definitely rejected this theory.‘ For 
example, in his book he says: 


Langmuir suggested that the molecules adhere by 
their chains, and that the heads are in a state of constant 
agitation, the chains tending to make a condensed, and 
the heads a gaseous film. But since it is not the length 
of the chains, but the nature of the heads, which deter- 
mine whether the films are liquid or vapor expanded, 
and since those heads which probably have the greatest 
attraction for one another laterally form liquid expanded 
films, a better suggestion seems to be the reverse of 
Langmuir’s, that the heads tend to hold the molecules 
together, while the chains try to disrupt the film. 


Adam considers, however, that “the detailed 


structure of the liquid expanded films is difficult 
to decide.”’ 


Rideal’ offers no explanation of expanded films 


° N. K. Adam, Proc. Roy. Soc. A101, 516 (1922). 
N. K. Adam, Physics and Chemistry of Surfaces, 
Clarendon Press (1930), page 75. 
SE. K. Rideal, Introduction to Surface Chemistry, 
Cambridge University Press, pp. 101-104. 
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and does not accept those previously proposed by 
the writer. 

The evidence in support of the theory proposed 
is, however, so striking that it needs to be 
presented in more detail than has been done 
previously. It is clear that Adam has mis- 
understood the theory and does not see its 
widespread application to his own experimental 
data. A recent more thorough examination of 
Adam’s published results has led me to extend 
the theory that I proposed in 1925 so that it now 
seems to offer a nearly complete theory of all 
observed characteristics of expanded films. 


CHARACTERISTICS OF TYPICAL 
MONOMOLECULAR FILMS 


The complete curve giving F as a function of a, 
the area per molecule, such as that represented 
by Fig. 1, represents essentially an equation of 
state for a 2-dimensional system which may 


Fic. 1. Diagram of F-a ome for typical monomolecular 
ms. 


contain gas, liquid and solid phases. The lines 
shown in Fig. 1 may not all be observable with 
any single substance for reasons that will be 
discussed later. Since only positive values of F 
are directly observable, all those portions of the 
curves which correspond to negative values in 
Fig. 1 can be inferred only from theoretical 
considerations. 

The line AGH corresponds to a typical con- 
densed film such as that of stearic acid on pure 
water at room temperature. When such fatty 
acids are placed on water containing dilute 
hydrochloric acid, a break occurs at G so that the 
curve below G is given by BG. Adam believes 
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that the area corresponding to B is that of the 
heads of the molecules. 

With a liquid expanded film such as myristic 
acid on dilute hydrochloric acid, the observed 
F—a curve leaves the line GB near R and gradu- 
ally flattens out as it approaches J where there is 
a sharp kink, and it then falls more rapidly along 
a line JK to a definite point M. It then extends 
horizontally along the line MN until at some 
point WN it begins to fall again approaching the 
horizontal axis asymptotically as a continues to 
increase. The theory proposed in 1925 indicated 
that the line JK is part of a rectangular hyper- 
bola in which the horizontal asymptote ST lies at 
a distance Fy below the horizontal axis OX, and 
the vertical asymptote is represented by the 
dotted line at V corresponding to an abscissa dp. 

It is the plan of the present paper to develop a 
theory which gives the mechanism underlying 
each part of the broken line shown in Fig. 1, and 
by which quantitative information regarding the 
fundamental properties of molecules may be 
obtained. 


LARGE FREE LENSES OF OIL ON WATER 


Let us consider the conditions for equilibrium 
between a pure water surface (I) and a thin 
layer of an oil (II) placed upon the water. We 
assume provisionally that the film, although of a 
thickness large compared to that of the oil 
molecules, is yet so thin that gravity plays no 
appreciable part in causing it to flow out over the 
surface. We are then surely justified in assuming 
that the oil film has an upper surface (surface 
tension y2) and a lower surface (interfacial 
tension 712) which do not influence one another. 

For equilibrium between this thin oil film and 
the water it is necessary that the surface tension 
71 of the water shall balance the sum of y2 and 
‘iz for the film. Let us put 


(1) 


If Fs>0, that is, if y1 is greater than yo+712, 
the oil spreads over the water until it covers the 
whole surface or until the film becomes so nearly 
monomolecular that ye2+vyi2 increases to the 
value 71. 

If Fs<0, the oil film contracts and gathers 
into one or more lenses of such thickness that 
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gravity can prevent any further coutraction. 
Harkins has called Fs the spreading coefficient. 

Harkins’ data® for normal chain hydrocarbons 
on water show that Fs is +3.3 for hexane, 1.7 for 
heptane, and 0.3 (dynes cm) for octane. No 
data seem to be available for yi2 or Fs for the 
higher liquid hydrocarbons. If, however, the rate 
of decrease of Fs per CH observed by Harkins 
continues beyond octane, the value of hexadecane 
should be F's = —11.7. For paraffin oil (Squibbs’ 
petrolatum) which probably consists mostly of 
saturated cyclic hydrocarbons, Harkins* has 
found Fs = — 13.69 at 20°C. 

A simple direct method for determining Fx for 
nonspreading substances consists in measuring 
the diameters D of very large lenses formed when 
known volumes V of oil are placed on surfaces of 
water. Except when the density of the oil is very 
close to that of water, the curvature of the 
interface between the oil and the water becomes 
negligible at a distance of 2 or 3 cm from the 
edge of the lens. For lenses of a diameter greater 
than about 6 cm, the central part of the lens is 
thus of a uniform thickness. For a series of lenses 
of increasing radius R, the central thickness / 
slowly approaches a limiting value ¢,, which may 
be readily calculated by the following method. 

Let z be a coordinate which measures vertical 
distances upward from the plane of the water 
surface (far from the boundary of the lens). Then 
since the oil is lighter than water (and only such 
oils can form large lenses), Z2, the level of the 
upper surface of the lens, will have a positive 
value. The pressure at the lower surface (the 
interface) at Z,. must be the same as in the 
water phase at the same level. This condition 


gives 
Z2=t(p1—p2)/p1 (2) 


and 
Z12= —tp2/pr, (3) 


where p; is the density of the water and p: that 
of the oil. 

Consider now a lens so large that the curvature 
of the boundary is negligible. The total force per 
unit length acting across a portion of a vertical 
plane cutting the lens far from the boundary 
consists of the surface tension forces yi1+712 a0 


¢W. D. Harkins, Colloid Symposium Monograph 6, 2° 
(1928). 
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OIL LENSES AND MONOMOLECULAR EXPANDED FILMS 


the hydrostatic liquid pressures. From the free 
surface down to the depth y in any liquid this 
latter force is the pressure integral /p dy and is 
thus equal to (1/2)gpy’. Thus down to the 
interface at Z1. the total spreading force per unit 
length in the central part of the large lens is 


Fy = (1/2)g pet? —y2— 12. (4) 


For equilibrium, this must equal the spreading 
force F; in the water, remote from the lens 
boundary, down to the same depth Zy». This 


force is 
(5) 


Equating and F,, eliminating by Eq. 
(3), and combining with Eq. (1) we have 


(6) 


From a knowledge of ¢,, we can thus calculate 
the spreading coefficient F's. The actual thickness 
t at the center of a lens, however, often differs 
from ¢,, by as much as one percent even for lenses 
as large as 60 cm in diameter. This is due to the 
fact that there is a linear tension f (a one- 
dimensional analogue of pressure) along the 
perimeter of the lens acting to keep the lens in a 
circular form. This force, which is of considerable 
magnitude (6.5 dynes for paraffin oil on water) 
must be balanced by an increase in the spreading 
force F,— F, by means of an increase in ¢, just as 
the surface tension in a soap bubble is balanced 
by the increased air pressure within the bubble. 
This effect can be taken into account by sub- 
tracting a term f/R from the 2nd member of Eq. 
(4) before equating to F;. In this way we obtain 


12=1,2(1—f/FsR). (7) 


CALCULATION OF THE LINEAR TENSION OF A 
CYLINDRICAL LENS BOUNDARY 


to? = —2F sp1/gp2(pi— p2). 


A plane through the vertical axis of symmetry 
of a large lens cuts the free surfaces and the 
interface in the manner illustrated in Fig. 2. As R 
is made to increase indefinitely, the surfaces 
represented by the meridian curves in this figure 
become noncircular cylinders with horizontal 
elements. 

Consider now the horizontal forces acting 
across a vertical plane which is perpendicular to 
the meridian plane, but near the lens boundary. 
Let a be the angle which the tangent to any 
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Fic. 2. Profile of a lens on water. 


meridian makes with a horizontal plane and let y 
be the vertical distance of any point from the 
horizontal plane at which a=0. 

If we take the spreading force F (associated 
with any surface) to be zero when the surface is 
horizontal (a=0), then we have for the spreading 
force F, in a radial direction 


(8) 


For equilibrium, with a truly cylindrical lens 
boundary, F, must be constant and hence equal 
to zero. Thus putting F,=0 we obtain 


F,=~(1—cos a) — 


(9) 


where a is the Laplace capillary constant defined 
by 


y? =a*(1—cos a), 


a?=2y/gAp (10) 


and Ap is the difference of density on the two 
sides of a phase boundary. Eq. (9) may also be 
written 

(11) 


An equation of this form holds for each of the 
three surfaces. 

Now consider the horizontal forces acting © 
across the meridian plane near a lens boundary. 
Here the surface tension force per unit of hori- 
zontal distance is y/cos a (instead of y cos @ in 
the former case) and thus we have 


y = 23a sin (a/2). 


Fr =y(1—1/cos a) — 3gpy?. (12) 


Eliminating y by using Eq. (8) after putting 
F,=0 we have 


Fr=y(cos a—1/cos a)=—ysin?a/cosa. (13) 


There is thus a negative spreading force 
(tension) in a horizontal direction which is 
perpendicular to that of the slope of the surface. 
The tension f corresponding to any one of the 
three surfaces that meet at the lens boundary at J 
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is obtained as an integral f{Frdx. We may put 
dx = (dx/dy) (dy/da)da, 
in which dx/dy=cos a/sin a and 
dy/da=a sin a/a(1—cos a)}. 
With the value of Fr from Eq. (13) we thus find 


f=(ya/2) f sin? a(1—cos a)~Mda 


or 


f =(2(2)4/3)ya[1—cos*(A/2)]. (14) 


Here the integration has been carried from a=0 
to A, the angle a at the junction J where the 
three surfaces meet. The total linear tension 
acting along the perimeter of the lens is the sum 
of three terms 


f=fitfetfie (15) 


found from Eq. (14) by inserting the appropriate 

values of a and A. Since A is a rather small angle 

in most cases, it is useful to expand Eq. (14) in 
powers of A obtaining 


f = (24/4) (16) 


THE MERIDIAN CURVES AT A LENS BOUNDARY 


The effect of the linear tension in modifying 
the profile curves shown in Fig. 2 for lenses of 
moderate size may now be determined. A 
horizontal plane through the junction J allows us 
to apply our equations separately to the three 
phase boundaries. Consider, for example, tne free 
water surface. For a lens of moderate size the 
spreading force Fr given by Eq. (8) is not 
everywhere zero as in the cylindrical case but at 
any point P where a differs appreciably from 
zero, it must have a value which balances the 
effect of the linear tension of the parts of the 
surface outside of P. We may thus put 


Fr=fi/Ri 


where R, is an effective radius at which the 
tension may be regarded as being concentrated. 
In Fig. 2 the positions of R:, Re, and Riz have 
been chosen so that the two shaded areas (shown 
only for Riz) are equal. This assumption gives for 
moderate values of A 


Ri =R-+a,/2}. 


(18) 


(17) 
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Combining Eqs. (8), (17) and (10) we thus 
obtain’ 
ai—fi/yiRi (19) 


in place of Eq. (9). It should be noted that the 
value of f; in this equation is to be obtained from 
Eq. (14) by replacing A by ai. Similarly for the 
other two surfaces we have 


(20) 


and 
(Zo —212)2/ar2? =1-—cos (21) 


where R, and Rj, are obtained from equations 
like (18) except that the sign of the last term is 
reversed. 

If we apply Eqs. (19), (20) and (21) to the 
junction J in Fig. 2, we obtain’ 


(22) 
(Z2—2%0)*/a2* = 1—cos (23) 
(9 —Z12)?/a122=1—cos (24) 


Expanding the first members of the last two 
equations, combining with the first equation so 
as to eliminate 2, and eliminating and by 
Eqs. (2) and (3) we obtain 


p2/pi1= — 12) 
cos A1—2 cos A2—712 COS A 12) 


For equilibrium between the horizontal com- 


ponents of the surface tensions at the junction J 
we must have 


(25) 


1 cos Ajo. (26) 


Using this relation and introducing Fo by Eq. (1) 
into Eq. (25) we find 


t?=[21/gp2(p1— p2) JL— Fs 
J; 
which agrees with Eqs. (6) and (7) except that 


the effective radii and Ri, are now n0 
longer identified with R. 


(27) 


7 If the values of f; and by Eqs. (14) and (18) at 
inserted in these equations they may readily be shown t 
be identical with a rather complicated equation given by 
Poisson in his book Nouvelle Theorie de I’ Action Capillait 
as quoted in the Handb. der Exp. Physik VI, 98 (1924): 
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VOLUME OF LENS 


If the lens is of sufficient size to give a uniform 
thickness ¢ in the central part then the volume V 


is 
V = rR*t— sin Ai2—d2? sin Ae). (28) 


This may be understood readily from the fact 
that the vertical component y sin A of the 
surface tension at the junction multiplied by the 
perimeter 27R must be the buoyant force 
pgAV of the displaced volume AV resulting from 
the deviations of the surface from the horizontal. 


ConTACT ANGLES AND THE VALUES OF A 


The relation of the angles A;, Az and A12 at the 
junction J to the contact angles is given in Fig. 3. 
The angles subtended by the three phases at the 


Fic. 3. Contact angles at a lens boundary. 


junction J are B,, By and 7 — B;. These angles 
are related to the three surface tensions by 


cos (29) 
B.o=B,+Bs3. (31) 


Combining these equations with Eq. (1) we 
have 


sin? (B,/2) = —F (272+ Fs)/4yiv12, 
sin? (33) 


Thus the angles approach zero as —Fs de- 
creases to zero. From Fig. 3 we see that 


A,=A,—B;, (34) 
Ay=A;+Bi. (35) 


The value of A; may be obtained from Eqs. 
(22) to (24). Taking A to be a sufficiently small 
angle so that (1/24)A? can be neglected in 
comparison with unity, we may expand 1— cos A 
in Eqs. (22) and (23). Introducing the values of f 
from Eq. (16) after inserting the value of Ry, etc., 
from Eq. (18) and eliminating Z, and Zi, by 
Eqs. (2) and (3), we may obtain from Eqs. (22) 
and (23) the following 


2420, (36) 


beA2= —25[(p1— p2/p1)t—20], (37) 
where 


and 


be = + (24/4)a2/R+ (3/16)ae?/R?]. (39) 


Eliminating zo we obtain a single equation in A, 
and A». By Eq. (34) we may then eliminate A» 
and obtain 

b2Bs—(p1— p2)2't/ 


bit+be 


1 


(40) 


MEASUREMENTS OF THE DIAMETERS OF LENSES 


By placing measured volumes (or weights) of 
oil on a water surface and measuring the diame- 
ters of the circular lenses that are formed, a 
series of values of V and R in Eq. (28) are 
obtained. By inserting the values of ¢ from Eq. 
(7) or (27) and dividing by 7R?, this equation 
may be put in the form 


(41) 


where C,; and C, can be expressed in terms of the 
angles A, and A>»,,, the values of A; and A, fora 
lens of very large diameter. Experimentally we 
may thus obtain ¢,, by plotting the average lens 
thickness V/7R? against 1/R and extrapolating 
back to the intercept ¢,, at 1/R=0. This extra- 
polation may be greatly facilitated by knowledge 
of the theoretical values of C; and C, that can be 
obtained by the foregoing methods, from even 
approximate values of y1, y2 and yi2. From the 
values of ¢,, thus obtained, Fs may be calculated 
from Eq. (6). If y: and yz are known, yi2 can 
then be determined from Eq. (1), and in this way 
new values of A;, C; and C; can be obtained for 
more accurate calculations of ¢,,. 
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The measurements of the diameters of lenses 

thus provide a method of determining the 
interfacial surface tension yi2. The method 
becomes increasingly accurate as F's approaches 
zero. Instead of calculating A: by Eq. (40) from 
surface tension data, it is very easy to make 
direct accurate measurements of A; and A: by an 
optical goniometer sighted at the junction J of 
large lenses. Of course such measurements by 
themselves give B; by Eq. (34) and thus by Eq. 
(30) give yi2 when y; and y2 are known, but when 
Fs is small, still greater accuracy can be obtained 
by combining these measurements with determi- 
nations of the lens diameters. 
- It is also a simple matter to measure directly 
the linear tension f,, for a boundary of zero 
curvature and thus determine the relation of ¢ 
and ¢,. This may be done by inserting two 
parallel vertically placed metal strips into a lens 
on water. One of these strips is attached to the 
vertical pointer of a balance (a balance to meas- 
ure horizontal forces). The second strip may then 
be moved away from the first one so that the 
lens becomes so distorted as to become approxi- 
mately rectangular with parallel sides. The force 
exerted on the other strip as indicated by the 
balance, is thus 2f,.. 

If such measurements of f,, are combined with 
direct measurements of ¢, ¢,, and thus Fs may be 
determined by Eqs. (7) and (6) without needing 
to know the angles A. By using two sharp points 
attached to micrometers arranged so that one 
point can be brought in contact with the upper 
and the other with the lower surface, accurate 
measurements of ¢ can be made. 

These various methods thus provide us with 
accurate determinations of yi2. They have the 
advantage over the drop weight method usually 
used, that they are purely static methods and 
that they permit us to study the effects of known 
amounts of adsorbed substances at the interface 
by the methods which will be described in a later 
section. 


EXPERIMENTAL DATA ON THE DIAMETERS 
oF LENSES 
Weighed amounts of Squibb’s Petrolatum were 
placed on a large tray (45X75 cm) full of 
distilled water at 24°C, which was then covered 
by a sheet of glass to eliminate air currents. The 
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diameters of the lenses were measured to 0.1 mm. 
The values of V/zR?, plotted as a function of 
1/R, are indicated in Fig. 4 by small circles. 

The family of curves in this figure represent the 
calculated relation of V/2R? to 1/R according to 
Eq. (28) based upon the various assumed values 
of yi2 which are given on the curves. It is seen 
that the experimental points, except for lenses of 
small diameter, agree with the curve for y12 = 53.0 
dynes cm~. The value of ¢,,=0.4806 cm, given 
by the intercept of this curve on the V/7R? axis, 
corresponds to F's = — 11.37. 

Harkins® found for this oil at 20°C, by the drop 
weight method, ye=31.12; y2=55.32, Fs= 
— 13.69. ‘This high value of yi2 is incompatible 
with our measurements on lenses. To find 
whether impurities were present in the oil which 
could lower yi2, a portion was shaken with 
Fuller’s earth and centrifuged, but the lens 
diameter remained unchanged. We have checked 
Harkins’ value of ye by Wilhelmy’s method 


obtaining y2=30.8 at 23.5°C, which should FE 


correspond to 31.1 at 20°C. 

The method of calculating the curves in Fig. 4 
may be illustrated by the curve for yi2=53.0. 
Taking p2=0.8849 (measured at 20°C), p; 
=0.9982, yi=72.75 we find by Eq. (1), Fs 
= — 11.37 and by Eq. (6) t,, = 0.4806. By Eq. (10) 
we get @;=0.3949; a2=0.2681; ai2=0.9768. Eqs. 
(32) and (33) give B,=0.3894; B;=0.7039 
radians. By Eq. (40) for R= we find A: 
=0.1682 and Eqs. (34) and (35) give A:= 
— 0.5337; A12=0.5576. By Eq. (16) we obtain 
fi=0.29; fo=0.81; fiz=5.44, and by Eq. (15) 
f=6.54. Substituting these values into Eq. (27), 
taking the square root and expanding as a series 
in powers of 1/R we have 


t=0.4806(1+0.288/R+0.127/R2). 
From this, by Eqs. (40), (38) and (39) we find 
A; =0.1676(1+0.0063/R) 


and similar expressions for Az and by 
(34) and (35). Finally Eq. (28) gives 


=0.4806 —0.4029/R+0.0602/R’, 


which is the equation of the uppermost curve | ~ 


Fig. 4. 


Some measurements were also made of lens | 
of purified tetradecane on water at 25°C. The 4 
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Fic. 4. Calculated curves giving V/7R* (in cm) as a 
function of 1/R (in cm") for lenses of petrolatum on water 
at 20°C taking y2=31.12 and the values of y:2 marked on 
the curves. The dashed line is calculated for tetradecane 
on water at 25°C with y2= 26.9, yi2=51.3. The experimental 
data for petrolatum are marked by circles and for tetra- 
decane by crosses. 


tetradecane from Eastman was found to spread 
actively on water indicating the presence of 
hydrophilic substances equivalent to about one 
percent of stearic acid. Purification was tried by 
shaking with fuming sulphuric acid and washing 
with water. About 6-8 treatments gave a sample 
that did not spread. Agitation with Fuller’s 
earth, and removal of the latter by centrifuging, 
was found to be far more convenient and 
effective. The purified material behaves markedly 
different when placed on water: it does not spread 
and a large lens shows extraordinary tendencies 
to break up into small lenses, for example, when 
touched by a clean wet stirring rod. The points 
marked by crosses in Fig. 4 give the results of 
these measurements. 

The surface tension 2 of tetradecane at 23.5°C 
by Wilhelmy’s method was found to be 27.0 
Corresponding to 26.9 at 25°C. By choosing 
various values for y12 and calculating curves for 
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V/7R?, it was found that best agreement with the 
data was obtained with yi2=51.3. This curve is 
shown as a dashed line in Fig. 4. The data for 
tetradecane on water at 25°C are: p;=0.9970; 
p2=0.7615; ¥1=72.0; 26.9; yi2=51.3; Fs 
= —6.2; a,=0.3840; a2= 0.2690; = 0.667; B, 
= 0.2846; B; = 0.5626; A; = 0.0908; Ao = — 0.4718; 
A\2=0.3754; f:=0.080; f2=0.53; 1.66; f 
= 2.30; Ai=0.0908 (1—0.004/R); t=0.265 (1 
+0.187/R+0.053/R?); V/rR?=0.265 —0.147/R 
+0.0141/R?. 


EFFECT OF ADSORBED HYDROPHILIC MOLECULES 
AT THE INTERFACE OF A FREE LENS 


It is possible by the following technique to 
introduce known numbers of molecules of certain 
hydrophilic compounds into the interface be- 
tween the lens and the water. A large lens (about 
10 to 20 cm in diameter) is placed on the surface 
of water in a wide tray equipped with a horizontal 
balance to measure surface forces.’ The movable 
barrier is pushed back and a definite amount of a 
benzol solution of a spreadable substance, for 
example, stearic acid, is placed on the water, 
taking care that the amount added is insufficient 
to cover the available water surface. By ad- 
vancing the barrier, the monomolecular stearic 
acid film can be compressed and pushed into the 
lens while the force F can be measured by the 
balance. If an equilibrium state is reached, 
measurements of the total area (determined by 
the barrier), the area of the lens, and the surface 
concentration of the film on the water give data 
for calculating the number of molecules intro- 
duced into the lens. 

Instead of attempting to obtain equilibrium 
between the hydrophilic molecules on the water 
surface and those in the interface, it is usually 
better after introducing a desired number of 
molecules into the interface to remove the excess 
of hydrophilic material from the water surface (so 
that y1=7yw). The escape of these molecules from 
the lens is usually such a slow process that it does 
not cause appreciable loss during measurements 
of the lens diameter. 

The molecules thus introduced into the lens, if 
their number is not too great, distribute them- 


8]. Langmuir, J. Am. Chem. Soc. 39, 1848-1906 (1917). 
See especially p. 1869. 
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selves uniformly over the interface under the 
lens, and lower the interfacial tension ie 
(without changing 2) and thus cause a change in 
the radius of the lens. From a family of curves 
such as that shown in Fig. 4, it is then possible 
from measurements of R and V to determine 712. 

The lowering of the surface tension of pure 
water by the presence of adsorbed molecules may 
be usefully regarded as being due to a spreading 
force F of the two-dimensional gaseous or liquid 
adsorbed film. Similarly at the interface between 
water, w, and a hydrocarbon oil, R, we may put 


=Ywr — (42) 


where yur represents the interfacial surface 
tension between pure water and pure hydro- 
carbon when no hydrophilic substances are 
present in the interface. 

In this way the two-dimensional equation of 
state of adsorbed molecules at a_ water-oil 
interface can be studied in much the same way as 
has previously been done at the air-water surface. 

If the applied surface pressure is sufficiently 
increased, the molecules can be crowded into the 
lens in such numbers that some of them leave the 
interface and become dissolved in the hydro- 
carbon phase. The distribution of the molecules 
between the hydrocarbon phase and the interface 
for changes in y12 should thus be in accord with 
Gibbs’ equation. In a later paper it is planned to 
report on these investigations. 

Let us consider a hydrocarbon lens which is 
surrounded by water having an uncontaminated 
surface, but at whose interface with water there 
is a known surface concentration o12 of adsorbed 
hydrophilic molecules. The spreading coefficient 
Fs, as given by Eq. (1) is then 


Fs =Yu—YR—YRwt Fiz. (43) 


Let us denote by F, the spreading coefficient 
for the lens on water when hydrophilic substances 
are absent so that 


Fo=Yw—YR—YRw- (44) 
Thus from Eqs. (43) and (44) 
Fs=Fot+ Fi. (45) 


We have seen that for a paraffin oil, lenses are 
formed which are about 5 mm thick. As oi: is 


increased by adding hydrophilic molecules to the 
interface so that Fj. increases, the thickness of 
the lens, as given by Eqs. (6) and (45), decreases 
and must approach zero when Fj. becomes 
nearly equal to — Fo. Thus at a certain critical 
concentration of adsorbed molecules, free lenses 
can no longer exist, for the oil spreads out as a 
film over the whole of the available water surface 
(unless it thereby becomes monomolecular in 
thickness). 


CONFINED OIL LENSES 


Let us consider a large free oil lens with an 
interface containing adsorbed hydrophilic mole- 
cules upon an area of clean water in a tray, this 
area being limited by a movable barrier and a 
piston. By advancing the barrier the free water 
surface can be made to disappear so that the lens 
covers the whole available area. The lens may 
then be said to be a confined lens. If the area is 
still further decreased, a force F (dynes cm™') is 
exerted against the piston. By the same method 
as was used in the derivation of Eq. (6) and by 
using Eq. (45) we see that 


F= (46) 


By this method of the confined lens it is thus 


possible to determine Fy: as a function of o12.To 


avoid having the oil creep up onto the piston, it 
has been found useful to construct a boat-like 
piston with its edge slightly raised above the 
water surface. This method for the study of the 
equation of state of adsorbed molecules at an 
interface has the following advantages over the 
free lens method: 

(1) It is not limited (as is the free lens method) 
to the determination of values of Fi_ which are 
less than — Fp. 

(2) It permits o12 to be easily and reversibly 
changed by the movement of a barrier instead of 
requiring separate experiments involving the 
production of lenses of different diameters oF 
containing different numbers of hydrophilic 
molecules. 

(3) Because it possesses one more degree 0! 
freedom than the free lens method, it permits a? 
independent variation of the film thickness / and 
so allows any effects dependent on ¢ to be 
measured. 
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DupLex Firms Fs>0 


If the thickness of the oil film used in the 
confined lens method is less than about 0.1 
mm, the last term in Eq. (46) becomes less than 
about 0.005 dyne cm~ and is therefore negligible. 
The gravitational effect of film thickness thus 
disappears. We shall call such films duplex films. 
For these films we then have by Eq. (46) 


F= Fo+ Fiz. (47) 


Since F becomes zero when Fi.< — Fo, we see 
that by the study of duplex films Fi. can be 
determined only when — Fo. 

The concept of duplex films must evidently be 
restricted to cases where we are justified in 
considering that the film has an upper surface 
and a lower surface (interface) which are distinct 
from one another. However, in view of the small 
range of molecular forces, it is to be expected 
that this condition will be fulfilled until the film 
becomes at least approximately monomolecular. 
In a later section we shall discuss the conditions 
under which monomolecular films may be usefully 
regarded as being duplex films. 


EQUATION OF STATE OF ADSORBED ACTIVE 
GROUPS AT THE INTERFACE 


The active groups or heads of the molecules of 
hydrophilic substances present at low concen- 
tration in a liquid hydrocarbon on water tend to 
occupy positions in the interface. The hydro- 
carbon chains of these hydrophilic molecules, 
however, being in an environment of other 
hydrocarbon chains, have no tendency of their 
own to go to the interface; they are free to move 
except that one end of the chain must remain 
attached to the active group at the interface. The 
thermal agitation of the heads should tend to 
cause them to become uniformly distributed over 
the interface. This tendency must manifest itself 
as a spreading force Fiz which could be measured 
by any agent, such as a piston, which interferes 
With this indefinite spreading. For sufficiently 
low surface concentration o12 of the heads in the 
interface, this force Fiz should be given by the 
ideal two-dimensional gas law 


Fi o12k (48) 


In studies of monatomic caesium films on 


tungsten? it has been shown that the finite area 
occupied by each atom or molecule can be taken 
into account by modifying this equation to 


Fy2=o12kT /(1—o12/om), (49) 


where oy is the maximum concentration corre- 
sponding to a close-packed film. If we rewrite 
this equation in terms of a, the area per molecule 
instead of the concentration ¢ we have 


1/o12 (50) 
and 
F\2(a —ao) =kT, (51) 


where do for our present purposes may be re- 
garded as an empirical constant. 

It should be noted that this equation is based 
upon the assumption that there are no attractive 
or repulsive forces between the heads except 
those that prevent any two from occupying the 
same place at the same time. 

Of course at higher concentrations (12) it is 
conceivable that attractive forces between the 
heads might produce large deviations from Eq. 
(51) and might even cause the heads to condense 
into a liquid two-dimensional phase within the 
interface. Let us therefore estimate the magni- 
tude of the forces between the heads and see if 
they are large enough to cause deviations from 
Eq. (51). 

The dipole moment of the carboxyl group on an 
aliphatic chain is'!® about e.s.u. The 
vertical component of the dipole moment per 
molecule for films of fatty acids on water has been 
found" to be about 0.2 e.s.u. 

Since the force between electrically charged 
bodies in a medium of dielectric constant ¢ is only 
1/eth as great as in free space, the forces between 
dipoles immersed in water are only 1/80th as 
great as those given by applying Coulomb’s law. 
This decrease is, of course, due to the orientation 
of the water molecules, each of which has a 
dipole moment 1.8X10~'%. It is the resulting 


®], Langmuir, J. Am. Chem. Soc. 54, 2798-2832 (1932). 
See Eq. (59). 

10 Smyth and Rogers, J. Am. Chem. Soc. 52, 1824 (1930). 

1N. K. Adam and J. B. Harding, Proc. Roy. Soc. A138, 
411 (1932). 

12 J. H. Schulman and A. H. Hughes, Proc. Roy. Soc. 
A138, 430 (1932). 
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large attractive force between the water mole- 
cules and the carboxyl group that causes the 
spreading of films of fatty acids on water. 
Similarly the decrease in the observed dipole 
moment of the carboxyl group from 1.4 to 0.2 
x 10-'§ when brought into contact with water is 
to be looked upon as the effect of the orientation 
of water molecules, rather than being due 
principally to an orientation of the carboxyl 
group itself with respect to the surface as has 
been suggested by Adam and Harding. 

In connection with studies of caesium films on 
tungsten a general equation has been derived for 
the spreading force F’ resulting from the dipole 
forces between adsorbed molecules. The equa- 
tion there given (Eq. (81)) relates to fictitious 
dipoles resulting from adsorbed ions and their 
induced images. For the real dipoles of oil films 
we should multiply the coefficients by 2 and 
divide by ¢ the dielectric constant so that the 
equation becomes 


F’ = (2/6) (3.340°/2M2+1.53 


X10-%o?T (52) 


where J<0.89. If we put o=3.3 
X10" corresponding to a=30A*, T=300 and 
«= 80, we obtain F’<0.7 dyne as the effect 
of the spreading force in modifying Fi. given by 
Eq. (51) as a result of the dipole attractions or 
repulsions between the heads. The value of Fiz 
for a=30 and ao=10 is 20 dynes cm~.. Con- 
sidering that the dipoles are probably not so 
oriented as to produce the maximum attractive 
or repulsive forces, it appears reasonably certain 
that the dipole forces in oil films have little 
effect in modifying Eq. (51) until the heads of the 
molecules are practically in contact. 


EQUATION OF STATE FOR DUPLEX FILMs 


If we ignore the distinction between the upper 
and lower surfaces of a duplex film and consider 
the spreading force F for the whole film as given 
directly by the observed force on the piston, we 
obtain by combining Eqs. (47) and (51) 


(F— Fo)(a—do) =kT. (53) 


Taking F and a as variables, this is the 
equation of a rectangular hyperbola with asymp- 
totes at F= Fy and at a=dp. It should be noted 
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that since negative values of F cannot be 
observed by the balance method, the portion of 
the hyperbola below the a-axis cannot be ob- 
tained experimentally in the case of duplex films. 
If, however, we are willing to add a sufficient 
thickness of pure hydrocarbon to the film to 
obtain the gravitational effects characteristic of 
confined lenses, we may determine values of F), 
less than — Fy by use of Eq. (46) and from these 
by Eq. (47) can calculate the negative values of F 
that we cannot observe directly. 


EXPANDED Fi_mMs or Myristic Acip As Ex- 
AMPLES OF DUPLEX FILMS 


Fig. 5:contains a family of F-a curves for 
myristic acid films on 0.01N HCl solution ob- 
tained by Adam and Jessop." The ordinates 
represent the spreading force F in dynes cm‘, 
and the abscissas are the areas per molecule in 
A’. The figures marked on each of the curves give 
the temperatures at which the measurements 
were made. 

Each curve shows a kink which corresponds to 
the point J in Fig. 1. The curves to the right of 
these points J are almost exactly rectangular 
hyperbolas. The constants Fo and do in Eq. (53) 
have been chosen separately for each of these 
curves to obtain as good fit as possible. The 
values of Fy so obtained are found to be practi- 
cally the same over the whole temperature range, 
the average value being 


—11.2 dynes (54) 


The values of ao are found to vary with the 
temperature in accord with the equation 


ao=12+0.178t, (55) 


where ¢ is the temperature in degrees Centigrade. 

Introducing the values of Fo and ap into Eq. 
(53), F has been calculated as a function of a for 
each temperature, and the resulting values have 
been plotted in Fig. 6. In view of the fact that 
only three adjustable parameters are involved in 
Eqs. (53), (54) and (55), the agreement between 
the curves in Figs. 5 and 6 is remarkably good. 

The largest discrepancies between the two 
families of curves are those for the curve at 34.4° 


13N, K. Adam and G. Jessop, Proc. Roy. Soc. All2, 
362 (1926). 
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20 30 40 50 


Fic. 5. Experimental F-a curves for myristic acid films 
on 0.01 N HCl. 


at values of F greater than 20; these may well be 
the result of experimental error for these high 
temperature measurements are difficult. For very 
low values of F, F< 2, the slopes of the curves for 
the observed data are definitely less than those of 
the calculated curves. With these exceptions the 
root-mean-square deviation of the observed and 
calculated curves is only 0.21 dyne em~. 

A striking feature of these curves is that the 
F-a curves cut the horizontal axis at a rather 
definite angle indicating that any horizontal 
asymptote must lie at negative values of F. 

On the whole the agreement seems to be 
sufficiently good to justify the belief that these 
expanded films of myristic acid are essentially 
duplex films. The thickness of the films varies 
from 14.6A at a=30 to 8.8A at a=50. These 
thicknesses are far greater than the range of 
molecular forces exerted by nonpolar molecules. 
Since the films are liquid and of a thickness much 
less than the length of the extended molecule, 
22A, the molecules must be inclined at an average 


Fic. 6. Calculated F-a curves for myristic acid films. 


angle of about 50° from the vertical when a= 30 
and 66° when a=50. The molecules are probably 
just as free to move as they are in the interior of 
any liquid except that one end of each molecule 
is constrained to remain at the interphase. It 
seems reasonable to assume that this constraint 
is a matter of minor importance which would 
lead, at most, to an alteration in Fo, but would 
not change the character of the phenomenon. 
This hypothesis that a liquid monomolecular 
film of molecules containing long hydrocarbon 
chains is essentially a duplex film having distinct 
upper and lower surfaces may be tested in many 
different ways. In the first place the value of Fo 
given by Eq. (54) is in reasonable accord with 
the value —5.7 which we found for tetradecane. 
Another test consists of adding to the myristic 
acid before it is mixed with benzol and spread on 
water, various proportions of a pure hydrocarbon 
such as tetradecane. If we are right in considering 
the film to be a duplex film, the shapes of the 
F-a curves should remain practically unchanged 


| 
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except for a small modification due to a slight 
alteration of Fo. Some preliminary experiments 
to test this point have shown that the addition of 
an equal weight of tetradecane to myristic acid 
produces no appreciable change in the F-a 
curves. Some variation in Fp upon the addition of 
tetradecane may, however, result from the fact 
that the CH; group at both ends of the tetradecane 
molecule can occupy positions on the upper 
surface, while in the myristic acid molecule there 
is only one CH; that could take a position in the 
upper surface. This should have the effect of 
reducing y, in Eq. (43). This conclusion is based 
on several lines of evidence that the field of force 
around a CH; group is less than that for a CHe 
radical. For example, the surface tensions of the 
isohydrocarbons is less than the m-hydrocarbons; 
the cleavage of crystals of stearic acid, etc., lies 
along the planes separating the CH; groups 
rather than those separating the CH: groups. 


Two PHASE DupPpLEXx FILMS 


Adam’s curves for myristic acid in Fig. 5 show 
a sharp kink at a definite area a, and at a force 
F,;, these points corresponding to the point J in 
Fig. 1. 

Such an abrupt change in the direction of a 
curve must be due to the appearance of a new 
element in the system such as a new phase. 

By measurements of the electric surface po- 
tentials Schulman and Hughes* have found some 
evidence that the film is often observably 
heterogeneous for areas less than ay, while for 
areas less than dp (see Fig. 1) the film is again 
homogeneous. 

With two phases of constant composition in 
equilibrium with one another, F should remain 
constant. It is seen, however, from Fig. 5, that 
when a is-less than a, each curve has a slope 
which is definitely not zero. At sufficiently low 
temperatures, for high values of F, the myristic 
acid forms condensed films giving curves such 
as HGB in Fig. 1. It would seem probable, 
therefore, that the new phase which appears at a 
is a phase of this type consisting of close-packed 
molecules. A difficulty with this hypothesis, 
however, is that this new phase should then be 
one of constant composition with definite 


* Reference 12, p. 443, 


properties. This would require that the curve RJ 
should be a horizontal line. 

Adam has demonstrated beyond any reason- 
able doubt that in condensed films of substances 
containing long hydrocarbon chains which give 
curves like HA-in Fig. 1, the chains are close- 
packed just as they are in crystals of stearic acid. 
Myristic acid, however, at low temperatures and 
especially on pure water gives a break in the 
curve at G so that the F-a characteristic follows 
the line GB. Adam regards the area at the point B 
as being that of the heads. 

The strong spreading tendencies of the fatty 
acids prove that the carboxyl group attaches 


itself very firmly to water. In the preceding 


section it was seen that the agreement of the 
data with Eq. (53) indicates that the forces 
acting between the carboxyl groups in the ex- 
panded films are very weak. This must mean that 
some of the water molecules become so oriented 
so as to neutralize the moment of the carboxyl 
group. Very strong forces must then exist 
between these carboxyl groups and the water 
molecules. From this point of view as the film is 
compressed to an area less than az, the decrease 
in area is due to the progressive squeezing out of 
these water molecules until at the point G all the 
water has been removed from the spaces between 
the carboxyl groups. A suggestion of this kind 
has also been made by Schulman and Hughes.* 

The tendency of fatty acid films to form con- 
densed films increases as the chain length in- 
creases. We may assume that the energy involved 
per molecule, when the chains become close- 
packed, as at the point G, is proportional to the 
length of the chain. This energy is a kind of heat 
of fusion and should be of the same order of 
magnitude as that involved in the solidification 
of such substances as stearic acid. 

With these facts in mind we may now under- 
stand the nature of the curve JR. At the point J 
the chains of the fatty acid molecules in the film 
begin to assume the close-packed arrangement in 
which the molecules have a cross-sectional area of 
20.5A2. The heads of the molecules, however, 
because of their hydrated carboxyl groups, 
cannot space themselves on the surface of the 
water closer than 25A?. The solidification process, 


* Reference 12, p. 449. 
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therefore, can only proceed to a small extent by 
forming groups of molecules or micelles of the 
type represented diagrammatically in Fig. 7. 


B='5 


Elevation Cross section 


Fic. 7. Diagrammatic representation of a micelle. 


The crowding apart of the lower ends of the 
molecules cannot prevent the upper parts of a 
few adjacent molecules from arranging them- 
selves as in condensed films, but no large number 
of molecules can gather together in this way 
because of the crowding apart of the heads. The 
new phase that appears at the point J, therefore, 
consists of micelles containing relatively few 
molecules. Because the molecular weight of these 
micelles is not indefinitely large, they therefore 
exert a pressure which adds itself to that of the 
phase of single molecules. Strictly speaking, the 
film between R and J is not a two-phase system, 
but is a mixture of two two-dimensional gases, 
one consisting of normal molecules and the other 
consisting of micelles, each of which contains B 
molecules. The kink at J is due to the sudden 
appearance of the micelles at this point. On the 
basis of these hypotheses we may readily calcu- 
late the shape of the curve that we should expect 
between R and J. 


THE SPREADING ForcE Fz DUE TO THE MICELLES 


Consider an interface having a total area S 
between a hydrocarbon liquid and water and let n 
be the total number of molecules in the interface 
which act as molecules in a two-dimensional gas, 
then the area per molecule is given by 


a=S/n. (56) 


Now let us consider that the ” molecules in the 
surface exist partly as single molecules and are 
Partly gathered into micelles, each of which 
contains 6-molecules. Let us assume further that 
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an equilibrium exists between the micelles and 
the single molecules such that the partial 
spreading force F, of the single molecules remains 
constant as long as the micelles are present. Let 
S,; be the total area at which the micelles first 
appear as the area is decreased. We may thus put 


ay=Sy/n, (57) 


where ay is the area per molecule at the point J 
(Fig. 1) when the whole surface is covered with 
single molecules giving a spreading force F;. Let 
Sr be the area of the film when it is compressed to 
such a point that the single molecules disappear 
so that all the molecules are present as micelles. 
The number of micelles is then 2/8, and each of 
these covers the area BSp/n=Bar. Let us now 
consider the intermediate stages where Se<S 
<S,. Let m; be the number of single molecules 
and mg the number of micelles. We then have the 
equation, 

mi +Bys=n. (58) 
We may derive a second equation involving the 
areas covered by the molecules and by the 
micelles. Since the partial pressure F; remains 
constant, each of the m, single molecules requires 
the same area a, as it does when the film area is 
Sz, whereas each micelle covers the area Bap. 
The total surface areas, according to Eq. (56), 
gives an, and this must be equal to the area 
required by the single molecules and the micelles. 
We thus have 


(59) 


Let o; and og be the surface concentrations of 
the single molecules and of the micelles. We then 


have 


op=ng/S; (60) 


Combining these equations with Eqs. (58), (59) 
and (56) and then dividing by S we obtain 


o1+Bos=1/a, 
= 1. 


(61) 
(62) 
Solving these equations for o; and og we find 


(63) 
(64) 
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The partial pressure Fs of the micelles in ac- 
cordance with Eq. (49) may be written 


(65) 


where agp is the concentration of micelles at the 
point R. Since the number of micelles is then 1/8 
we have by Eq. (56) 


opr =n/BSr=1/Bar. (66) 


Substituting the values of org and og from Eqs. 
(66) and (63) into Eq. (65) we obtain 


(67) 


Thus between J and R the curves should be given 
by 


F=Fj+Fs, (68) 


where Fs is given by Eq. (67). 

Examination of the curves in Fig. 5 shows that 
they are of the type required by this theory. 
There is a sharp kink at J as the theory would 
indicate but there is no kink at R where the 
curve joins the line GB, the approach to this line 
being asymptotic because of the denominator in 
the second member of Eq. (67). 

The values of 8 have been calculated from 
Adam’s data of Fig. 5. For this purpose the 
values of ag were assumed to lie along a straight 
line like GB in Fig. 1 with a value of ag of 25.A?. 
Thus 


ar=25—0.156F.. (69) 


The values of 8 obtained in this way from the 
experimental data on F and Fy; by Eqs. (67) and 
(68) varied from about 9 to 18. For values of a 
only slightly greater than ap the accuracy of these 
determinations of 6 is naturally not large for it 
depends too much on the rather arbitrarily 
assumed value of ap. 

If we may assume that the molecules tend to 
arrange themselves in a hexagonal lattice like a 
honeycomb, each molecule should be surrounded 
by six others and the next ring of molecules 
should also contain 6. The number of molecules in 
a micelle may thus naturally be assumed to be 7 
or 13. The value 6=13 is approximately the 
average of the values found from the data of 
Fig. 5. 


Two-DIMENSIONAL VAPOR PRESSURE OF 
THE MICELLES 


The points J in Fig. 5 thus correspond to the 
two-dimensional vapor pressures of the micelles, 
that is, to the value Fi, the spreading force due to 
the single molecules in equilibrium with micelles. 
If we examine Fig. 5 we see that the points 
corresponding to J decrease steadily as the 
temperature is lowered and become zero at a 
definite temperature of about 6°. This is a very 
unreasonable behavior for a vapor pressure. 
When we consider, however, that the true vapor 
pressure is exerted at the interface, we see that 
the observed spreading force F does not measure 
this pressure. The value of F; is thus given by 


F,\=Fy—Fo. (70) 


With the data for F, in Fig. 5 and by using the 
value of Fy from Eq. (54) the logarithms of 
(F;— Fo) have been plotted against the recipro- 
cals of the absolute temperature. As is usual for 
vapor pressure curves, a straight line was 
obtained; this corresponded to the equation 


logio (Fy — Fo) =8.058—1950/T. 


The fact that these vapor pressures measured 
from Fo) agree with the Clapeyron equation, 
whereas the uncorrected values do not, affords 
further justification for regarding these myristic 
acid films as duplex films. 

The values of F; calculated from Eq. (71) have 
been used to fix the points J given in Fig. 6 
which are marked by small circles. Comparing 
Figs. 6 and 5, it will be seen that the agreement is 
probably nearly within the experimental error. 

With these values of Fy and ap, as given by 
Eq. (69) and by taking 6 = 13, the curves in Fig. 6 
to the left of points J have been calculated by 
Eqs. (67) and (68). The general agreement 
between the slopes of these portions of the curves 
and those given by the experiments indicates 
that the value of 8 is approximately constant and 
equal to 13. The reasonableness of this value of 8 
thus helps justify the hypothesis that in this 
two-phase system the condensed phase consists 
of micelles. 

Comparing the curves of Figs. 6 and 5 we see 
that there are considerable discrepancies 4 
temperatures of 5 and 2.5° and in general at very 
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TABLE I. Expanded films of fatty acids on dilute HCl at F=1.4 dynes cm~. 


1 2 3 4 5 6 7 8 
No. of 
Catoms to ty — Fy) 
Acid n a4 < (obs.) do — Fo calc.) (calc.) 
Lauric 12 47.6 14.7 <0°C 14.6 10.4 10.8 —9.0 
Tridecylic 13 44.6 13.5 <0 14.4 11.4 12.1 +0.2 
Myristic 14 43.7 21.0 9 15.7 13.0 13.4 9.6 
Pentadecylic 15 42.1 23.9 19.5 16.2 16.2 14.7 18.8 
Palmitic 16 41.5 34.8 28.5 18.1 16.5 16.0 28.2 
Margaric 17 42.2 42.2 37.5 19.4 17.5 17.3 37.4 
Stearic 18 43.0 52.8 46 21.3 19,2 18.6 46.6 
Heneicosoic 21 46.7 78.0 65.5 25.8 21.6 22.5 74.2 
Behenic 22 46.0 75.0 72.5 25.3 21.6 23.8 83.2 


low values of F. The probable reason for these 
deviations will be discussed later. 


ApAM’s DATA FOR OTHER Fatty AcIps 


For the other saturated fatty acids Adam has 
not published full curves such as those of Fig. 5 
for myristic acid, but gives merely the area a;.4 of 
the expanded films at a standard small com- 
pression of F=1.4 dynes cm and at some 
temperature f distinctly above the range at 
which expansion at this compression occurs. He 
also gives the temperature ¢; at which half the 
expansion from the condensed to the expanded 
film has occurred at F=1.4. The data of the first 
four columns in Table I are those of Adam." The 
value of a1,4 for myristic acid is about 9 percent 
lower than the value given in his later papers, as 
shown in Fig. 5, so that exact comparison with 
our foregoing analysis is not possible. 

Taking do as given by Eq. (55) and as tabu- 
lated in column 5, we may calculate Fy by means 
of Eq. (53) from Adam’s data on a and T after 
placing F=1.4. These values, which are given in 
the 6th column, when plotted against m, lie 
roughly along a straight line corresponding to 


Fy=4.8—1.3n. (72) 


The 7th column gives the values of — Fo 
calculated by this equation. It is seen that — Fo 
increases 1.3 dynes cm= for each additional 
carbon atom. This agrees approximately with the 
increment of 1.5 which Harkins® found for — Fs 
between hexane and octane. 

The temperatures f; as given in the 4th column 
may be used to generalize Eq. (71) for all 


“N. K. Adam, Proc. Rov. Soc. A101, 521 (1922). 


saturated normal fatty acids. Referring to Fig. 5, 
the point marked by a cross just below the curve 
for 7.2° represents half expansion at F= 1.4 for it 
lies about halfway between the intersections with 
the curves corresponding respectively to JM 


- and GB in Fig. 1. These data thus give about 7.0 


for the temperature of half-expansion. The value 
of F,; at this temperature is evidently about 0.4 
lower than F=1.4, so that the spreading force F; 
of the micelles, as defined by Eq. (70) is to be 
taken as 1— 

Eq. (71) may be written 


logio F,=8.058—B,/T, (73) 


where B, is a constant which depends upon 1, the 
number of carbon atoms in the fatty acid 
molecule. Taking the values of — Fy from the 7th 
column of Table I, adding 1 to get — Fi, substi- 
tuting this in Eq. (73) together with the value of 
T corresponding to ¢, enables us to determine B,, 
for each acid. These values plotted against m give 
roughly a straight line corresponding to 


B,,=1208.+53.n. (74) 


To show how well these values agree with the 
experimental data, they have been used in Eq. 
(73) to calculate (from Fy in column 7) the 
temperatures of half-expansion which are given 
in the 8th column. The agreement with the data 
of the 4th column is good except for the last two 
acids where, as Adam points out, the experi- 
mental accuracy is much less. The values of B, in 
Eq. (74) correspond to a heat of “evaporation” 
for the micelles of 5.52+0.242n kg cal. per 
mole or (3.82+1.67”)10-“ erg per molecule. 

The large value for the constant term in these 
expressions (5.52 kg cal.) suggests that more than 
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half of the energy of formation of the micelles is 
due to attractive forces between the carboxyl 
groups when they are in contact. The increment 


_ of 0.242 kg cal. per mole for each CH: should 


correspond roughly with the heat of fusion of the 
hydrocarbons which is 0.60 kg cal. per CHe.'** 
Referring to Fig. 7 we see that when 13 hexagons 
are closely packed, 30 out of the original 78 
hexagonal faces are still free. The contribution to 
the heat of formation of the micelle should thus 
only be 48/78ths of 0.60 or 0.37 kg cal. per CHa, 
even if the solidification were as complete as in 
a large crystal. The experimental increment 0.242 
which is thus 2/3 of the heat of fusion helps 
justify our hypothesis of the cause of the micelle 
formation. 

It is possible by the methods we have given for 
myristic acid, applied to the data contained in 
Eqs. (55), (72), (73) and (74), to construct the 
whole family of F-a curves for any saturated 
normal fatty acid. When experimental data are 
available, it will be interesting to compare them 
with these predicted curves. 

For unsaturated fatty acids on dilute HCl, 
Adam" gives the data contained in the first four 
columns of Table II. For oleic acid on dilute 


TABLE II. Unsaturated fatty acids. 


1 2 3 + 5 6 7 
Acid nas ty a B, 


Iso-oleic 18 46 35°C 25 18.6 24.9 2016 
Oleic(trans) 18 52.5 10 <O (11.2) 21.7 (1700) 
Elaidic(cis) 18 47 10 —0.5 14.0 (21.7) 1875 
Erucic(trans) 22 47 10 45 164 25.2 1892 
Brassidic(cis) 22 45 40 34. 19.2 22.6 2073 


HCI at 15° Adam and Jessop have given a 
complete curve!® which agrees within a standard 
deviation (S.D.) of 0.24 dyne cm~ with Eq. (53) 
if we take a9=19.9 and Fy=—11.2. For the 
corresponding saturated acid (stearic) at this 
temperature, we should have, by Eqs. (55) and 
(72), a9 =14.7 and Fy= — 18.6. Thus the effect of 
introducing the double bond in the middle of the 
hydrocarbon chain has been to increase do by 
5.2 and to decrease — Fy by 7.2. The increase in 


4s W. E. Garner, K. Van Bibber and A. M. King, 
J. Chem. Soc. (1931), 1533. 

18N. K. Adam and G. Jessop, Proc. Roy. Soc. A120, 
473 (1928). 
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do indicates that the double bond tends to 
occupy a position on the water surface in accord 
with the known fact that the effect of a double 
bond is to increase the solubility of a hydrocarbon 
derivative in water. 

Although for the other unsaturated acids, since 
only a;,4 is given, we cannot determine dp» and F, 
separately, we may draw certain tentative 
conclusions. 

Taking Fy for oleic acid to be —11.2 from the 
more complete later data, we calculate, by Eq. 
(53), 21.7.The difference in for oleic and 
elaidic acids is probably due to a change in Fj 
rather than to one in do. We thus take a)= 21.7 
for elaidic acid and find the value of Fo given in 
the 5th column. For the corresponding acids 
with 22 carbon atoms, assuming the increment of 
1.3 per CHe as in Eq. (72), we take a value of 
—F, 5.2 units greater than for the 18 carbon 
acids. For iso-oleic acid with the double bond 
between the a, B-atoms, Fy has been taken to be 
the same as for stearic acid. The variations in the 
resulting values of ad) are probably of no signifi- 
cance, since the data for 1.4 are uncertain by 
about 2 units (Adam), but that they are on the 
average 7.5+2.5 larger than the values given for 
the saturated acids by Eq. (55) indicates that the 
general effect of the double bond is to raise ap. 

The values of B, in the 7th column have been 
calculated from t; and 1— Fp by Eq. (73). They 
should be compared with the values of 2162 for 
stearic and 2374 for behenic acid given by Eq. 
(74) for n=18 and 22, respectively. The effect of 
the double bond near the carboxyl group in iso- 
oleic acid is to lower B, by 146 which is equiva- 
lent to a decrease in of 2.8. A double bond near 
the middle of the chain lowers B by an average 0! 


294 (equivalent to a decrease of 5.6 in m) in the | 


case of the cis-compounds elaidic and brassidic 
acids, while in the trans-compound erucic acid 
the decrease is 482 which corresponds to the 
effect of shortening the chain by 9.1 CHe groups. 
These very large and specific effects of the 
double bonds near the middles of the chains are it 
striking contrast with the small changes produce 
on the JK part of the curve as derived from the | 


data for a;,4. This latter part depends only on Fi a 


and do, whereas the value of Fy; at any temper 
ture (which determines B) involves all the 


factors that determine the stability of th . 
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micelles. Any irregularity in the chain naturally 
affects the ease with which the molecules can 
pack into a micelle. 


COMPOUNDS WITH MULTIPLE CHAINS 


In condensed films of such substances as 
tristearin, glycol dipalmitin, etc., which have 
two or more chains of equal lengths, each chain 
occupies the area 20.5A? so that the area for a 
whole molecule increases in proportion to the 
number of chains. With expanded films, however, 
the area depends only on dp and Fp and thus, 
according to Eq. (53), should not be expected to 
increase in proportion to the number of chains 
even if a) does so. Adam finds that the area at no 
compression is 135 for triolein at 15°C, while for 
tristearin and tripalmitin at 60°C it is 115. He 
notes that these areas are about 2.4 instead of 
three times the areas of the “corresponding 
compounds with single chains.”’ 

Adam" has given a complete curve for glycol 
dilaurate at 15°C which agrees within about 
SD=0.2 dyne cm=! with Eq. (53) if we take 
a)=39 and Fy= —8.5. Comparing this with the 
data for lauric acid in Table I we see that the 
value of ap has been increased very greatly, from 
15 to 39, while — Fy has decreased by 2 units. 
Adam and Jessop’s curve for triolein" gives best 
agreement with Eq. (53) (rather poor) with 
ao=90 and Fy= —8.5. 


THE Errect oF UNEVEN UPPER SURFACE 


Films prepared from mixtures of cholesterol 
and myristic acid'® on dilute HCI give areas less 
than the sum of those obtained from the con- 
stituents separately. 

Cholesterol consists of large compact molecules 
that form a condensed film with an area of 40.8 
and a thickness of 15.34. An expanded myristic 
acid film of an area a=47 has a thickness of 9.3. 
In the mixed film the large molecules thus extend 
through the film and tend to produce humps in 
the upper surface. If the total area of the upper 
surface is increased, for example, 10 percent, 
from this cause, and the surface tension of the 
upper surface is 50 dynes cm“, the effect is to 
increase —Fy by 5 dynes cm. This may 


“N. K. Adam, Proc. Roy. Soc. A126, 366 (1930). 
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perhaps be the main cause of the decrease in area 
of the expanded films containing cholesterol. 
Similar effects were observed when tripalmitin 
and other large molecules were used instead 
of cholesterol. The three chains in tripalmitin 
being attached at one end must lie closer than if 
unattached and produce a local thickening of the 
film. 

The micelles present in the JR portion of the 
curves of expanded films should also tend to 
cause an uneven upper surface. This effect 
should be especially marked when Fy, is small, 
that is, at low temperatures, for then the film of 
single molecules is so thin that the micelles rise 
above it. It appears probable that the deviations 
between the curves in Figs. 5 and 6 at low values 
of a and at low temperatures are due to this 
cause. 


THE EFFECT OF LARGE HEAD GROUPS IN DE- 
CREASING 8 


For para-dodecyl-phenol a family of curves 
have been obtained" like that for myristic acid. 
The area at no compression is, however, unusu- 
ally small (a=36). The JK part of the curves at 
30°C gives ap= 16.3, Fo= —20. This large value 
of — Fo explains the small value of a at F=0. 
The JR portions are far steeper than those of 
myristic acid and are found to correspond to 
micelles for which B=5+1. This low value is in 
accord with the presence of the phenyl group in 
the head which should restrict the number of 
chains that can unite in a micelle. The value of 
B, is 1970. 

A similar example of the effect of a large head 
group in reducing the value of 8 is furnished by 
the data'® on the a-bromo-fatty acids. For the 
acid containing 17 carbon atoms, at 15° the JR 
portion of the curve gives a value 8=7. The JK 
portion gives @)=17.1; Fo= —13.7. 


FILMs OF ESTERS AND VAPOR EXPANDED FILMS 


A family of curves is also given for mono- 
palmitin.'” The JK portions agree well with Eq. 
(53) (SD=0.3) when ao=22 and Fy= —8.8. The 
area at no compression, in accord with these 


17N. K. Adam, W. A. Berry and H. A. Turner, Proc. 
Roy. Soc. A117, 532 (1928). 
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values, is large (a=67). The values of Fy— Fo 
furnish, by Eq. (74), the value B, = 2052. When 
the values of Fy are then calculated by this 
equation at the four temperatures (from 16 to 
37°) at which the experiments were made, the 
calculated and observed values agree with an 
S.D. of 0.13 dyne cm , although only one ad- 
justable parameter is available in these calcu- 
lations. 

The JR portions of the curve are fairly flat. 
Calculation shows that B=13.54+2. The magni- 
tude of 8 seems to depend very largely on 
whether or not the F—-a curve for the condensed 
film possesses a GB branch (Fig. 1). We have 
concluded that the line GB is due to the hydration 
of the heads. If there is little or no hydration, the 
F-a curve follows the line GA. The fatty acids 
and particularly the a-bromo acids show a very 
distinct GB line with values of ag of about 26 
which indicate that the carboxyl group has a 
strong tendency to be hydrated and then 
occupies an area distinctly greater than 20.5. 
With the alcohols the point B lies at about a= 22. 
Thus the larger size of the head in the mono- 
palmitin as compared to the acids is offset by a 
lower degree of hydration. 

The effect of replacing the hydrogen atom in 
the carboxy] group of the fatty acids by a methyl 
group is to Jower the expansion temperature ¢; by 
one or two degrees, while an ethyl group produces 
a lowering of 15°. This indicates that the 
tendency to form micelles is decreased very 
greatly by the presence of an ethyl or higher 
group in the ester chain. The condensed films of 
these esters show, however, an area of 20.5 so 
that the fact that the—COO—group must be 
pulled away from the water in the formation of 
the micelles (as in the condensed films) explains 
the lesser stability of the micelles. 

Another interesting feature of these ester films 
is that the GB section of the curve is wholly 
absent in the condensed films and as a result of 
this lack of hydration the values of 8 are very 
large (the JR lines are flat). This means that the 
expanded film passes to the condensed state more 
nearly as a definite change of phase instead of by 
the growth of micelles. The curves for propyl 
palmitate give a)=12, Fy = —4.8, B,= 2000, and 
8=60. For ethyl palmitate’ similar data give 
ao=12, Fo= —5.7, B, = 2059 and B=55. 
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The small values of — Fo for these esters are 
presumably due to a lack of sufficient flexibility 
of the molecule in the neighborhood of the—COO 
—group to allow both chains to assume a steep 
angle to the water surface. With glycol dilaurate 
and triolein the larger value of — Fyo=8.5 results 
from the additional linkages between the two 
chains by which greater flexibility can be 
attained. 

The marked effect of attachment to the water 
in the middle of a chain in decreasing — Fp is best 
illustrated by the data which Adam and Jessop 
have given" for films of oleic acid on acidulated 
water containing potassium permanganate. The 
oxidation of the double bond gives strong 
attachment to the water near the middle of the 
chains. The curve at 15° yields a=27.1 and 
— Fy)=2.2, which is the lowest value derivable 
from any of Adam’s curves. 

The films which Adam has called vapor 
expanded films are evidently merely those which 
are characterized by very low values of — Fp. At 
low compressions the curve then extends to such 
large values of a that the film becomes so thin 
that the molecules lie nearly flat on the surface 
before F falls to zero. Thus the films can no 
longer be considered duplex films and it is not 
surprising that they should merge into true gase- 
ous films without the abrupt change of slope that 
characterizes typical (liquid) expanded films. 


FACTORS WHICH DETERMINE do 


In the foregoing analysis of the experimental 
data a» has been treated as an empirical constant, 
although its derivation from Eq. (49) would 
suggest that ao should have the value of a which 
corresponds to a close-packed film. The empirical 
values of ad are usually less than this. 

The fact that no two molecules in a mono- 
molecular film can be in the same place at a given 
time indicates the existence of repulsive forces of 
very short range. Such forces yield a term 
which should correspond to the value of a for 4 
condensed film (at least for fairly concentrated 
films).° 

We have seen, however, particularly in our 
discussion of Eq. (74) that in the formation of the 
micelles there is evidence of short range attractive 
forces between carboxyl groups. Such forces 
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should tend to counteract the effect of the 
repulsive forces at contact and thus tend to 
decrease the value of ado. We should therefore not 
attach too much significance to the actual values 
of ay but may well interpret the differences in do 
for different substances as indications not only of 
the relative sizes of the active groups, but also of 
the magnitude of the forces between them. 


FLEXIBILITY OF HyDROCARBON CHAINS 


Adam has summarized" the evidence that the 
chains are flexible. The fact that the upper and 
lower surfaces of duplex films of a thickness of 
only 5-10A? can act nearly independently just as 
in a thicker film of liquid is perhaps still more 
convincing evidence. 

In a series of papers Staudinger" and his 
coworkers have developed a theory of the 
viscosity of solutions of high polymeric com- 
pounds which is based on the assumption of long 
rigid molecules. It would seem that such theories 
require revision. Perhaps it is the shearing stress 
in a liquid during viscosity measurements that 
pulls a very long molecule into a nearly straight 
line and thus, in part, justifies the theoretical 
derivation of the equations. 


NoTE ON ANTONOFF’S RULE 


Antonoff’s rule!® states that for true equilib- 
rium Fs=0. As Harkins and others have pointed 
out, this cannot be true in general. For example, 
for a pair of liquids A and B there are two 
spreading coefficients: for lenses of A on B, 
Fss=ys—ya—vYap, while for lenses of B on A, 
Both spreading coefficients 
can be zero only when yas=0, i.e., when the 
liquids are miscible. 

It is nevertheless a fact that Fs is extremely 
close to zero for very large numbers of organic 
liquids placed on water and in equilibrium with 
it, so that as a useful empirical rule Antonoff’s 
rule is justified. 

If the molecules of a vapor of an organic liquid, 
when they strike a water surface, form an 
adsorbed film which nearly completely covers the 
surface, then some molecules must condense on 

“* H. Staudinger and E. Ochiai, Zeits. f. physik. Chemie 
A158, 35 (1931). 


* Antonoff, J. Chim. Phys. 5, 372 (1907). Bartell, Case 
and Brown, J. Am. Chem. Soc. 55, 2769 (1933). 


top of those already adsorbed. We then have two 
cases to consider: 

(1) The molecules in the 2nd layer may have a 
shorter life ¢ (before evaporating again) than the 
life t) of molecules on the surface of the organic 
liquid. 

(2) The life ¢ of the molecules may be equal to 
or greater than fo. 

In the ist case the 2nd layer of molecules will 
not become complete so that F's <0, while in the 
2nd case the 2nd layer and therefore presumably 
a 3rd and 4th layer, etc., will grow in presence of 
saturated vapor until the water is covered by a 
liquid film and therefore Fs=0. 

According to the principle of independent 
surface action” which the writer has so often 
employed, the properties of a surface depend 
primarily on the atoms or groups of atoms 
actually forming the surface and depend on 
the underlying layers only insofar as they alter 
the actual surface atoms. Thus in general the 
properties of the upper surface of an adsorbed 
monomolecular film are nearly the same as those 
of the substance in bulk. This may be illustrated 
by some examples. 

Carbon tetrachloride vapor forms an adsorbed 
monomolecular film on water which tends to 
become nearly complete before the vapor be- 
comes saturated. The properties of the upper 
surface of the film thus approximate to that of 
liquid carbon tetrachloride, so that the life ¢ of 
molecules of CCl, which may condense on the 
film is about the same as that of molecules on 
liquid CCl,. Thus the film grows in thickness 
which means that Fs=0. 

Orientation of the molecules in the adsorbed 
film may make ¢ for the second layer very 
different from to. Thus lenses of molten palmitic 
acid on water, although in equilibrium with it, do 
not spread indefinitely. This is due to the fact 
that the orientation and tight packing of the 
molecules in the film on the water has the effect 
of preventing the carboxyl groups of any vapor 
molecules incident on the surface from coming 
into contact with other carboxyl groups. On the 
other hand a molecule striking a surface of liquid 
palmitic acid can more easily penetrate the less 
tightly packed surface layer and thus its carboxyl 

207, Langmuir, Chem. Rev. 6, 467 (1929); also Colloid 
Symposium Monographs 3, 48 (1925). 
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can meet other carboxyls. Thus ¢ for molecules 
in the 2nd layer of the adsorbed film is less than 
to and therefore Fs <0. 

A corollary of this point of view is the con- 
clusion that if, for two immiscible liquids A and 
B, A is strongly adsorbed (from a saturated 
vapor) on B then the vapor of B cannot be 
strongly adsorbed on A (since both Fsa and Fszp 
cannot be zero). 


The frequency with which strong adsorption 
occurs and the relatively small effects on the lives 
t produced by orientation is thus the fundamental 
reason that Antonoff’s rule is so often applicable. 

Of course there are numerous cases such as 
that of the higher hydrocarbons and water where 
only a negligible adsorption of either vapor by 
the other liquid occurs. Such cases are usually 
those yielding large negative values for Fs. 
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Molecular Symmetry and the Reduction of the Secular Equation 
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Henry EyrinG, ARTHUR A. Frost AND JOHN TuRKEVICH, Frick Chemical Laboratory, Princeton University 


A method is given for expressing a bond eigenfunction 
of any multiplicity in terms of a linear independent set. 
The independent sets for the multiplicities corresponding 
to one bond and to two bonds for an arbitrary number of 
electrons are given. A method recently given for obtaining 
matrix components for singlet bond eigenfunctions is 
generalized to include all multiplicities. Group theory is 


(Received September 15, 1933) 


applied to the bond eigenfunctions for symmetrical 
molecules to reduce the secular equation and examples 
are worked out for all multiplicities. The possible reasons 
for the approximate additivity of bond energies for 
molecules is examined in connection with the chemist’s 
custom of ignoring all but the strongest bonds. 


INTRODUCTION 


N the calculation of activation energies for 
reactions involving complicated molecules it 
is important to know at least approximately how 
nearly the experimental bond strength corre- 


sponds to a simple exchange integral between the | 


paired electrons. In the many reactions for which 
this is approximately the case the electrons not 
actually changing partners may be ignored or 
considered as causing only small steric effects. 
The binding energy for a few typical molecules is 
given here. We hope later to extend the cases 
examined to include more complicated molecules 
and to supplement these considerations by calcu- 
lations of the integrals involved. 


REDUCTION OF BoND EIGENFUNCTIONS TO A 
LINEARLY INDEPENDENT SET 


Rumer' has shown how a linearly independent 
set of singlet bond eigenfunctions may be chosen. 
Points representing the atomic eigenfunctions are 
drawn as a polygon without reentrant angles. 
Straight lines representing bonds are then drawn. 
Only bond eigenfunctions for which no bonds 
cross are included in the linearly independent set. 
To express any singlet eigenfunction in terms of 
this linear independent set we require but one 
relationship. This may be expressed diagram- 
matically as follows: 


*G. Rumer, Nach. d. Ges. d. Wiss zu Gottingen, M. P. 
Klasse, p, 337 (1932). 
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This relationship holds not only for bond 
eigenfunctions involving four electronic eigen- 
functions but for those involving an arbitrary 
number providing the bonding not indicated in 
our diagram is identical for all three bond 
eigenfunctions. This same relationship holds even 
when an arbitrary number of the electronic 
eigenfunctions are connected by no bonds, i.e., 
for all multiplicities having as many as two 
bonds. For multiplets higher than singlets one 
additional relationship is required. If m bond 
eigenfunctions are identical except for a single 
bond in each and if these » unequal bonds added 
together vectorially have a resultant equal to 
zero then the sum of the bond eigenfunctions is 
zero. Diagrammatically this is simply 
+ + f- @ 
These two relationships enable us to express any 
bond eigenfunction in terms of a linearly inde- 
pendent set. Both propositions are readily veri- 
fied by expressing the bond eigenfunctions in 
terms of Slater functions. A recent article by 
Rumer, Teller and Weyl'* is of interest in this 
connection. 

If our 2m electronic eigenfunctions (m an 
integer) be thought of as arranged in the polygon 
described above, a linearly independent set of one 
bond eigenfunctions is formed by 2n—1 of the 2” 
eigenfunctions around the perimeter which then 


requires but one bond to form a closed figure. 


1 Rumer, Teller and Weyl, Nach. d. Ges. d. Wiss. zu 
Géttingen, M. P. Klasse, p. 449 (1932). 
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Since any other bond will form a closed figure it is 
clear all other bond eigenfunctions can be 
expressed in terms of the set selected, and as 
there are (2n—1) independent eigenfunctions we 
have a linearly independent set. An equally 
satisfactory set corresponds to the 2n—1 bonds 
connecting a particular electronic eigenfunction 
with each of the others. 

A linear independent set of two bond eigen- 
functions corresponds to the ”(2n—3) ways that 
two bonds can be drawn on the perimeter, never 
of course using the same electronic eigenfunction 
in both bonds of the pair. Our two propositions 
enable us to express any other two bond eigen- 
function in terms of this set. The independent 
set for more than two bonds will involve bonds 
not falling on the perimeter and is, therefore, 
slightly more difficult to describe. The two 
propositions given are, however, all that are 
required to single out a satisfactory set for any 
number of bonds. 

Kimball and Eyring? defined an operator which 
acting on all the antisymmetric Slater functions 
for which M, (the 2 component of spin) is zero 
gives the particular bond eigenfunction with 
M,=0. To form a bond function with M,+0 
proceed exactly as before. Take the linear combi- 
nation of all Slater functions, with the proper M,, 
which associate opposite spin eigenfunctions 
with each pair of electronic eigenfunctions bound 
together. Give the Slater function the positive 
sign if it associates negative spins an even 
number of times (or not at all) with the first 
letters in each bond as written, otherwise give it a 
negative sign. The proof that the functions so 
formed are eigenfunctions for the operator S’, the 
total angular momentum, proceeds exactly as in 
the case M,=0.* 

The secular equations for a particular inde- 
pendent set of bond eigenfunctions is the same 
for each value of M, except for a constant factor 
multiplying every term in the equations. This 
factor is smallest for the largest possible M, for 
each multiplet. This factor of course divides out 
of the secular equations and we get the same 
eigenvalue repeated in a multiplet for each M,. 


2G. E. Kimball and H. Eyring, J. Am. Chem. Soc. 54, 
3876 (1932). 

3H. Eyring and G. E, Kimball, J. Chem, Phys. 1, 239 
(1933). 
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CALCULATION OF MATRIX COMPONENTS BeE- 
TWEEN BOND EIGENFUNCTIONS OF 
ALL MULTIPLICITIES 


Bond eigenfunctions with a different number 
of bonds belong to different multiplets and are 
noncombining for the energy operator considered 
in which magnetic terms are neglected. Consider 
two eigenfunctions with the same number of 


bonds 


Now to calculate the matrix component write 
down ‘the electronic eigenfunctions in order 


abcdefghigkIimn 
Bi Bi a1 Bi aj a; a as Be as Bs 


and below each letter an a or 6 is written. We 
shall call this array of a’s and #’s a sequence. 
Assign an a or 8 to an electronic eigenfunction in 
the following way. If there is a bond in the 
first eigenfunction between two letters neither of 
which form part of a bond in the second eigen- 
function write ay under one letter and 8; under 
the other. The letters 77 are an example. For a 
similar bond in the second eigenfunction, as &/, 
do the same except use the subscript s. If either 
eigenfunction has more than one such bond the 
matrix component is zero to the approximation 
of first order exchange integrals. If of two letters 
forming a bond in one of the eigenfunctions for 
example e—g only one, e, is bonded in the second 
eigenfunction write a; under the unbonded letter 
and £; under the bonded one. Now continue 
assigning successively a; and B; to the chain of 
letters connected to these by bonds (connected 
first in one function and then in the other) until 
you again come to one which forms part of a 
bond in one function but not in the other. Here 
the chain terminates. We shall use the word 
chain for such a series of atoms to distinguish 
them from a cycle defined below. The chain 
gef is an example. If there is another chain of 
this sort do the same using the subscript /: 
Continue until all chains have been assigned. | 
Now if a letter forms part of a bond in both | 
functions, for example a, write a; below it and 

write 6; below the letter to which it is bound, } 

in our example, and below the companion of } li 
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in the second function write a; alternating thus 
around the cycle until you return to the initial 
letter. Such a cycle in the example is abcd. A 
cycle may contain as few as two letters or as 
many as all of them. Now if there is a second 
cycle proceed exactly in the same way using the 
subscript 2 and so on until all the cycles have 
been assigned. Finally below letters which are 
not bonded in either eigenfunction write a 
without any subscript. We next write a + or — 
sign before our sequence of a’s and f’s depending 
on whether to the totality of first letters of bonds 
in the two functions an even or odd number of 
6's have been assigned. In this example four 8's 
have been assigned so the sign is +. In general if 
6's have been assigned p times this sign is 
(—1)?. 

From our sequence we may now tell immedi- 
ately how many Slater eigenfunctions are com- 
mon to the two eigenfunctions. If the subscript f 
or s occurs obviously there are none. If these 
subscripts do not occur there are 2” where m is 
the largest subscript used for numbering the 
cycles. This is because for each cycle there was 
really an arbitrary choice of a or £8 to be assigned 
to the first letter while outside of cycles there was 
no choice. The coefficient of the coulombic 
integral is zero if the subscript f or s occurs 
otherwise it is (—1)?2". The corresponding 
coefficient of the energy W in the secular 
equation is always just the negative of the 
coefficient of Q. 


Case I 


If either the subscript f or s occurs on an @ as 
many as two times in a sequence the matrix com- 
ponent is identically zero. 


Case II 


If both f and s occur as subscripts on but 
one bond only the four exchange integrals occur 
which may be formed by taking one of the 
letters from each of these two bonds. The two 
integrals corresponding to letters to which like 
spins have been assigned occur with a plus sign 
the other two with minus signs. It is to be 
understood in addition that these quantities are 
multiplied by the factor (—1)2”. Thus for the 
two eigenfunctions written down the matrix 
component is 2(ik+jl—kj—il) since (—1)?2"=2. 
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Case III 

Next consider the case where either the 
subscript f or s does not appear while the other 
subscript appears only on a single bond. If the 
bond k—/ in eigenfunction y2 is replaced by g—h 
we have such an example. The sequence would 
then differ from the former one only in having 8; 
assigned to h while to k and / are assigned a’s 
without subscripts. The value of (—1)?2™ is 
unchanged. For this case the only terms ap- 
pearing will be exchange integrals between one of 
the two letters in this bond (when the subscript f 
has been assigned) and a second letter from a 
chain. All such integrals where like spins have 
been assigned in the sequence enter with the 
coefficient 1 while the others have the coefficients 
—1. All are multiplied by the factor (—1)?2” as 
before. For the two eigenfunctions just described 
the matrix component is 


—ie—th—jf —jg). 
Case IV 


Finally consider the case where neither the 
subscript f or s appears. An example is obtained if 
in eigenfunction Y2 we replace the bond k—/ by 
4—j. The sequence would then differ from the one 
written down in that to both k and / would be 
assigned a’s with no subscripts and the subscripts 
associated with 7 and j would be changed from f 
to 2. Again it is to be understood that the 
coefficients given below are all to be multiplied by 
(—2)?2™. The coefficients of the coulombic 
integral Q and all exchange integrals between 
letters to which opposite spins with identical 
subscripts are assigned is then 1; integrals 
between letters with like spins and like subscripts 
have coefficients —2; those between letters with 
positive spin and no subscripts have the coeffi- 
cient —1; those between letters in cycles and any 
letter outside the cycle have the coefficient —3; 
any integral between a letter in one chain and a 
second letter outside this chain but not in a 
cycle with spins alike has the coefficient —1, if 
the spins are unlike the coefficient is zero. 

This gives the general scheme for all multi- 
plicities. For singlets the sequence contains only 
cycles and the rules reduce to the very simple 
form already given.*: 4:5 The proof of the rules 
given follow from the sequence by inspection 
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using considerations of the type found in the 
article*® just referred to. 

We now consider the simplifications arising in 
secular equations from molecular symmetry. 


GENERAL GROUP THEORETICAL CONSIDERATIONS 


We shall state the few group theory propo- 
sitions required without proof and refer the 
reader to standard treatises for more detail.®: 7: § 
If we have a molecule containing like atoms 
symmetrically placed certain bonds are equal and 
the matrix components in the secular equation 
are entirely unchanged by any rotation or 
interchange of atoms which simply changes like 
bonds into like bonds. These operations can be 
expressed as a group of permutations which 
acting on any bond eigenfunction changes it into 
another bond eigenfunction which may or may 
not be a member of our linear independent set. 
It can in any case be readily expressed in terms of 
such a set. The effect of such a permutation on 
our linear independent set of bond eigenfunctions 
can be expressed as a matrix and the sum of 
the terms on the principal diagonal, called its 
character, is denoted by x;. Now the group of 
operations breaks up into classes such that all 
members of a class C;, have the same character 
xi. We denote by 4; the number of members in 
this class. If all the members of a class C; be 
multiplied into the members of a class C; we get a 
relationship of the type C;C;= jC; (3) where 
Cije iS some integer or zero. The identity, c 
always constitutes a class by itself; so that x; is 
the degree of our representation, being simply 
the number of functions operated upon. For a 
particular group there is only a limited number 
of irreducible representations possible, i.e., repre- 
sentations for which no linear combinations of 
the functions can be found such that they break 
into two or more sets not transforming into each 
other for any of the group operations. These 
may be found by using the relationship h;h;x:x; 
= x124C;;x4xxx (4) in connection with the set of re- 


4H. Eyring and G. E. Kimball, J. Chem. Phys. 1, 626 
(1933). 

5L. Pauling, J. Chem. Phys. 1, 280 (1933). 

6 A. Speiser, Theorie der Gruppen, 2nd edition. 

7E. Wigner, Gruppentheorie. 

8H. Weyl, The Theory of Groups and Quantum Me- 
chanics, tr. by H. P. Robertson. 
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lationships (3).°:* The possible choices of the 
characters in (4) which are consistent with (3) 
fix the characters of the irreducible representa- 
tions. Examples will be given later. 

Now in molecular problems for each multi- 
plicity we have a linear independent set of bond 
eigenfunctions as a representation of the group 
corresponding to the molecular symmetry. This 
representation, unless especially chosen, is in 
general reducible for a symmetrical molecule. By 
taking suitable linear combinations we may break 
up our representation into equivalent irreducible 
representations of the group any one of which 
may fail to appear or may appear once or more 
times. There is a simple way of determining how 
often each irreducible representation will appear, 
however. If for each class we multiply the 
character, x;’, of our reducible representation by 
the character of the irreducible representation in 
question, x;, and multiply this by the number of 
members of the class, h;, the sum of all such 
quantities divided by the number of mem- 
bers in the group, N, gives the number of 
times, , the irreducible representation will 
appear 1/N2Zjhix;’x;=n (5). One other theorem is 
frequently useful. The average value over all 
members of the group of the square of the 
character of the irreducible representation equals 
the sum of the squares, G, of the number of 
times each irreducible representation will appear 
in the reduced representation. 1/N,hix;" =G (6). 

The advantage of using irreducible repres- 
entations is that they are non-combining and if a 
particular one of degree x; occurs but once all x: 
roots are equal; so that only a linear equation 
must be solved to get the energy. If a particular 
irreducible representation occurs / times we must 
solve an equation of the /-th degree. 

Having determined a particular irreducible 
representation is to appear, very simple means 
suffice to obtain it. Beginning with any bond 
eigenfunction from the reducible representation 
we apply an operation to it and add to it the 
eigenfunctions which give the proper character 
for this class. Should this turn out to be impos- 
sible this eigenfunction does not occur in this 
representation and we start with a new one. 


Having found eigenfunctions which transform 


®H. Bethe, Ann. d. Physik 3, 133, 19 (1929). 
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properly for a member of one class we adjust it to 
have the proper character for a member of a 
second class adding functions which do not spoil 
the first transformation. This is repeated for a 
member of each class. In some cases it was found 
easier to start with the most general repre- 
sentation which had the proper character for the 
first class and specialize it for successive classes. 
For the identical representation we begin with 
any eigenfunction and add together all the 
quantities obtained by operating on it with each 
member of the group. The linear independence of 
two representations of the same irreducible 
representation (equivalent representations) is 
readily investigated by expressing them in terms 
of the linear independent set. 


AMMONIA 


We assume the molecule has one three-fold 
axis of symmetry and three planes of symmetry. 
We take the three electronic eigenfunctions on 
the nitrogen atom as b, d, f and assume they are 
bonded respectively to the three hydrogen 
eigenfunctions a, c and e. Our symmetry opera- 
tions are then C; : identity; C2 : (ac) (bd), (ae) (Of), 
(ce)(df); Cs : (ace)(bdf), (aec)(bfd). The relation- 
ships (3) are C2 =3C,:+3C;; and 
C:C;=2C2. These with Eq. (4) give for the 
characters of the classes of the three irreducible 
representations: 


Ci | C2 | Cs 
rj 
1 | -1] 1 
r; | 2 | 0 | -1 


We take as the linear independent set for the 
singlets the Rumer set of bond functions: a—f 
b-—cd—e, a—bc—de-f, a—fb—ec—d, a—d 
b-—ce—f and a—b c—f d—e. The group of 
Operations transform these five eigenfunctions 
among themselves with only one new eigen- 
function appearing d—a b—ec—f which by the 
method of resolving crossed bond functions is 
readily shown to be just the negative of the sum 
of the five Rumer functions. We write for the five 
Rumer eigenfunctions yi to Ws; respectively and 
¥s for the other eigenfunction. The members of 
the group are represented by the letters E, A, B, 
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C, D, F respectively. When the operations act on 
the eigenfunctions in the first line of the table we 
get the eigenfunctions in the line with the 
operation. 


E || | v2 | vs | | 
Ve | vo | vs | | vs 
v2 | vs | vs | 
We | vo | | | 
| v2 | va | | vs 
Foi | ws ws | 2 


Here the operations in the first column acting on 
the first row of eigenfunctions give the subse- 
quent rows with the character given in the last 
column. Using proposition (5) we find there are 
two irreducible representations of T; one of T's 
and one of I3. These are found almost by 
inspection to be respectively: YitwWe, ve; 
and ¥3;—¥ with ¥;—y, the latter of course being 
a two-dimensional representation. The five 
energies are, respectively, Q+a—$8—3y—6 
+[(2a— 26)*?+ (36 —$y —36)?]!; O—3a; and two 
equal energies with the value Q—a— 26. Here the 
interchange integrals ab=cd=ef=a; ac=ae=ce 
=B; bd =bf =df ad =af = bc = be =cf =de=5. 


> 


TRIPLET STATES 


Next consider the triplet terms. These are two 
bond eigenfunctions and the method of selecting 
the linear independent set of nine eigenfunctions 
has been described. In specifying the bond 
eigenfunctions only the bonded electronic eigen- 
functions are written down. Also to further 
simplify the notation a bond will be indicated by 
simply writing the two bonded eigenfunctions 
next to each other leaving a space between 
atoms not bound together. Thus ad cd indicates 
the bond eigenfunction with a bonded to } and 
with c bonded to d it being understood that e and 
f are unbound. A line will be used, however, for 
one bond eigenfunction to prevent confusing this 
notation with that for interchange integrals. The 
linear independent set is then ab cd, ab ef, cd ef, 
ab de, af bc, af cd, af de, be de, bc ef. By going 
through precisely the same sort of considerations 
as for singlets these eigenfunctions break up into 
two representations of ¢: =ab cd+ab ef+cd ef 
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and ¢.=af bc+af bd+ch de+ch df+ed fb+ed fa; 
one representation of ¢;=ab de+dc be+ab cf 
+dc af+ef ad+fecb; and three representations 
of T3, ¢s:=abcd+abfe with $¢;=abcd+cd fe, 
o6=ac bd+fb ae with ¢;=ac bd+ec df, and finally 
gs=ab cft+ab de+ef ch+ef da with ¢ be 
+dc af+ef ch+ef da. 

These energies lie considerably higher than the 
lowest singlet and so are not particularly inter- 
esting. The matrix components can be readily 
calculated by the scheme outlined. 


QUINTET STATES 


For the quintet states we can conveniently 
take for the five linearly independent single bond 
functions [a—b], [a—c], [a—d], [a—e] and 
[a—f]. Proceeding as before we get one represen- 
tation of T;, Q,=[a—b]+[c—d]+[e—f]; two 
representations of T3, Q.=[a—e]+[a—f]+[b 
—f]+[b-—e] with 
+[d-—e] and Q;=[e—f]+[b-—a] with 2;=[e 
—f]+{¢d-c]. 


SEPTET STATES 


These states correspond to no bonds and there 
is a single eigenfunction for each M, correspond- 
ing to this multiplicity, namely the sum of all 
Slater functions for each M,. For the largest M, 
this eigenfunction is simply the single Slater 
function with all spins positive and the energy is 
Q-—3(a+6+7+28). This is the highest possible 
energy for such a system if all the interchange 
integrals are negative. 

It is interesting to consider under what 
circumstances the lowest energy for the system 
corresponds to the energy which would be 
calculated from y2, the bond function chemists 
ordinarily assume. This happens when 6 = 3(8+7) 
and of course for the much more specialized case 
when these three integrals are zero. The results 
given for ammonia apply to all cases of three 
equivalent bonds, such as might be expected in 
C=N, unless in addition we have symmetry about 
the plane perpendicular to the triple bond as for 
example in the case of the nitrogen molecule. 


NITROGEN 


The additional symmetry in nitrogen over 
ammonia just doubles the number of group 


operations and also the number of classes; so 
that proceeding exactly as in the case of ammonia 
we find the characters of the six irreducible 
representations given in the following table. 


Ci | Co | Cz | Ca | Cs | Co 
Tr; 1 1 1 1 1 1 
1 1 1 | -1 
i j-ii 1 ~i 
i 1 1 | 1 
rT; 2 0 2 @ | 
Ts |} 2 0 | -—2; -1 1 
h; 1 3 1 2 3 2 


The last line of the table gives the number of 
operations in each class. The operations of 
classes Ci, C2, Cy are the same, respectively, as 
Ci, C2, C3 for ammonia. The other classes may be 
specified by one operation taken from each; 
C3(ab) (cd) (ef) ; Cs(ad) (bc) (ef) ; Ce(adebcf). Each of 
the irreducible representations of ammonia have 
split into two representations; the first one into 
the first pair and so on. Exactly the same 
eigenfunctions form irreducible representations 
for nitrogen as for ammonia but with the 
difference that certain eigenfunctions now belong 
to different representations and are non-com- 
bining. The energies are exactly as written for 
ammonia if we remember that now B=y. The 
following is the assignment of eigenfunctions to 
representations for nitrogen. For the singlets 
Ty, and Yo; Ts, We; Ys with Ys— 
Triplets $1, and $2; Ts, $3; T's, with $9; 
T¢, with ¢; and with ¢;. Quintets %; Ts, 
Q; with Q;; I's, Qe with Q3. 


METHANE 


The group of twenty-four symmetry opera- 
tions includes all the rotations of a regular 
tetrahedron into itself plus the interchange of 
two hydrogen atoms, the rest of the molecule 
remaining fixed. The method of procedure is pre- 
cisely the same as before. The five irreducible 
representations have the characters given in the 
table. 


i 
i 
a 
3 
n 
~ 
el 


REDUCTION OF THE 


C3 
‘Ti 1 1 1 1 1 
Ts 2 0 21-11 
r; 3 1 
hj 1 6 3 8 6 


The four hydrogen eigenfunctions are indicated 
by the letters 0, d, f, h and they are supposed 
bound, respectively, to the four carbon eigen- 
functions a, c, e, g. The five classes are identified 
by the five permutations given in order: identity; 
(eg) (fh) ; (ae) (bf) (cg) (dh) ; (ceg)(dfh) ; (aceg) (bd fh). 
The representations for the fourteen singlet 
states are then: three representations of Ti, yi 
=II and y2=HI+V+XIX and ¥;=IV+VI 
+VIII+-XI+XII+XIV; one of Tx, 


3a—(9/2)(B+y) -—3W/4 


—9a—18(8+7) +95—9W/2 
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two of T;, 
—III with and y,=XI+XII-IV 
—VI with one rep- 
resentation of Ty, ¥y=VI—IV, VII—XIV, 
¥i1= XII—XI; and finally one representation of 
yu=XVI 
— XVIII. The bond eigenfunctions indicated by 
the Roman numerals are: I=ahbedegf; II 
=abcd ef gh; Il=ahbgcfde; 1V=abchdg fe; 
V=ad be eh fg; Vi=af be cd hg; VI1=ah be dg ef; 
VIII =ah bg cd fe; 1X =ab ch de fg; X =af be de gh; 
XI =ab cd eh gf; X11 =ad bc ef hg; XII1 =ah becd 
fg; XIV=ab cf de hg; XV=<ah be cf dg; XVI 
=ad be chfg; XVIl=af bgchde; XVIII =af bc 
dg eh; XI1X=af bechdg; XX=ad bg cf he. The 
last six can be expressed in terms of the first 
fourteen but nothing is gained by this after linear 
independence of the eigenfunctions for each 
representation has been established. 
The lowest energy lies in the symmetrical 
eigenfunctions, [';, which involves the solution 
of the cubic: 


12a0—(27/2)(B-+7) 
9a— (63/2) (8-+y) +185—9W/2|= 0. 


36a—72(8+y7)—27W/2 


a is the energy of one of the four exchange 
integrals corresponding to principal bonds; 6 
corresponds to the interaction between hydrogen 
atoms; y to the interaction between two bond 
electrons on carbon and 6 to the interaction of 
an electron on hydrogen with one on carbon to 
which it is not bound. Examining the cubic it can 
be seen that if 6=(8+y)/2 one root of the 
equation, the lowest, is the upper principal 
diagonal element. This is the energy correspond- 
ing to ¥; the bond function regularly assigned to 
methane. This eigenfunction is, therefore, prob- 
ably a rather good approximation. The letter W 
as written includes only the interchange energy. 
The energy corresponding to Tz is Q—2a—68 
where Q is the coulombic binding. The four 
energies for T'; correspond to taking the plus and 
minus sign each twice in Q—(2/3a+5/3(6+y7) 
+46) £2/3[a?+ 
—3(8++7)6]!. For Ty we have a triple root with 
energy Q—28—2y—45 and for T's; we have a 
triple root Q—2(a+6+y+8). 
These results offer some justification for 


Woods" use of the lowest bond eigenfunction in 
an approximate calculation. For more accurate 
results a cubic should of course be used. 
Van Vleck’s!"' very interesting papers on methane 
treat the whole question from a somewhat 
different point of view. 

Instead of using the symmetry group on the 
bond functions one could first apply it to get 
representations in terms of Slater functions with 
subsequent separation of the multiplets. This 
was carried through for certain of the cases given 
here but there seemed to be no real advantage 
and the labor of calculating the final matrix 
components is greatly increased. The method as 
outlined seems readily applicable to all sym- 
metrical systems of bond functions. 

F. Seitz, A. Sherman and A. E. Stearn have 
considered this problem from a somewhat differ- 
ent point of view and we wish to thank them 
for helpful joint discussions. 

 H. J. Woods, Trans. Faraday Soc. 28, 877 (1932). 


J. H. Van Vleck, J. Chem. Phys. 1, 177 (1933); 1, 219 
(1933). 


a 
e 
of 
iS 
1; 
of 
0 
e 

| 
or 
ts 
95 
» 
- 
of 
le 
le 


NOVEMBER, 1933 


JOURNAL OF CHEMICAL PHYSICS 


VOLUME 1 


On the Quantum Mechanical Stability of a Benzol Molecule 
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The Heitler-London-Rumer treatment of valence bonds 
is applied to the benzol molecule, considered as a system 
of six trivalent nitrogen-like CH groups. Partial solution 
of the secular equation of the thirty-fourth degree shows 
that this model of the benzol ring is thermodynamically 


unstable relative to three acetylene molecules, in disagree- 
ment with experiment. The method does not take cogni- 
zance of the directional properties of valence bonds 
(Slater-Pauling). Editors. 


T the first attempt to approach the interpre- 
tation of the stability of a benzol molecule 
from the standpoint of Heitler-London’s homo- 
polar bending it is natural to start from the 
simplified model of a benzol molecule, namely, 
to assume that the benzol ring consists of six 
nitrogen-like atoms; that is, to take each of the 
CH groups that enter symmetrically into the 
structural formula of benzol as a single nitrogen- 
like trivalent atom. 

The indication of the nitrogen-likeness of the 
CH molecule may be found, for instance, in the 
work of R. S. Mulliken,! who gives for the 
reaction of formation of CH the following 
formula: 


C(1s?2s22p2) +H (1s) + D, 


where D is the heat of dissociation and 1s”, 2s?, 

++, etc., are the current notations in the 
systematics of spectra. Thus, we replace all CH 
in the benzol molecule by nitrogen atoms in 
their fundamental states. 

Following the usual method for solving such 
problems it is necessary to start from ‘‘eigen- 
value problem” in the space of valency functions: 
(A—W)¢=0. In case of a cyclic chain of six 
trivalent atoms 1, 2, 3, 4, 5, 6 the operator W 
is equal to Ty+T7o3+T34+Tast+Tsct+To1, and 
T., is a differential operator of the kind? 7,, 
=D,,D,:—3. In other words the ‘exchange 
integrals” between the neighboring atoms are 


1R. S. Mulliken, Phys. Rev. 40, 55; 41, 49 (1932). 
*M. Born, Ergebnisse der Exakten Naturwissenschaften 
Vol. 10, 1931, 
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put equal to 1 and the interaction between the 
noncontiguous atoms is neglected. 

Further, it is necessary to select independent 
valency functions ¢;; “to build a basis’’; i.e., to 
select a complete set of independent functions 
gi each of them being a product of symbols of 
the kind ([xy]; [xz] ---). The symbol ([xy]}) 
corresponds to a valency bond between the atoms 
x and y and, moreover, to a valency bond having 
a definite sense (from x to y); i.e., to a valency 
arrow. 

If one uses the methods developed by G. 
Rumer the building of the basis does not present 
any difficulty even in the case of such a compli- 
cated system as a benzol molecule. To build the 
basis one has according to Rumer’ to proceed as 
follows. Atoms are arranged arbitrarily along a 
circumference‘ and linked up by valency arrows; 
each of the figures of this kind having no inter- 
secting arrows then corresponds to an_ inde- 
pendent valency state, and all other figures with 
intersecting arrows (as shown for the general 
case by Rumer) represent superpositions of 
independent valency states. For example: 


If the brackets of the kind ([xy]) correspond 
to the valency arrows directed from X to Y, 


3G. Rumer: Zur Theorie der Spinvalenz, Nachrichten 
von der Gesellschaft der Wissenschaften zu Gottingen. 
1932. No. 27. 

‘ Generally along the convex side of the closed curve. 
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we can write the same in symbols: [xy ][zu]=[xu ][zy]+[«z][yu]. 
In our case there are thirty-four independent valency states. They are given in Fig. 1. 
Operating on the valency functions by the operator W, using the notations of the valency 
functions of Fig. 1, and expressing the result by independent valency states, we obtain: 
Way, = 20a, + 18b,4+ 31; 
We, = 16c; — 2i2 — 6d2— 6d5— 18m2—4a4, 
Wd, = — — 6b2+ 
Wb, — 8m2+4 2d3— 2d2+ — 2dg, 
W1, = 221; — 8d¢+8d3+ 2c; — 
Wk, = 22k; — 8d2+ 8d; + 2c; — 
= 30); —18a4+ 18a), 
Wi = 182; 92, 
Wm, =6m,+4b; C2. 


Thus, generally speaking, we get a matrix of the thirty-fourth order. But as the operator W 
commutes with the cyclic operator Z, there is the possibility of reducing the order of the matrix. 
By the cyclic operator we understand the operator which acting on any function ¢;(a;; bi; ¢; «+ -) 
represented in Fig. 1 changes it into if i<m (é=1, 2 m) and into +4 or — if i=n. 
Transforming the operator Z to the principal axes, we get 5 


0506) 
+0705) 


046) 


Denoting linear combinations of the kind a; 
by A;, C;, etc., we have in case of c; and d; also: 


ZQ,=0°C;, ZD,=0'D,, 
ZC2=0'C2, ZD2=0'Ds, 
ZC3=0°C3, ZD3=0°D3, 
ZC;=0C;, ZD;=aD;, 
ZCg=1Cs, 
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In case of 0;: 


= 


Denoting the corresponding combinations of 
the kind 0; by capital letters B; etc. for 1;; Ri; ji 
we get: 


ZJ\= —o?J;; 
ZJ3= —1J3; 


ZL,=—0°L; 
ZL2 = 
ZL3 =— 1L3; 


ZK,= 
ZK2=0K3; 
ZK3= —1K3; 
and finally: 

+42) =1(ti-+t2); Z(mi+me) =1(m+me) ; 
Z(t: —t2) = ; m2) = —1(m,— m2). 


Thus our matrix of the thirty-fourth order 
may be decomposed into six matrices: one of the 
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ninth, two of the seventh, two of the third, and 
one of the fifth order. Each of these matrices 
corresponds to a definite eigenvalue A; of the 
operator W. Acting with the operator W on the 
functions denoted by capital letters we obtain 
the following matrices: 


6-r» 0 0 1 0 
54 12 | 
24 863 0 4 
20-r» -8 0 
-4 6 (2) 
0 
-8 
6 +» 3) 
10-» 0 0 0 1 -2 
0 0 30-r’ 0 18 0 0 
12 0 0 22-r~ 3 2 |=0, 4) to 
36 6 2 6 20-r’ 8 0 
0 0 0 0 4 va 
: 
10-rA 0 0 0 1 0 —2 the 
12 22-v’ 0 0 3 2 8? 
Cas 
0 0 30-r’ O 18 0 0 es 
12 0 0 3 2 (5) pla 
. 36 6 2 6 2-rA 8 0 " 
0 0 0 0 be: 
-12 -4¢ 402 O -2 16-2 cate 
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0 
0 
1 
3 


When using the mirror symmetry K; and L; the order of the latter matrix may be lowered by a 


unit, and we obtain 


As the energy of three N,-like molecules is equal 
to —9 we are interested only in matrices whose 
roots lie in the region of negative \ and especially 
in the region lower than —9, for only these 
values for the energy secure mechanical and (in 
case A< —9) thermodynamical stability of the 
benzol molecule as well.® 

It is quite evident that all the matrices with 
the exception of the last one have no roots in 
the region A<—9. In fact, the proper values 
Ae +++ Of matrices in the general case (in 
case of complex matrix elements ajx) lie only 
within a certain closed region (in a complex 
plane) made up by circles K; (i=1, 2, «++ 1) 
with centers in a,; and radii: 


|ax|, 


‘The minus sign before the “exchange integral” indi- 
Cates attraction. 


where >’ means the sum over all k with excep- 
tion of k=.’ Transforming our matrices in a 
suitable way (by multiplying rows and columns 
by certain numbers) we can slightly narrow the 
region of proper values of these matrices deduced 
directly from 1, 2, 3, 4, 5, 6. 

The proper values of the transformed matrix of 
the fifth order lie between —6<d;<36, of the 
third between 6 <A; < 28, of the seventh between 
—8<);<40, of the eighth between —17<); 
<54. Matrices (2) and (3), (4) and (5) have the 
same roots. The root in question hence can arise 
only from the last matrix (6’), which when 
transformed may be written: 


7S. Gerschgorin: Uber die Abgrenzung der Eigenwerte 
einer Matrix, Bulletin de |l’Academie der Sciences de 
U.S. S. R., 1931. 
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d 6—-rv 0 8 0 0 0 1 0 
0 18-rA 0 0 0 0 9 0 
24 10-rA O 0 0 0 4 
12 0 12 22-A oO 0 2 
0 0 0 0 30-rA 0 18 0 0 |=0. (6) 
12 0 12 0 ae 2 8 
0 0 36 6 2 6 20-r 8 0 
1) 54. 0 0 0 0 16-r 12 
48 3 24 4 0 4 0 4 16-2 | 
2) 6-r» 0 8 0 0 0 1 0 
0 18-r 0 0 0 0 9 0 
24 0 10-vA 0 0 1 0 4 | 
24 0 24 22-r O 6 4 16 
=0 (6’) 
3) 0 0 0 0 30-r’ 18 0 0 
0 oO 36 6 2 2-r» 8 0 
54.6 0 0 0 16-r 12 
48 24 4 0 0 and \y=22. 
(4) 
(5) 
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The direct substitution of definite values of \; 
shows that the lowest proper value of this 
matrix lies between —8 and —9; hence in the 
model chosen by us the benzol molecule is stable 
mechanically but not thermodynamically. 

We may suppose that this thermodynamical 
instability is the result of the incompleteness of 
our model of the benzol molecule. We have 
assumed the nitrogen-like atoms replacing CH 
in our model of the benzol ring to be in the state 


4S; we have further assumed the benzol model to 
be flat while we had to refer to it as to a space 
one, in fact it is known that N does not enter 
NH, in its fundamental state and that NH; (as 
the result of the absence of spherical symmetry 
of N) has not a plane structure, but a pyramidal 


one. 

In conclusion I wish to offer my best thanks 
to Professor Rumer for the suggestion of the 
problem and for his continued help in this work. 
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A short discussion of the nature of distribution equilibria 
between mixed crystals and aqueous solutions is given. 
From the general theory of the thermodynamical potential 
and the activity theory of strong electrolytes, equations 
are deduced giving the influence of the composition of the 
liquid phase on the distribution ratio and the dependence 


of the partition constant on the properties of the pure 
components. The relative importance of the various factors 
in these equations, depending on the concentrations and 
properties of the components, is discussed. The deduced 
equations are in agreement with the experimental data. 


HE conception of solid solution was first 

introduced into science by J. H. van’t Hoff 
in the year 1890.! In support of such a view of 
solid mixtures he had already pointed out the 
distribution between solid and liquid solutions. 
If solid mixtures are really solutions the distribu- 
tion of a substance between a solid and a liquid 
phase must take place according to the Berthelot- 
Nernst distribution law—C,/C,=Kwy, where C, 
and C, are the concentrations of the micro- 
component in the solid and liquid phases, 
respectively, in case the following three require- 
ments are fulfilled: (1) the solutions must be 
dilute enough with respect to the distributing 
substance in order that the variations in its 
content do not affect the thermodynamic prop- 
erties of the solutions; (2) true equilibrium 
between phases and throughout each phase must 
be attained; and (3) the distributing substance 
must be in the same molecular state in both 
phases. 

Since then a great number of investigations 
have been published concerning the question, 
but with very discordant results.2 The cause of 
the discrepancies lie in the difficulty of attaining 
equilibrium, because of the extreme slowness of 
diffusion in crystals at ordinary temperatures; 
for the same reason inhomogeneous crystals once 


' J. H. van’t Hoff, Zeits. f. physik. Chemie 5, 322 (1890). 
*G,. Bruni, Feste Lésungen und Isomorphismus, Leipzig, 
1908, where the early literature is given. See also G. 
Tammann, Zeits. f. anorg. Chemie 87, 1 (1919); G. 
Tammann and A. Svorykin, Zeits. f. anorg. Chemie 173, 
i G. Meyer, Rec. Trav. Chim. Pays-Bas 42, 304 
23). 


formed ‘‘freeze in’’ and do not spontaneously 
become uniform. Many investigators have paid 
much attention to this circumstance and have 
even considered the equilibrium conception, as 
well as the thermodynamical laws, inapplicable 
to the crystalline mixtures.* 

But thermodynamics is not bound to any 
particular mechanism that fulfills its require- 
ments. Diffusion is only a very convenient 
method for attaining equilibrium and does not 
affect the fundamental thermodynamical prop- 
erties such as internal energy, entropy, etc. 
Therefore it is clear that the practical absence of 
diffusion in crystals at ordinary temperatures is 
only a matter of experimental method of attain- 
ing equilibrium and not of the inapplicability 
in principle of thermodynamics to the crystalline 
mixtures. 

In fact in recent years Professor V. Chlopin 
and his collaborators* have worked out several 
methods of attaining equilibrium between mixed 
crystals and solutions, based on the recrystalli- 
zation phenomenon. During the process of re- 
crystallizaton each molecule of salt passes many 
times through the solution and, if this process is 
repeated often enough, an equilibrium between 
crystals and solution is established, and all the 
resulting crystals grow up from a solution of 
constant composition. Such crystals are quite 
homogeneous, and the Berthelot-Nernst distri- 


3 See Tammann, reference 2. 

‘For a summary see: V. Chlopin, A. Polessitsky, A. 
Ratner and P. Tolmatscheff, Ber. 64, 2653 (1931); see also 
R. Mumbrauer, Zeits. f. physik. Chemie A156, 113 (1931). 
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bution law is strictly valid for them. The true 
character of equilibrium attained is proved by 
the fact that the value of the distribution 
constant is the same independently of the side 
from which the equilibrium is attained, i.e., 
whether the distributing salt (usually a radio- 
element in concentrations of the order of 10-!° 
mol./cc) was added to the solution, or to the 
crystals at the beginning. The same result may 
be attained also more conveniently, especially 
at ordinary temperatures, by means of a rapid 
crystallization from a supersaturated solution 
with powerful stirring. Under such conditions, 
the recrystallization of minute crystals goes on 
very fast, and the equilibrium is reached in the 
‘process of growth of the crystals from the 
submicroscopical dimensions. In order to attain 
the equilibrium by crystallization from a super- 
saturated solution without such stirring, the 
crystallization must take place very slowly for 
otherwise instead of the simple form of the 
distribution law, the differential form or the 
so-called Doerner-Hoskins formula is followed. 
This formula is Inx= In y, where x and y are 
the percents of the micro and macrocomponents, 
respectively, remaining in the solution. In that 
case, the concentration ratio remains constant 
for each surface layer of the growing crystals, 
but not for the crystal, as a whole. The same 
formula holds also for crystallization by slow 
evaporation, but the value of \ varies depending 
on the rate of evaporation and stirring, reaching, 
when the latter is good enough, the value of the 
partition factor D= KyC,z'/d, where Ky is defined 
as above, and C,’ is the concentration of the 
macrocomponent in its saturated liquid solution. 
_ These experiments have shown the possibility 
of attainment of an equilibrium state between 
a solid and a liquid solution, and have furnished 
a basis for further theoretical treatment of this 
subject. Such treatment is required in the first 
place in order to find the influence of the compo- 
sition of phases on the value of the distribution 
constant. Such an influence was found experi- 
mentally by Professor V. Chlopin and his co- 
’ workers,’ but the relations are too complicated to 
be made clear by purely experimental methods. 

5 V. Chlopin and A. Polessitsky, Zeits. f. physik. Chemie 


A145, 57 (1929); V. Chlopin and A. Ratner, Comptes 
Rendus Acad. Sciences U. S. S. R. 723 (1930). 
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Since the work of van’t Hoff and Roozeboom,*® 
the question has not been subjected to theoretical 
treatment, except in two works, one by K. 
Herzfeld on thermodynamics of the formation of 
mixed crystals’ and the other by L. Imre on the 
absorption view of the distribution.*® 

In the first work for which the distribution 
was only of secondary importance, the influence 
of the activity coefficients, so important in 
concentrated salt solutions, were lost sight of 
completely. Also the equation obtained is very 
complicated, because of the correction terms 
which have no importance for the distribution 
phenomena. 

The second work deals only with the cases 
where the Doerner-Hoskins equation is valid, 
and L. Imre thinks, that in this case we are 
dealing with quite another phenomenon than 
the distribution attained by recrystallization, 
and that, due to the greater speed of the process, 
the adsorption equilibrium ‘‘freezes in’’ and that 
the work of ‘‘building on’’ (Anlagerung) and not 
that of “building into’ (Einlagerung) is the 
more important. But these assumptions seem 
not to be true, for the range of rates of crystalli- 
zation is the same for both cases. Which of the 
two forms of the distribution law holds depends 
not on the rate of crystallization, but on other 
conditions? If the number of crystals finally 
present is small compared with the number of 
crystallization centers originally formed, the 
crystals are homogeneous and the simple form 
of the distribution law is valid. But if the number 
of crystals is in the main the same during the 
entire process of crystallization, as in the evapo- 
ration experiments, or if their number diminishes 
by binding together of subcrystallites, and not 
by repeated recrystallization, the differential 
form of the distribution law holds, and if diffusion 
is eliminated, the value of \ is the same as that 
of the partition factor D. When the crystalli- 
zation takes place without stirring, the equi- 
librium is not attained due to the slowness of 
diffusion and the observed values of \ lie nearer 
to 1 than corresponding values of D. Thus we 
have no reason to use different laws for different 


cases of crystallization and to take into account 


6 B. Roozeboom, Zeits. f. physik. Chemie 8, 504 (1891). 
7K. Herzfeld, Zeits. f. Physik 16, 84 (1923). 
8 L. Imre, Zeits. f. physik. Chemie A164, 364 (1933). 
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the work of “building on’’ instead of “building 
into”’ even when the Doerner-Hoskins formula 
is followed. L. Imre himself takes as the base of 
his calculation the values of heat of solution of 
the components when calculating the value of \ 
from his equation. But the heat of solution 
corresponds to the work of bringing the salt into 
the solution from the volume of the crystal and 
not from its surface, in other words, to the work 
of ‘‘building into’? and not of “building on.” 
Obviously the right order of the value found 
for \ based on those heats of solution speaks 
against L. Imre’s ideas of “freezing in’ of 
absorption equilibrium rather than for them. 
Because of these considerations I took as a 
starting point for my work the thermodynamical 
treatment of the equilibrium state based on the 
general theory of chemical potentials and on 
modern activity theory of strong electrolytes 
and not the kinetics of crystallization. This 
allowed me to find definite connections between 
the value of D, the composition of the solution 
and the properties of the pure components. 


THEORETICAL PART 
Notation 


In this work the following notations are used: 
x=mol fraction; c=concentration in grams per 
liter; m= molality; Z=solubility product; a,, a_ 
=activities of the ions; a=activity of the 
electrolyte as a whole=a,’+a_’-; y4, y-=ac- 
tivity coefficients of the ions; y=activity co- 
efficient of the electrolyte; =Gibbs’s chemical 
potential; y=chemical potential in the standard 
state chosen in the usual way for a solute;° 
d=specific weight of the crystals; D=partition 
factor = KyC’/d=C,C'/Czd. The subscripts s and 
L refer to solid and liquid phases, respectively. 
The superscript ’ refers to the macrocomponents. 


DEDUCTION OF THE FORMULAE 


According to the main statements of the 
theory of the chemical potential its value for 
the distributing substance at equilibrium must 
be the same in both phases: 


(1) 


*Lewis and Randall, Thermodynamics, McGraw-Hill 
Book Co., New York, 1923, p. 259. 


As in the solid phase we have a very dilute 
solution we can write for ys 


In x,"*. (2) 


The x,’— is omitted for the anion is common for 
both components and the only one in crystals. 
From (1) and (2) we obtain: 


In =(ur—ws)/RT. (3). 


For the saturated solution of the microcompo- 
nent in pure water, we can write 


Mos =ULt+RT In do. (4) 
Substituting wz in (3) from (4) we obtain finally: 
(5) 


The right-hand side of this equation does not 
depend on the composition of the liquid phase; 
therefore we have in it the main expression for 
the distribution constant of an electrolyte be- 
tween a solid and a liquid phase. 

For the most important case practically, i.e., 
when the solid phase is a pure salt isomorphous 
with the microcomponent, the Eq. (5) can be 
transformed into another form more convenient 
for comparison with experiment. This may be 
done by multiplication of both sides of the 
Eq. (5) by the activity of the microcomponent, 
the same in all the solutions and equal to ap, 
the activity in pure saturated solution since they 
are all in equilibrium with the pure solid macro- 
component. Taking into account that a_’ and a_ 
is the same quantity, for it refers to a common 
ion in a common solution, and expressing the 
activities through molalities, we obtain: 

Ay’ 
= — e(Hos—ps)/RT — Do. (6) 
ao 


The quantity Dy) may be called the ‘“‘true 
partition constant,’’ for on one hand, it does not 
depend on the composition of the liquid phase 
and on the other hand it is connected with the 
partition factor D, which is usually employed 
to express the experimental results, and as a 
characteristic of the course of the fractional 
crystallization. We can divide the Eq. (6) into 
two: 


(6a) 


| 
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and 


the first of which expresses the dependence of 
the concentrations of the components in both 
phases on one another, and the second the 
dependence of the value of Do upon the properties 
of the pure components. 


DISCUSSION OF THE RESULTS 


The left-hand side of Eq. (6a) consists of two 
factors: x,’+m,/’+/m,’+ represents the law of 
distribution of an electrolyte as it would be 
if the solutions were ideal, and the second 
y+/"+/y+’+ containing the activity coefficients, 
gives the correction due to deviations of the 
solutions from ideal. We can see from it, that in 
ideal or very dilute solutions the first term must 
’ be constant regardless of the composition changes 
of the solution. For the case v,=1, as in all 
hitherto studied systems, the first factor becomes 
identical with the partition factor D. 

Up till now in all the works on the distribution 
of radioelements, between solid and _ liquid 
phases, it was admitted, that because of the very 
minute concentrations used (of the order of 107° 
mol!/cc or even less) we can consider their 
solutions as ideal and neglect the activity co- 
efficients. But that supposition is only partly 
true; namely, in that the variations in the content 
of the radioelements do not influence the activity 
coefficients of the components, and this is the 
reason for good constancy of the partition factor 
at constant composition of the solution obtained. 
But as the activity coefficient of an electrolyte 
depends on the total ionic strength of the 
solution and not upon the concentration of the 
given electrolyte, and we usually deal with more 
or less concentrated salt solutions, the y is very 
far from unity and depends very much on the 
composition of the liquid phase. 

But, as it is seen from the Eq. (6a) the second 
factor consists not of ys only, but of the ratio 
7+'/y+, and that ratio can be considered as equal 
to unity, or at least as independent of the compo- 
sition in a limited range, even in fairly concen- 
trated solutions, for the ions of the components 
are of the same valency and do not differ very 
much in dimensions for otherwise they could 
not form isomorphous crystals. 
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The range of the probable constancy of D 
becomes even larger if we take into consideration 
that y+ of the Eq. (6a) is the activity coefficient 
of the microcomponent when added to the 
solution of the macrocomponent and not in a 
solution of the same concentration of itself only. 
As some of the factors that determine the 
deviations of the behavior of electrolytes from 
the ideal depend on all the ions in solution 
(electric field in the solution, dielectric constant) 
and some upon the particular ion under con- 
sideration (polarizability, degree of ionization, 
side chemical reactions) we may expect that 
generally the value of the activity coefficient of 
an ion, when added to a much larger quantity 
of another ion would lie nearer to the value of + 
for the latter at the same ionic strength, than 
the value of pure microcomponent.!® Also, the 
changes of the activity coefficient of micro- 
component in such solutions lie nearer to the 
corresponding changes of the activity coefficient 
of the macrocomponent than they would to such 
changes in a pure solution of the microcompo- 
nent. Obviously the difference of the behavior 
of the components must be the smaller, and the 
constancy of the ratio of the y’s the better, the 
more dilute the solution is and the more alike 
the ions of the components are. 

Turning now to the comparison of the Eq. (6a) 
with the experimental data, we meet with the 
lack of knowledge of the activity coefficients in 
such mixed solutions. But the above considera- 
tions show a way for approximate proof of this 
equation, namely, by taking D instead of Dy and 
comparing then the deviations from constancy 
caused by this substitution with what may be 
expected from the electrochemical properties of 
the components. 

In the Tables I, II and III are collected all 
the available data on the variation of the 
distribution constant with the composition of 
the liquid phase. 


10 So it is generally, but not always for the above factors 
change the activity coefficients in different directions, and 
we may find sometimes that as a result of their interaction 
the activity coefficient of the microcomponent would get 
farther from that of macrocomponent instead of ap- 
proaching it. Such is the case for HCl in concentrated 
solutions of LiCl (see J. Hawkins, J. Am. Chem. Soc. 54, 


4480 (1932). 
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TABLE I. Pb(NOs)2—Ra(NOs)2 at 25°.* 


Composition of the solution Kn D 
1.82 m Pb(NOs)> 18.2 2.19 
4.55 (329 26.34.55 


* V. Chlopin and A. Ratner, reference 5. 


TABLE II. BaBr.—RaBr, at 0°.* 


Composition of the solution Ky D 
3.11 m BaBr: 52.8 12.4 
0.27 m HBr 57.7 13.2 
“ 59.1 12.8 
0.64 58.3 11.9 
™ 116.0 15.4 
iar 146.0 14.2 
— * 4.00“ 173.4 9.8 
0.42“ “ a * 243.0 7.8 


*V. Chlopin and B. Nikitin, Zeits. f. anorg. Chemie 166, 
311 (1927). 


TABLE III. Ba(NO;)2—Ra(NOs)2.* 


Composition of the solution a © D 
0.411 m Ba(NOs)2 25° 46.0 1.49 
0.367 0.530 m BaCl._- 21.7 1.56 
0.195 “ 85.2 1.31 
0.674 } 25.5 1.34 
0.400 ‘ 0.209 “ HgCl. 44.7 1.39 
0.416 “ 1.18 ‘ glucose 56.1 1.60 
0.189 “ 5 0° 154.0 2.31 
0.080 “ 0.423 HNO; “334.0 2.12 
0.185 “ 0.202 BaCl, 64.0 2.24 


*V. Chlopin, A. Ratner and P. Tolmatscheff, unpub- 
lished. 


As it may be seen from the tables, the constancy 
of D is in full accord with the above considera- 
tions, the worst in the case of lead nitrate as 
macrocomponent, because of the large solubility 
of lead nitrate and the great difference in 
properties between lead and radium ions and 
the best for the sparingly soluble barium nitrate 
whose cation is very much more similar to the 
radium ion. In the last’case the constancy of D 
is fairly good (within 20 percent) though the Ky 
values vary many times. 

But not only the degree of constancy of D, 
but also the direction of deviations from it is 
in full accord with what could be expected from 
the properties of the components. Thus in the 


case of lead nitrate which is much more associ- 
ated electrolytically than barium and radium 
salts" and gives nondissociated basic acetates 
on addition of HNO; or NaC:H;0O: the activity 
coefficient of lead vy,’ falls a great deal faster 
than that of radium, the ratio y,’/7 falls too, 
and as Dy» must remain constant, D must rise, 
which is the case. For barium salts the conditions 
are reversed, for radium salts must be more 
associated than those of barium because of the 
larger polarizability of radium, and we see in 
Tables II and III that D in these systems 
decreases slightly but systematically in all cases 
on addition of corresponding acids. So it seems 
that the Eq. (6a) is wholly confirmed by the 
experiment. 

Expressing the activities through molalities, 
we can write Eq. (6b) 


Do = (Lo! /Lo) (Yo! 


The first factor on the right-hand side of this 
equation, which would have been the only one 
if both phases were ideal mixtures, gives the 
connection between the partition factor and the 
solubilities of the components. The existence of 
such a connection was admitted by all chemists 
who had worked with the fractional crystalliza- 
tion and it was thought to be so direct that 
conclusions in regard to the identity of the 
solubility of the components” were made from 
an unsuccessful attempt of fractional crystalli- 
zation. But an exact investigation, experimental 
or theoretical, on the existence and form of such 
a connection was still lacking, and in recent 
years there has been accumulated material 
contradicting it, at least in such exact form as it 
was assumed before.'* 

It follows from Eq. (6b) that this connection 
consists of a proportionality to the ratio of the 
solubility products but is masked by the influence 
of the non-ideal behavior of the solutions ex- 
pressed in the liquid by the activity coefficients, 
and for the solid in the term e“*-**)/®7, These 
influences, especially that of the liquid phase, 


1 Plake, Zeits. f. physik. Chemie A162, 24 (1932). 

12M. Curie, Traite de Radoactivite, Gauther, Villard, 
Paris, 1910. V.I. p. 174. 

13 Chlopin, Zeits. f. anorg. Chemie 143, 97 (1923); 
O. E. Erbacher, Ber. 63, 141 (1930); A. Polessitsky, Zeits. 
f. physik. Chemie A161, 325 (1932). 
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may change the value of Dy many times and thus 
the above-mentioned experimental facts are 
explained. 

The comparison of the Eq. (6b) with experi- 
ment is even more difficult than that of Eq. (6a) 
for we meet again with the absence of data for 
the ratio y’/y+, but here we need the constancy 
of that ratio as well as its absolute value. 
Further the term e“*-#*)/27 js also undetermin- 
able as yet, for the theoretical calculation of 
lattice energies is not exact enough and its 
computation from the heats of solution of the 
mixed crystals also fails because of the absence 
of such determinations through a sufficiently 
large range of concentrations. We could pass 
over all these difficulties by means of a determi- 
nation of the partition constant for a weakly 
soluble macrocomponent at different tempera- 
tures, but even with data now available we can 
judge of the validity of the Eq. (6b) by com- 
paring the values of D with the ratios ao’ /ao for 
all the systems, for which the solubilities of the 
components are known. 

These data for t= 20°C are given in Table IV 
in the first column of which the systems con- 
sidered are given, the microcomponent being 
placed in brackets. The solubilities of radium 
salts are taken from the work of O. Erbacher," 
and the values of activity coefficients necessary 


TABLE IV. 

System m'o/mo a'o/ao D 
Ba(Ra)Bre 3.50 1.91 60 11 
BRN 1.71 2.07 16.5 4.5 
Ba(Ra)(NOs3)2 0.33 0.83 1.14 1.6 
Ba(Pb)(NOs)2 0.33 0.20 0.12 0.16 
(NOs)2 1.66 4.15 11.9 
Pb(Ba)(NOs)2 1.66 5.03 8.32 2.4 


“4 Q. Erbacher, Ber. 63, 141 (1930). 


for the figures of the fourth column were evalu- 
ated approximately by a graphical extrapolation, 
using Lucass’!® data on alkaline earths chlorides 
for the halide systems, and Gellbach and 
Heppke’s'* on barium bromide, and the author’s 
determinations at 70° on the nitrates (by vapor 
pressure methods, unpublished). 

From Table IV we see that the partition 
factors show an unmistakable parallelism to the 
activity ratios, though they are by no means 
equal to the theoretical values, the discrepancy 
increasing with increased solubility of the macro- 
component. At the same time, it is difficult to 
speak even of a parallelism between the solubility 
ratios and the partition factors. 

We may go slightly further for the system 
Ba(Pb)(NO3)2 due to the determination of the 
heats of solution of mixed crystals of barium 
and lead nitrates made by G. Tamman,!’ which 
enables us to compute at least the sign and order 
of magnitude of uo,—us. This computation gives 
the value of 400-600 cal., from which 


Do do’ BT = 0.30. 


The y4’'/y+ for this system is equal approxi- 
mately to 1.30 and the experimental value of Do 
is approximately 0.20. The agreement between 
these values may be considered as quite satis- 
factory, taking into account the very approxi- 
mate character of the calculations. 

At this point I wish to take the opportunity 
to acknowledge my gratitude to Professor V. 
Chlopin for his advice and continued attention 
throughout this work. 


18 Lucass, J. Am. Chem. Soc. 47, 743 (1925). 

16 Gellbach and Heppke, J. Am. Chem. Soc. 48, 1504 
(1926). 

17 G, Tammann and Krings, Zeits. f. anorg. Chemie 130, 
229 (1923). 
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Correlation of Heat Capacity, Absolute Temperature and Entropy 


GeorcGE A. LinHART, Riverside Junior College, Riverside, California 
(Received August 7, 1933) 


Correlation of heat capacity, absolute temperature and 
entropy: Starting with the two assumptions that at very 
low temperature the heat capacity is proportional to the 
entropy, and that at relatively high temperature it 
approaches a limiting value (C,,), very simple relations 
are obtained between the heat capacity and the absolute 
temperature, and between the entropy and the absolute 


temperature of a given substance, by the introduction of 
a simple algebraic function which satisfies the two as- 
sumptions. The validity of these relations is substantiated 
by the excellent agreement of the observed and the 
calculated values for the heat capacities and for the 
entropies of a considerable number of substances, ten of 
which are given in this paper. 


N establishing the relation of heat capacity, 
absolute temperature and entropy of a given 
substance at constant pressure, two assumptions 
are made: first, the heat capacity C is propor- 
tional to the entropy S as the absolute tempera- 
ture approaches zero; second, the heat capacity 
C approaches a limiting value C,, at relatively 
high temperature. These assumptions may be 
expressed analytically as follows: 


(a) (8C/dS),=K(T-0), 
d 
(b) 


Making use of the thermodynamic equation, 


(0S/8T) >= C,/T, 


gives 


Arithmetical experience with a large number 
of heat capacity data shows that the left-hand 
member of Eq. (1) may be replaced by a very 
simple function, namely, 


(2) 
so that Eq. (1) becomes, 
(3) 


The right-hand member of Eq. (2) obviously 
satisfies the two assumptions (a) and (b). 
The integrated form of Eq. (3) is, 


log C/(C,,— C) =K log T+log k, (4) 


which may also be written in the forms, 
(5) 
C=C.,ke™ 7/(1+keX 7), (6) 


In testing the validity of Eq. (4), the following 
method was used. 

First, the experimental values for C were 
plotted against the corresponding logarithms of 
the absolute temperatures, and the curves 
representing Eq. (6) sketched in. 

Second, the best value for C,, was in each 
case located on the curve. This was easily 
attained, since the curve is symmetrical with 
respect to the inflection point, as may be readily 
shown by placing the second derivative of C with 
respect to log T, of Eq. (6) or (7), equal to zero; 
at this point, C=C,,/2. 

Third, the best values for K and for logk 
were obtained from the straight line plot of Eq. 
(4), where K is the slope and log k is the intercept 
on the log C/(C,,—C) axis. 

Fourth, with the aid of these quantities, the 
values for C “calculated” may be readily ob- 
tained by means of Eq. (4). 


or, 


DISCUSSION OF C,, 


From the analogy of the specific heat of a 
“perfect monatomic gas,” the older school of 
physical chemists reached the conclusion that 
C., for all solids at constant volume should be 
3R (5.966), calories per gram atom of substance. 
Any excess of this value was attributed to the 
work done by virtue of the expansion of the 
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substance. Modern investigators have found, 
however, that the rather large excess of 3R 
calories near the melting point cannot be ac- 
counted for in this manner. Some suggest that 
at the upper temperature, or near the melting 
point of the substance, the electrons themselves 
may acquire an appreciable heat capacity; but 
the specific heat data in the region of the melting 
point are still meager and of insufficient accuracy 
to warrant calculation of the heat capacity of 
the electrons. Investigations in this field are 
now being conducted by Jaeger and Bottema,' 


1 Jaeger and Bottema, Proc. Acad. Sci. (Amsterdam), 
1932. 


and their results will very likely lead to definite 
conclusions. 

Similar investigations in the neighborhood of 
the absolute zero are now being carried on by 
Keesom and his associates.” When sufficient data 
have been accumulated in these two fields they 
will be analyzed in a similar manner. 


THE GRAPHS 


The graphs are self-explanatory. The letters 
M.P. denote the melting point; the letters S.P. 
denote the sublimation point. It is of interest 


2 Keesom, Proc. Acad. Sci. (Amsterdam), 1932. 
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TABLE I. Summary of entropies from 0 to 298° absolute, and constants. 


Copper 8.10 8.18 5.97 2.900 —5.49550 
Tin (white) 12.22 11.17 6.40 2.074 —3.41445 
Tin (gray) 10.82 9.23 6.70 1.700 — 3.03844 
Silicon 4.25 4.70 6.30 2.316 — 5.15068 
Graphite 1.45 1.30 6.40 1.820 —4.79222 
Magnesium 8.55 8.30 6.60 2.050 — 3.94848 
Mercurous chloride 23.88 23.20 6.50 1.480 —2.52140 
Potassium chloride 20.50 19.70 6.40 2.380 — 4.23445 
Lead chloride 34.86 33.20 6.40 1.769 — 2.99838 
Lead bromide 39.02 39.75 5.50 1.840 — 2.96575 


to note that the experimental points follow the 
calculated curves very closely, and that the 
average deviations do not exceed the average 
experimental variations. 


CONCLUSION 


In conclusion it may be mentioned that the 
heat capacities of all substances, measured by 
the thermoelectric method, conform statistically 
to Eq. (6). The ten substances given in this 
paper were taken at random from a group of 
about three times that number. It is also 


interesting to find that the entropies calculated 
from the integrated form of Eq. (2) between 0 
and 298° absolute are practically identical with 
those given in the thermodynamic literature, 
which have been obtained by the mechanical 
integration of the areas under the C—log T 


curves. A comparison of these values is given in 
Table I. In the second column are given the 
values for the entropies from 0 to 298° absolute 
calculated from the integrated form of 2, 


S=2.3(C,,/K) log C,/(C,-C), (7) 


and in the third column are given the values as 
recorded in the thermodynamic literature; these 
have been obtained by mechanical integration. 
It is seen that in every one of the ten sub- 
stances the observed values agree with the 
calculated values within experimental error. 


Note: References to the data used in this paper may be 
found in the International Critical Tables, Vol. V. The 
data recorded in that volume are not, however, the original 
values and are therefore not suitable for investigational 
purposes. For an extensive and interesting discussion of 
specific heat data at constant volume, the reader is referred 
to an article by F. Schwers, Phys. Rev. 8, 117 (1916). 
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The Thermal Decomposition of Propyl Mercaptan 
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The thermal decomposition of normal propyl mercaptan 
over a temperature range from 405 to 435°C and pressures 
from 100 to 300 mm is homogeneous on a poisoned glass 
surface. The reaction is complex showing an induction 
period with auto-accelerating rate, ending finally in a 
pseudo-equilibrium involving propylene, hydrogen sulphide 


and isopropyl mercaptan. The energy of activation of the 
induction period is about 40,000 calories and that of the 
later reaction 39,000 calories. A mechanism, has been 
proposed which involves primarily the formation of propyl 
sulphide reacting further to form a sulphonium hydro- 
sulphide as an unstable intermediate. 


ROM a recent study of the decomposition of 

ethyl mercaptan Trenner and Taylor? con- 
cluded that the mechanism involved a series of 
consecutive reactions leading to an equilibrium in 
which ethylene and hydrogen sulphide pre- 
ponderate. The marked similarity in the de- 
composition of ethyl sulphide suggested a parallel 
mechanism and the same intermediate was 
tentatively assigned to both reaction schemes. 
The isolation of this intermediate could not be 
accomplished as definite proof of the conclusions 
drawn. Propyl mercaptan being homologous with 
the ethyl derivative a similarity in their de- 
compositions might logically be expected whilst 
the individual rates of the consecutive reactions 
might be such as to allow the isolation of the 
intermediate compound for further study. Al- 
though this has not been achieved the work has 
served to verify, amplify or in some ways modify 
the conclusions previously drawn. 


EXPERIMENTAL 


The sample of normal propyl mercaptan used 
was from the Eastman Kodak Company and was 
carefully dried and fractionated to give a product 
boiling in agreement with the findings of Ellis and 
Reid’ from 67.4 to 67.6°C at 763 mm. 


! Abstract from a thesis presented in partial fulfillment 
of the requirements for the degree of Doctor of Philosophy 
at New York University, June, 1933. 

? Trenner and Taylor, J. Chem. Phys. 1, 77 (1933). 

3 Ellis and Reid, J. Am. Chem. Soc. 54, 1674 (1932). 


A. Static method 


The apparatus used in the work was practically 
identical with that used by Trenner and Taylor. 
The reaction vessel of Pyrex glass had a capacity 
of 279 cc. The total free volume outside the 
furnace was less than two percent of this value 
and was maintained by auxiliary heaters at 65°C 
to prevent condensation of mercaptan. As was 
found with ethyl mercaptan the induction 
periods in the decomposition of the propy| 
mercaptan were somewhat erratic in duration on 
a clean surface but after three or four samples 
were successively decomposed, the products of 
each being removed by five minutes pumping 
with a “hyvac”’ oil pump, the induction periods 
as well as the entire decomposition rates could be 
reproduced at will. The results for various 
pressures and temperatures will be presented 
later. 


B. Dynamic method 


To obtain a qualitative knowledge of the 
products of reaction particularly during the early 
stages, as substantiation of the proposed mechan- 
ism, several dynamic experiments were made. 
The propyl mercaptan was maintained at a 
temperature of 78°C with ethyl alcohol constantly 
refluxing, to exert a constant vapor pressure on a 
calibrated capillary connected with the reaction 
tube in the furnace. The temperature was 
measured by a platinum resistance thermometer 
located in the center of the reaction tube. Passing 
out of the reaction tube the gases were carried 
through three traps. The first trap was immersed 


798 


DECOMPOSITION OF 


in a Dewar flask containing ice and salt at 
— 19°C; the second was immersed in solid carbon 
dioxide and toluene and the third in liquid air. 
This latter trap was connected to a mercury 
manometer and gas burette. The pressure on the 
reservoir side of the calibrated capillary was thus 
constant at 1066 mm, the rest of the system being 
maintained at atmospheric pressure. That the 
system should be entirely free from oxygen, it 
was first swept out with helium. During a run as 
the helium was displaced it was allowed to blow 
off at atmospheric pressure, whilst any small 
quantities remaining would not interfere with a 
qualitative analysis of other gases produced by 
reaction. 
RESULTS 


Since there are no data available on the vapor 
pressure of propyl mercaptan and since the 
apparatus for the static experiments lent itself 
admirably to this determination the values given 
in Table I were obtained. 


TABLE I. 

Temperature Vapor pressure 
284.6°K 318.1°K 84.0 mm 351.0 mm 
289.9 323.3 97.5 425.0 
295.8 328.2 141.0 509.5 
301.7 333.3 182.4 600.0 
307.5 340.57 (calc’d) 233.0 760 
313.0 284.2 


Fig. 1 shows these results graphically where 
the logarithm of the pressure is plotted against 
the reciprocal of the absolute temperature. The 
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Fic. 1. Vapor pressure curve of propyl mercaptan. 
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straight line may be represented by the following 
equation: 


log p= —1647/T+7.7190. 


A. Static method 


By this method data on the rate of decompo- 
sition were obtained at varying pressures from 
100 to 350 mm and at temperatures from 405 to 
435°C. Fig. 2 illustrates graphically the rates 
obtained at initial pressures of 101, 152, 199, 300 
and 351 mm at 405°C. It will be observed that 
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Fic. 2. Decomposition rate curves of propyl mercaptan at 
pressures of 101, 152, 199, 300 and 351 mm at 405°C. 


the curves are identical in form to those obtained 
by Trenner and Taylor for ethyl mercaptan. 
Since it was there shown that the point of 
inflection was of considerable importance, its 
exact location was made experimentally by 
taking pressure readings every fifteen seconds in 
the neighborhood of the approximate position of 
maximum rate previously determined. Table II 
gives the location of the points of inflection at the 
various pressures and temperatures studied. 
The last two columns will be referred to later in 
the discussion. 

The peculiar behavior of the induction period 
previously mentioned, suggesting heterogeneity, 
the effect of an increased surface was tested. 
Results again identical with those for ethyl 
mercaptan were found. The first reaction after 
the chamber had been filled with short lengths of 
Pyrex tubing and very thoroughly evacuated 
showed practically no induction period. In 
subsequent reactions this became progressively 
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TABLE II. TABLE III. 
Init. Time to Press. Init. 
press. inflection Press. unreacted press. Press. Percent tos tso 
(mm) (min.) increase Py P;/P | (mm) increase decomp. (min.) 
Temperature 405°C Temperature 405°C 
351 8.75 132 219 0.624 | 351 276 78.6 1.13 2.55 
300 10.00" 107.6 192.4 0.641 300 243 81.0 1.26 2.91 
199 12.75 71.8 127.2 0.638 199 164 82.0 1.42 3.30 
152 17.50 60.7 91.3 0.601 152 126 82.8 1.72 4.27 
101 25.00 41.5 59.5 0.589 101 86 85.3 1.85 4.35 
Temperature 415°C Temperature 415°C 
351 5.50 134 217 0.618 | 351 279 79.5 0.76 1.75 
304 6.12 107.5 196.5 0.646 | 304 247 81.3 0.80 1.84 
202 9.00 85 117 0.585 202 166 82.3 0.97 2.32 
149 13.00 63.5 85.5 0.574 149 125 83.9 1.30 2.75 
99 19.00 44.7 54.3 0.554 99 86 86.0 1.56 3.42 
Temperature 425°C Temperature 425°C 
351 3.75 132 219 0.624 351 283 80.7 0.50 1.15 
297 4.00 111.9 185.1 0.624 | 297 242 81.5 0.55 1.29 
200 5.75 88 112 0.560 | 200 165 82.5 0.62 1.62 
147 9.00 63.6 83.4 0.569 147 125 85.0 0.92 2.47 
101 12.00 45.5 55.5 0.550 101 87 87.0 1.12 2.78 
Temperature 435°C Temperature 435°C 
347 2.50 133 214 0.617 347 281 81.2 0.36 0.85 
298 2.75 109 189 0.634 298 244 81.7 0.37 0.85 
202 3.75 85.9 116.1 0.581 202 168 83.2 0.42 1.02 
150 6.00 67.2 82.8 0.552 150 130 87.0 0.60 1.77 
100 8.00 43.8 56.2 0.562 100 89 89.0 0.91 2.1 


longer eventually reaching a constancy to dupli- 
cate that found for the empty vessel. This same 
“aging” effect was shown also when finely 
powdered glass was used to increase the surface. 
It is apparent therefore that the reaction is one 
which may be catalyzed under suitable conditions 
but that under the conditions prevalent in this 
work is strictly homogeneous. The slope of the 
rate curves subsequent to the inflection was 
always constant even for reactions in a new vessel 
before aging. This fact would suggest that the 
secondary reactions are always homogeneous. 
That the primary reaction becomes homogeneous 
may be connected with the observed darkening of 
the reaction bulb by a dark brown deposit which 
could be distilled slowly under vacuum and which 
therefore could not be carbon but in all proba- 
bility was a complex polysulphide. 

An examination of the end-points of the 
reaction shows a variation with temperature and 
pressure. The greatest percentage decomposition 
occurred at the highest temperature and lowest 
initial pressure; the least percentage decompo- 
sition at the lowest temperature and highest 


pressure. This variation is apparent from the 
data in Table III and would suggest the existence 
of an equilibrium similar to that established by 
Trenner and Taylor for ethyl mercaptan, ethyl- 
ene and hydrogen sulphide. In an attempt to 
verify the existence of this equilibrium for 
propylene and hydrogen sulphide a mixture of 
350 mm of the latter and 355 mm of the former 
were allowed to react at 435°C. Such a mixture 
should be equivalent to 350 mm propyl mer- 
captan which decomposed under these conditions 
with a pressure increase of 81.2 percent. If the 
suggested equilibrium existed, the mixture should 
react with a pressure decrease of 18.2 percent. 
Actually a maximum decrease of only 10 percent 
was obtained after twelve hours, after which 
time the pressure began to increase again very 
slowly. Duplicate runs always showed the same 
tendency though quantitative reproducibility 
was small suggestive of some heterogeneity in 
the complicating reactions. 

The effects of foreign gases were determined in 
eleven cases with oxygen, ethylene, propylene, 
sulphur dioxide, carbon dioxide, nitrogen, helium. 
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hydrogen sulphide, ethane, hydrogen and am- 
monia. In each instance, oxygen excepted, 100 
mm of gas were added to 200 mm of mercaptan 
vapor. With oxygen only 5 mm were used 
because of the pronounced effect it had on the 
length of the induction period. Fig. 3 shows eight 
of the eleven cases studied, the normal rate curve 
being included for comparison. The curve for 
sulphur dioxide is practically coincident with 
that for carbon dioxide, whilst helium was 
apparently without effect. The striking feature of 
Fig. 3 is the parallelism of the curves beyond the 
inflection in each case, ethylene and propylene 
excepted, the sole differences lying in the lengths 
of the induction periods. The more specific 
effects of the products of decomposition were 
determined by allowing a sample of mercaptan 
partially to decompose and then adding a fresh 
sample and following the reaction to the end. By 
knowing from previous experiment the shape of 


the reaction rate curve for the initial sample, the 
rate curve for the injected sample may be 
synthesized by subtracting the ordinates of the 
curve for the initial sample alone from those of 
the actual observed curve. Three such tests were 
made and checked: first, the fresh reactant was 
added before the original sample had reached the 
maximum reaction rate, that is, early in the 
induction period; second, the addition made at 
the point of inflection; and third, injecting the 
second sample when the first had practically 
completely decomposed. The results of these 
three tests are condensed on one diagram in Fig. 4 
for the sake of simplifying comparison. The 
points where the injections were made have been 
taken as origins for axes for the injected sample. 
The continuous curve marked “normal” repre- 
sents the course of decomposition of 181 mm of 
mercaptan at 405°C. The first injection of 172 
mm was made when the original 181 mm had 
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Fic. 4. Composite diagram of the effects of three separate additions of propyl mercaptan to 181 mm 
already partially decomposed at 405°C. Sample 1 of 172 mm injected after 7 min. Sample 2 of 177 mm 
injected after 14 min. and sample 3 of 181 mm, after 30 min. 


reacted for seven minutes. The curve labelled 
“actual’”’ represents the course of decomposition 
of the combined samples. By subtracting the 
ordinates of the normal curve from those of the 
actual curve the synthesized curve for the 
injected sample is found and is shown dotted in 
the diagram. In the second case 177 mm were 
injected after 181 mm had decomposed 14 
minutes and in the third case 181 mm were 
injected at the thirtieth minute of decomposition 
of the original 181 mm. It should be noted that 
the induction period of the injected sample is in 
each instance shorter than that of the original 
sample. The most marked shortening occurs 
however in the sample injected at the point of 
inflection of the original sample. 

It was shown by Trenner and Taylor that the 
decompositions of ethyl mercaptan and ethyl 
sulphide proceeded analogously and were in all 
probability intimately connected. Since propyl 
sulphide boils at 142°C it was impossible with the 
present apparatus to make any extended study 
of its decomposition due to the inability of 


keeping the stopcocks and capillaries at this 
elevated temperature. Nevertheless by main- 
taining the capillaries at 75°C a few observations 
at initial pressures between 30 and 60 mm and at 
temperatures between 435 and 455°C were made. 
Examples of the curves obtained are shown in 
Figs. 5 and 6. The lower curve in Fig. 5 is for the 
decomposition of 31 mm of propyl sulphide at 
435°C. The shape is essentially the same as that 
found for low pressures of mercaptan. Although 
the marked flattening at these pressures appears 
to mask the presence of the induction period the 
fact is that the curve is practically all induction 
period, the rate of reaction never becoming 
sufficiently great to produce a noticeable point of 
inflection. There is no doubt that at higher 
pressures the curves would show the same 
hysteresis shape that mercaptan shows. This can 
be seen from Fig. 6 where curve 1 is for 30 mm of 
mercaptan and curve 2 for 31 mm of sulphide 
both at 435°C where the mercaptan curve is s0 
flattened as to indicate an apparently smooth 
decomposition free from induction. 
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Fic. 5. Decomposition rate curves of propyl sulphide (lower) and propyl sulphide and hydrogen 
sulphide (upper). 
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Fic. 6. Rate curves for the decomposition of (1) propyl mercaptan, (2) propyl sulphide, (4) a mixture 
of both. Curve (3) for comparison is the result of adding the ordinates of (1) and (2). 


Since with ethyl sulphide it was shown that the 
presence of hydrogen sulphide completely re- 
moved the induction period a comparison was 
made in this case with propyl sulphide. Fig. 5 
shows the result at 435°C, the upper curve 
representing the pressure change with time of 30 
mm each of propyl sulphide and hydrogen 
sulphide. The lower curve for comparison is that 
of 31 mm propyl sulphide alone. The shape of the 
former curve is entirely different from that of the 
latter; there is no induction period, the reaction 
beginning at maximum rate. It should be noted in 
passing that hydrogen sulphide alone is stable at 
435°C, 

Since mercaptans are in reality substituted 
hydrogen sulphides it might be expected that 
mixtures of sulphides with mercaptan would 
behave as did mixtures with hydrogen sulphide. 
In Fig. 6 are given the results observed at 435°C 


for a mixture of 32 mm propyl sulphide and 31 
mm propyl mercaptan indicated by curve 4. 
Curves 1 and 2 as already stated are for 30 mm 
mercaptan and 31 mm sulphide alone, respec- 
tively. The dotted line, curve 3, is simply the sum 
of the ordinates of curves 1 and 2 and is thus the 
pressure change which would be expected for 
the mixture if each constituent reacted indi- 
vidually. The marked difference between the 
observed curve 4 and the ‘‘theoretical’”’ curve 3 
must surely indicate an interaction in some way 
between the sulphide and mercaptan as between 
the sulphide and hydrogen sulphide. Similar 
results were found also at 445° and 455°C. 


B. Dynamic method 


Experiments were made at temperatures of 
- 405°, 435° and 440°C. The contact times under 
the conditions previously given were calculated 
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to be 0.97, 0.93 and 0.92 minute, respectively, at 
these temperatures, all estimated to be well 
within the induction period of a reaction with an 
initial pressure of 760 mm. 

In the first run at 405°C all of the reaction 
products condensed in the first trap at — 19°C as 
a colorless liquid which distillation proved to be 
unchanged mercaptan. There had been therefore 
practically no reaction. At 435°C the majority of 
' the products condensed again in the first trap, 
there being, however, a small amount of liquid in 
the carbon dioxide trap. This latter proved to be 
only liquid hydrogen sulphide by the usual tests. 
The liquid in the first trap was washed with 
sodium hydroxide, most of it dissolving. The 
alkali-insoluble, oily residue was separated by 
extraction with ether and the ether removed by 
distillation. An amount of a brownish liquid 
remained, too small to be tested but having the 
disagreeable odor of propyl sulphide. The alkali- 
soluble portion was acidified with dilute hydro- 
chloric acid liberating a lemon colored liquid 
which practically all distilled in the range for 
propyl mercaptan leaving a small amount of a 
tarry residue which showed no boiling point but 
seemed to decompose at about 200°C. 

At 440°C the first trap at —19°C contained a 
large amount of a brown liquid, the second at 
—78°C, a considerable quantity of a colorless 
liquid, the liquid air trap was empty and no 
uncondensed gases were collected in the gas 
burette. The contents of the first trap were 
treated as before. The alkali-insoluble portion, 
again a brownish oily liquid with an odor of 
propyl sulphide, was treated with mercuric 
chloride solution to yield a finely divided 
crystalline precipitate of buff color, which was 
taken to be the white addition compound of 
propyl sulphide and mercuric chloride colored 
with traces of polysulphides, the presence of 
which was indicated by a distillation with 
decomposition above 200°C. The alkali-soluble 
portion proved again to be unchanged mercaptan 
contaminated also with polysulphides which 
decomposed about 200°C. The colorless liquid in 
the carbon dioxide trap proved to be hydrogen 
sulphide and hydrocarbons. The hydrogen sul- 
phide was allowed partially to distill off and the 
trap then connected to an Orsat gas analyzer. A 
large fraction of the gas was absorbed by sodium 
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hydroxide, another large fraction by bromine 
water and the small remainder by absolute 
alcohol. The gas removed by sodium hydroxide 
was undoubtedly hydrogen sulphide. The bro- 
mine water would remove propylene whilst 
absolute alcohol would take up propane. Boiling 
points will substantiate these conclusions since 
that of hydrogen sulphide is —59°C, that of 
propylene —47°C and of propane —44.5°C. 
Lower paraffins and olefines are excluded since 
they would not condense at —78°C. It must be 
admitted that higher hydrocarbons would be 
retained but their probability of formation seems 
remote. 


DIscUSSION OF MECHANISM 


The marked analogy between the foregoing 
data and those found by Trenner and Taylor for 
ethyl mercaptan would suggest a similar mecha- 
nism of decomposition for propyl mercaptan. The 
presence of the induction period is indicative of at 
least two reactions occurring during the de- 
composition and most obviously a distinction can 
be made between the two portions of the reaction 
up to and succeeding the point of inflection. 
That this latter is a comparable point for all the 
cases studied can be seen from a perusal of the 
last column of Table II in which the ratio of the 
pressure of mercaptan undecomposed to the 
initial pressure is constant and independent of 
the temperature. This fact makes it possible to 
study separately the course of the reaction in its 
later stages with the point of inflection as origin 
for this portion. It was on this basis that the 
fractional lives given in Table III were calcu- 
lated. Because of the variable end-point of the 
reaction at different pressures and temperatures, 
an approximate correction has been applied to 
make the times more comparable from case to 
case. To do this the end-point of the reaction at 
the highest pressure for each temperature was 
chosen as standard and all the other end-points 
at this temperature corrected to this value by 
multiplying each by the ratio of the percentage 
increase of this standard reaction to that of the 
particular reaction to be corrected. The quarter 
and half lives so found will be seen to show 
approximate constancy with, however, a defi- 


-nitely rising trend with decrease in initial pres- 


sure at each temperature. This would lead one to 
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suspect that this secondary portion of the re- 
action constituted a pure unimolecular reaction 
which deviated from unimolecularity around 300 
mm. It is probable however as will be shown that 
even this portion of the overall rate curve is not 
typical of a simple reaction and that therefore 
these fractional lives are themselves complex. 

The energy of activation of this secondary part 
of the reaction calculated from the above 
fractional lives by means of the Arrhenius 
equation is found to have an average value of 
39,000 calories. In a similar manner using the 
times taken to reach the points of inflection the 
energy of activation of the induction period may 
be calculated to be on the average 40,000 
calories. The agreement between these two values 
is striking when it is considered that an error of 
one-hundredth of a minute in the fractional life 
will change the value of the energy of activation 
by almost 1000 calories. The agreement also has 
as a consequence the impossibility of separating 
the primary and secondary reactions by a 
temperature variation. The energy of activation 
of the reaction of mixtures of propyl mercaptan 
and propyl sulphide was calculated, using the 
quarter life, to be 39,500 calories, also in striking 
agreement with the previous values. 

On the basis of these facts it is possible to 
postulate a probable mechanism for the reaction 
which will account for all the observations and at 
the same time suggest where complicating side 
reactions might mask a complete understanding 
of the reaction. It will be convenient to consider 
the reactions postulated individually, for since 
the mechanism formally resembles that of 
Trenner and Taylor for ethyl mercaptan some 
amplification is possible. 

In the first place the induction period is 
viewed as being due to the bimolecular reaction 
of two propyl mercaptan molecules to give 
propyl sulphide and hydrogen sulphide without 
volume change: 


C3H; 
2C;H;SH> 
C;H; 


S+H3S. 


Well known as a heterogeneous reaction on 
metallic sulphide surfaces it is probable that a 
clean glass surface would catalyze it and cause 
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the observed erratic values in the length of the 
induction period. That hydrogen sulphide is 
formed early in the reaction has been shown by 
the results from the dynamic method at 435°C 
where it was further shown (the second case) to 
be present without any accompanying propylene 
or other hydrocarbon. Furthermore it may be 
recalled a small amount of an alkali-insoluble 
portion was also found. To split hydrogen 
sulphide from a mercaptan without forming any 
olefine would seem to require the above reaction 
which is in keeping with numerous previous 
findings. 

For the next step, the observation that both 
hydrogen sulphide and propyl mercaptan were 
capable of removing completely the induction 
period of the decomposition of propyl sulphide, 
must mean that a definite chemical reaction 
occurs in which an addition compound is formed. 
Trenner and Taylor suggested such an inter- 
mediate formed from ethyl sulphide and hydro- 
gen sulphide. Since propyl mercaptan has been 
shown to remove the induction period of the 
sulphide and since it must still be present in 
considerable quantities early in the reaction, 
there would appear to be two competitors for the 
propyl sulphide formed in the first step as shown 
by the following equations: 


C;H; H 
S (a) 
C3H; C3H, SH 
C;H; C3H; 
\ 
S+C;H;SH> S (b) 
C3;H; C;H; SH 


The intermediate in Eq. (a), analogous to that 
which Trenner and Taylor postulated and at- 
tempted unsuccessfully to prepare in quantity, 
appears at first sight to be improbable. It is 
however, only by comparison with the compound 
in Eq. (b) that its possibility becomes apparent, 
for this latter compound would simply be 
tripropyl sulphonium mercaptan. Now although 
no such compounds are definitely reported in the 
literature there is every indication that they do 
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exist. There has been reported‘ an analogous 
selenium compound (C2H;)2Se.SeH. The sul- 
phonium halides and hydroxides are well known, 
the latter being strong bases. Bernthsen and 
Sudborough® state specifically ‘‘trimethyl sul- 
phonium hydroxide . . . gives salts with acids 
even with hydrogen sulphide, these latter closely 
resemble the alkali sulphides.’’ Crum-Brown and 
Blaikie® succeeded in preparing trimethyl sul- 
phonium sulphide [(CH3)3S ]eS from equivalent 
amounts of trimethyl sulphonium hydroxide and 
hydrogen sulphide. Attempts to concentrate an 
aqueous solution by drying over phosphorus 
pentoxide were successful only to a certain point 
however, beyond which decomposition to di- 
methyl sulphide occurred. 

In view of these facts a synthesis of triethyl 
sulphonium hydrosulphide was attempted.’ Tri- 
ethy! sulphonium bromide was treated with 
sodium hydrosulphide in absolute alcohol at 
35°C. After some time the mixture assumed a 
yellow color and smelled strongly of mercaptan. 
A fractional distillation yielded ethyl mercaptan 
and diethyl sulphide. The presence of these com- 
pounds is in reality positive evidence that the 
desired sulphonium hydrosulphide had first been 
formed but decomposed rapidly. A second 
attempt using liquid hydrogen sulphide as sol- 
vent was unsuccessful because of the small solu- 
bility of sodium hydrosulphide. A substitution of 
ammonium hydrosulphide for the latter, however, 
despite its solubility in liquid hydrogen sulphide, 
showed no apparent reaction with the sulphonium 
bromide. Finally an attempt to prepare the 
analogous compound (C;H;);35—SC2H; by treat- 
ing sodium mercaptide with triethyl sulphonium 
bromide in absolute alcohol was without success. 
This inability to isolate the suggested compound 
indicates, in view of the existence of other parallel 
compounds, its marked instability, such that its 
immediate decomposition might be expected 


4 Pieroni and Coli, Gazz. chim. ital. 44, II, 353 (1914). 

5 Bernthsen and Sudborough, Organic Chemistry, pp. 
93, 752. D. Van Nostrand Co. (1930). 

®Crum-Brown and Blaikie, J. prakt. Chim. 23, 400 
(1881). 

™The experimental work was kindly performed by 
Dr. Nelson R. Trenner to whom thanks are due. 
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particularly under the experimental conditions 
used in the mercaptan decomposition studies. 
The third step in the series would then involve 
the decomposition of the sulphonium hydro- 
sulphide into propylene and hydrogen sulphide: 


(C3H7) 35 2H2S. 


The rapidity of this reaction with its large pres- 
sure increase would most certainly counter- 
balance the pressure decrease in the second step 
of the scheme, namely, the formation of the 
intermediate. No pressure decrease was ever 
observed. 

On the basis of this mechanism the production 
of polysulphides which are always found as a 
“tarry” residue receives a plausible explanation. 
The sulphonium hydrosulphide referred to above, 
being essentially a complex mercaptan can 
undergo the reaction in Step I to form complex 
sulphides by loss of hydrogen sulphide from two 
molecules. These complex sulphides in the pres- 
ence of more mercaptan or of hydrogen sulphide 
would unite to form new mercaptans still more 
complex, a cyclic chain of reactions which could 
produce large molecules of high sulphur content. 
There is evidence from the dynamic runs that 
such tarry materials contain both mercaptan and 
sulphide linkages since part is soluble in alkali 
and part insoluble. 

The absence of an induction period in the 
decompositions of propyl sulphide-hydrogen sul- 
phide and propyl sulphide-propyl mercaptan 
mixtures as also its presence in propyl mercaptan 
alone is well accounted for. The initial step in the 
latter involves the formation of propyl sulphide 
and hydrogen sulphide; a reaction with no 
volume change giving therefore an induction in 
the rate of pressure increase. With the mixtures 
the formation of the intermediate and _ its 
subsequent rapid decomposition can immedi- 
ately occur. In fact under the experimental con- 
ditions the rate of pressure increase must actually 
measure the rate of the bimolecular association 
with a pressure decrease in terms of the larger 
pressure increase of the much more rapid, 
probably unimolecular, decomposition of the 
intermediate. It is this counterbalancing of effects 
that probably gives to the fractional lives of the 
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later portions of the reaction the falling values 
with increasing pressure. 

The decreased induction periods in the in- 
jection tests also receive explanation in this same 
way. Since propyl sulphide or some substances 
with similar thio-ether linkages, possibly of the 
polysulphide type, are produced in the de- 
composition of the initial mercaptan sample, the 
rate of the association reaction, Step II, would be 
increased by the addition of the injected 
mercaptan, and hence the overall rate of pressure 
increase, giving a decreased induction for the 
injected sample. The suggested mechanism would 
necessitate that the concentration of propyl 
sulphide formed from propyl mercaptan by Step 
I should be greatest at the point of inflection. 
This is in line with the observation that the 
decrease in the induction period is greatest for the 
sample injected at the point of inflection. 

Finally, with regard to the equilibrium existing 

at the ‘‘end-point”’ of the decompositions, it has 
been pointed out that mixtures of propylene and 
hydrogen sulphide did not yield results in line 
with expectations from propyl mercaptan. This 
was shown to be due to the instability of 
propylene at these elevated temperatures. Duffey, 
Snow and Keyes have recently*® studied the 
equilibrium between isopropyl mercaptan, pro- 
pylene and hydrogen sulphide in presence of 
various catalysts at 275-300°C by a dynamic 
method. The rule of addition to unsaturated 
compounds would suggest that the —SH group in 
hydrogen sulphide would add to the least 
hydrogenated carbon atom in propylene thus 
giving isopropyl mercaptan. It is this same 
equilibrium then, that should eventually be 
reached in normal propyl mercaptan decompo- 
sition. Keyes and coworkers experienced diffi- 
culty in overcoming propylene polymerization in 
presence of many of their catalysts, such that the 
same end-point was not reached from the two 
sides of the equilibrium. The amount of isopropyl 
mercaptan present at equilibrium from the 
decomposition reaction was always greater than 
that from the combination reaction. 


® Duffey, Snow and Keyes, Paper presented at the 
Washington meeting of the Am. Chem. Soc. March, 1933. 
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To ascertain whether this propylene polymeri- 
zation would account for the difficulties in the 
present work, a quantity of propylene was 
introduced into the reaction vessel at 435°C. A 
very slow reaction with a pressure decrease 
occurred. After twenty-four hours this decrease 
having reached a maximum, the pressure began 
to increase indicating depolymerization and de- 
composition. This is in keeping with the findings 
of Wheeler and Wood® who find similar poly- 
merization followed by cracking. It would seem 
therefore, that a shift in the previously observed 
“end-point” of the propyl mercaptan decompo- 
sition should be noticeable after a long period of 
time due to the propylene reactions. To this end 
300 mm of mercaptan were maintained at 435°C. 
After half an hour the pressure had increased 244 
mm changing only one millimeter in the next 
half hour with no further change during the 
second hour. This is the value previously quoted 
as the ‘‘equilibrium end-point.” After seven hours 
the pressure had decreased five millimeters 
corresponding to propylene polymerization. No 


further change was noted after twelve hours but 


about twenty-four hours from the beginning of 
the reaction the pressure again slowly increased 
continuing thus for four days to apparent 
constancy at a value twenty-eight millimeters 
more than the previous ‘‘end-point,”’ thus con- 
firming the belief that the propylene decompo- 
sition was responsible for the lack of agreement in 
the equilibrium studies from the two sides of the 
reaction. From the experiments the decompo- 
sition rate is five times as rapid as the combi- 
nation reaction, consequently the time of contact 
of propylene-hydrogen sulphide mixtures must be 
five times that for the mercaptan thus giving 
opportunity for considerable propylene decompo- 
sition. This does not invalidate, however, the idea 
that the ‘end-points’ quoted in this paper 
represent pseudo-equilibrium values since all the 
reactions showed a constancy for at least two 
hours; it merely explains why the same “end- 
points” cannot be reached from the reverse side. 
The main steps therefore in the complete 
mechanism may be represented as follows: 


® Wheeler and Wood, J. Chem. Soc. 1819 (1930). 
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Step I C3H, H 
H Step III S 
Step II (C;H;)2S+H:S > SH 
C3H; SH C3H; 
C;H C3H and S 
and (C3H;)2S+C3;H;SH> S C3H; SH 


Step V 
The mathematical difficulties in the solution of 
the differential equations for the kinetics of such faster than either I or II, that Step I mainly 
a complex series of reactions prevent elucidation controls the length of the induction period and 
of individual rate constants. It may be noted, that Step II controls that portion of the reaction 
however, that Step III must be very considerably arising at the point of inflection. 
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The Dielectric Constant and Absorption of Several Organic Fluids at 1.82 m 


W. T. SzyMANowskI1, Department of Physics, University of Pittsburgh 
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For the further test of Debye’s theory of anomalous 
dispersion, dielectric constants and absorptions were meas- 
ured at \=1.82 m for the following fluids: ethyl alcohol, 
n-propyl alcohol, isopropyl alcohol, isobutyl alcohol, n- 
butyl alcohol, amyl alcohol, chlorobenzene, quinoline and 
nitrobenzene. The resonance method for undamped waves 
was used. The results are generally in good agreement with 
Debye’s theory as well as with the results of Mizushima 
and others for fluids whose polarizations for infinite dilu- 


tion do not differ greatly from those at normal concentra- 
tion, as for the alcohols and chlorobenzol. For fluids like 
quinoline and nitrobenzene, the agreement is far from 
satisfactory. For these, however, when Malsch modifica- 
tion of Debye’s theory is taken into account, a somewhat 
improved check is found, that for the dielectric constant 
being satisfactory, but that for the absorption of energy, 
while improved, being still far from satisfactory. 


INTRODUCTION 


EBYE’S dipole theory of anomalous dis- 

persion in the region of short electric waves 
has been lately confirmed by several authors. 
However, certain discrepancies have been found 
as in the case of glycerine by Bock ' and Mizu- 
shima.? These have been attributed to three 
possible causes. (1) Debye’s equations neglect 
the possible effects of associations of molecules. 
(2) The Debye equations are based on Mossotti’s 
hypothesis for the internal field. (3) The calcu- 
lation of the relaxation time is based on the 
assumption of a spherical molecule and therefore, 
on the validity of Stokes’ formula. The third 
assumption does not seem, however, to cause 
very great discrepancies between theory and 
experimental data since Johnston and Williams* 
obtained for nitrobenzene, a highly associated 
fluid, a value of the right order for the molecular 
radius calculated from the relaxation time 
obtained by extrapolation to infinite dilution. 

The purpose of this paper is the collec- 
tion of further data on dielectric constants for 
wave-lengths of the order of a few meters, 
especially for associated fluids for the sake of a 
further test of Debye’s dipole theory. 


*R. Bock, Zeits. f. Physik 31, 534 (1925). 

*San-Chiro Mizushima, Bull. Chem. Soc. Japan 1, 47, 
83, 115, 143, 163 (1926); Phys. Zeits. 28, 418 (1927). 

*J. H. L. Johnston and John Warren Williams, Phys. 
Rev. 34, 1483 (1929). 


EXPERIMENTAL METHOD 


A simple resonance method was used. Two 
UX 852 Radiotron tubes, as in Fig. 1, generated 
short electric waves. The resonance circuit con- 
sisted of a small loop, a variable condenser C,, 
and, connected in parallel with the latter, a 
Pyrex glass vessel C, with short platinum or 
tungsten lead-in wires, containing the liquid 
under investigation. The aperiodic circuit con- 
sisted of a loop of wire, a Weston thermocouple 
and a galvanometer. To eliminate any systematic 
error, preliminary measurements were made: 
(1) with two different condensers, one a variable 
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Fic. 1. Wiring diagram. 77, transformer; Ch, choke coils; 
Gr, grid leak; Th, thermocouple; C,, Pyrex vessel; C,, 
variable condenser. 
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glass condenser, the other a variable cylindrical 
air condenser. Both had linear characteristics and 
the two sets of measurements gave the same 
results within the limits of experimental error, 
which amounted to a few percent. (2) The 
measurements were made under different coup- 
ling conditions. (3) Vessels of different capacity 
were used, their capacities being as 1: 1.4: 2. 
In this way, it was possible to measure the 
dielectric constants of the investigated fluids 
through a wide range, and also to eliminate 
systematic errors by repeating single sets of 
measurements with different vessels. 
Calibration curves showing condenser setting 
for each vessel, as a function of the dielectric 
constant, were obtained with the aid of benzene- 
acetone and acetone-water solutions of known 
dielectric constants. These curves were approx- 
imately linear. A curve showing absorption of 
energy by the liquid in the cell C, of Fig. 1 was 
determined by a method similar to Mizushima’s 
in which KI solutions in acetone of known con- 
centration and conductivity were used. The 
deflections a of the galvanometer at the point of 
resonance were noted for each vessel containing 
- a fluid of a given conductivity at a given tem- 
perature. The ratios of these deflections to the 
deflections at resonance point without the fluid 
in the vessel k=a/ao were the ordinates of the 
absorption curves, whose abscissae were the cor- 
responding conductivities in ohm™ cm~. To 
correct for the apparent increase in capacity due 
to the conductivities of the fluids, correction 
curves for each vessel were drawn in a way 
similar to Mizushima’s. At the wave-length used 
in this investigation, acetone had not yet shown 
any effect of anomalous absorption at room tem- 
perature. Distilled water, however, showed a 
marked anomalous absorption.* A corresponding 
* Although this seems to be surprising, many measure- 
ments made in different coupling conditions and with 
acetone-water solutions of different concentrations showed 
always this absorption effect. A marked heating of distilled 
water was also observed at this frequency, if the vessel 
were put between two plates in a resonating circuit. 
That this absorption was not a normal absorption effect 
but an anomalous one must be concluded from the galva- 
nometer deflection data at different temperatures. These 
deflections showed an increase with temperature elevation 
which expressed in equivalent conductivity ranged from 


about 11.4 10-* to 6.5 X 10-* ohm~ cm~ for a temperature 
interval of about 15-30°C. 
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smaller effect was also observed in the water- 
acetone solutions used for the calibration of the 
vessels, therefore, the calibration curve of the 
vessels with greater capacity had to be corrected 
for the small apparent increase in capacity due 
to this effect. 

A wave-length of }\=1.82 m was chosen. By 
this choice, difficulties involved in using a 
thermostat with a resonating circuit for short 
electric waves were avoided, and data at a 
wave-length other than those already reported 
(Mizushima used \=3.08, 9.5, and 50 m) were 
obtained. Further, it was expected that the 
anomalous dispersion and absorption effects 
would be marked enough at this wave-length for 
the substances under investigation, without 
using extremely low temperatures. The wave- 
length was measured in a parallel wire Lecher 
system. The fluids were cooled and then allowed 
to warm up slowly. In the meantime, the tem- 
perature of the fluid ¢°, the condenser setting S, 
and the galvanometer deflection a were noted. 
The same observations were made during the 
cooling of the heated fluids. Before and after 
each set of measurements, and without fluid in 
the cell, the condenser setting Sp and the gal- 
vanometer deflection a were carefully deter- 
mined. Although the condenser allowed a sharp 
determination of the resonance point, the effects 
of the anomalous absorption in the investigated 
fluids flattened the resonance curve appreciably 
and caused errors ranging from one to three 
percent. This accuracy, however, was thought to 
be sufficient for the purpose of this investigation, 
since it was of the same order of Mizushima’s. 
The resonance curve did not have any irregu- 
larities. 


RESULTS 
Alcohols 


All alcohols used in this investigation were 
C.P. products of the Eastman Kodak Research 
Laboratories, with the exception of ethyl alcohol 
which was purified and distilled from a technical 
product. 


Ethyl alcohol 


Fig. 2 shows the dielectric constant and 
absorption curves for an interval of 0—50°C. 
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Fic. 2. Variation with temperature of dielectric constant 
) and absorption k& (----) for ethyl alcohol at 
\=1.82 m. 


These data may be compared in Table I with the 
values computed from Debye’s equations *4 


€1 — €9 
Cente. = eo Path (1) 
€ cate. = (€1—€0 L(er+2)/(eo+2) Jot 
1+[(a+2)/(eo+2) 
t=4rna®/KT 


where: 
T=absolute temperature; =viscosity; 


t°=temperature in °C; 


K=Boltzmann constant; 
t=relaxation time; w= X frequency; 
a=molecular radius; 
« =dielectric constant at high frequency =square of 
the refractive index for the D lines; 
«, =dielectric constant for zero frequency; 
€'ovs,=real part of the complex dielectric constant inter- 
polated from Fig. 2 at \=1.82 m; 
calc, = calculated value‘for the real part of the dielectric 
constant for \=1.82 m from Eq. (1); 


cops. =imaginary part of the complex dielectric constant 
=4na/w where o is the electrical conductivity 
corresponding, as shown by a calibration curve, 
to the observed ratio k of Fig. 2; 


e" calc, =imaginary part of the generalized dielectric con- 
stant calculated from Eq. (2). 


*P. Debye, Polar Molecules, p. 97, 1929, Chem. Cat. Co., 
New York, 
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TABLE I. Ethyl alcohol.* 
a=1.8X10-§ cm; \=1.82 m. 


” 
€ cale. 


” 
€cale. € obs. 


€'obs. 


€0 

50 0.00701 1.8 21.2 21.3 21.11 1.20 1.47 
40 0.00831 1.82 22.7 22.6 22.52 1.69 1.97 
30 0.00992 1.83 24.25 23.9 23.93 2.57 2.74 
20 0.0119 1.85 25.8 25.15 25.15 4.48 3.80 
10 0.0145 1.86 27.0 26.45 25.86 6.88 5.24 
J 1.88 27.72 26.28 9.10 7.10 


* Debye’s Eqs. (1) and (2) were used in obtaining e'catc, 
and eé’’eaic, in Tables I through IV. 


The computations were made under the 
assumption of a molecular radius a=1.8X10-* 
cm. The values for « and ¢; at 30°, 40°, and 
50°C were extrapolated from Mizushima’s data. 
The experimental values of 1s. and k at 0°C 
were also extrapolated from the curves in Fig. 2. 
The same extrapolation procedure was used in 
all the following tables. From Table I, we can 
see that the theoretical values check well the 
experimental ones as was found by Mizushima 
for other wave-lengths. However, considerable 


‘discrepancies are noted when the computations 


are carried out under the assumption of a=2.0 
10-8 cm as calculated from critical point data,° 
the van der Waals constant } being taken as 
16/37Na’. 


n-Propyl alcohol 


The absorption curve shows a maximum at 
around 5°C. The dielectric constant starts to 
decrease in the region of 20°C (Fig. 3). For an 
assumed molecular radius a=2.210-* cm as 
calculated from critical data, a comparison of 
experimental and theoretical values is shown in 


TABLE II. n-Propyl alcohol. 
a=2.2X10-§ cm; A=1.82 m. 


(poise) €0 € obs. calc. 
50 0.011305 1.87 18.35 16.2 27.74 1.80 2.96 
40 0.01405 1.89 19.65 16.85 18.57 2.95 4.26 
30 0.0178 1.90 20.9 17.42 18.82 4.48 5.67 
20 0.0226 1.92 22.2 17.4 18.35 7.54 7.92 
10 0.0292 1.94 23.5 15.35 16.73 10.1 9.97 
0 0.0388 1.95 24.8 12.6 13.73 10.6 11.42 


5 Landolt and Bornstein, Tables, 5th Edition. VI. 
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Table II. The check is not as good as for ethyl 
alcohol. The general characters of the theoretical 
and experimental curves, however, are the same. 


at: 


050 


_ Fic. 3. Variation with temperature of dielectric constant 
e ( : - and absorption k (----) for n-propyl alcohol at 
A=1.82 m. 


10 20 


0 lo 20 to 
Fic. 4. Variation with temperature of dielectric constant 
e ( , a and absorption k (----) for isopropyl alcohol at 
A= m, 
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case. 


Isopropyl alcohol 

We do not have sufficient data for «9 and e¢; 
at different temperatures to compare the theo- 
retical and the experimental values; however, 
since the viscosity of isopropyl alcohol is but 
little different from that for n-propyl alcohol, 
we expect, according to Debye’s theory, the 
maximum of absorption and the decrease of the 
dielectric constant to be in the same region as 
for n-propyl alcohol. Fig. 4 shows this to be the 


Isobutyl alcohol 

In Fig. 5 we find a maximum for the absorption 
at about 15°C and a decrease of the dielectric 
constant at about 30°C. Calculating the molec- 


TABLE III. Isobutyl alcohol. 
a=2.3X10-§ cm; A=1.82 m. 


e ( 
A=1.82 m. 


n 

(poise) €cale. € obs. € calc. 
50 0.0161 1.88 17.4 13.25 16.19 2.14 4.12 
40 0.0212 1.90 18.25 13.7 15.92 3.54 5.71 
30 0.0287 1.93 19.15 13.6 14.74 5.56 7.51 
20 0.0391 1.95 20.0 11.9 12.33 8.00 8.92 
10 0.0555 1.97 21.0 8.6 8.92 8.45 9.20 
0 0.0804 1.99 21.8 5.4 5.88 6.82 7.88 
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Fic. 5. Variation with temperature of dielectric constant 
) and absorption & (----) for isobutyl alcohol at 
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ular radius from critical data, we obtain 
a=2.3X10-* cm. With this value, Table III 
was computed and the experimental values 
compared with the theoretical ones. In general 
the check is good except for the higher tem- 
peratures. 


n-Butyl alcohol 


Here again we do not have sufficient data for 
€) and ¢; at different temperatures in order to 
compare theoretical and experimental values as 
for isobutyl alcohol; however, since the viscosity 
of n-butyl alcohol is less than that of isobutyl 
alcohol, we expect a shift of both the absorption 
and the dielectric constant curves towards lower 
temperatures. This is the case, as may be seen 
from a comparison of Figs. 5 and 6. 


Fic. 6. Variation with temperature of dielectric constant 
. ae, and absorption k (----) for n-butyl alcohol at 
=1.82 m. 


1 


Amyl alcohol 


Here is greater than that of the above- 
mentioned alcohols; therefore, a shift towards 
higher temperatures for the maximum of ab- 
sorption and for the bend in the dielectric con- 
stant curve must be expected. Fig. 7 shows this 
to be the case. The maximum absorption occurs 
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Fic. 7. Variation with temperature of dielectric constant 
4 C a and absorption k (----) for amyl alcohol at 
=1.82 m. 


at about 20°C and the bend in the dielectric 
constant curve starts at about 40°C. We are not 
in possession of sufficient critical data in order 
to compute the molecular radius. It is expected, 
however, to be much greater than that of the 
above-mentioned alcohols. Assuming a=2.5 
X10-* cm, Table IV was computed. The check 
is, in general, good. 


TABLE IV. Amzyl alcohol. 
a=2.5X10-§ cm; A=1.82 m. 


€' obs. 


11.6 


n 
(poise) 


0.01864 
0.02434 
0.03235 


13.85 
14.6 
15.3 
16.0 
16.8 
17.5 


Chlorobenzene 


A C.P. product from the Eastman Kodak 
Research Laboratories was used. For chloro- 
benzene the molecular radius computed from 
critical data is a=2.4X10-* cm. From Eq. (1) 
we compute for 18°C and for }=1.82 m 
€'calc. = 5.7. On the other hand, the experimental 


COTA 
4 5 
ac 
13 ~ 80 
12 Fond 70 
Bt ee 
| 
‘ 
50 1.91 12.58 2.13 3.69 
J + 40 1.94 11.55 12.16 3.60 4.95 
30 1.96 10.7 11.00 5.46 6.24 
20 0.0437 1.98 8.4 9.02 6.84 7.00 
10 0.0610 2.00 6.1 6.56 6.11 6.84 
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value is €’ obs. =5.61. The discrepancy is still not 
great. No anomalous absorption was observed 
at this temperature. 


Quinoline 

A C.P. product from the Eastman Kodak 
Research Laboratories was used. For quinoline 
sufficient critical data are not known for the 
computation of the molecular radius. We have 
here to compute it from the experimental die- 
lectric constant in a manner similar to that used 
by Mizushima for glycerine and to consider 
whether or not this value is reasonable. Putting 
in Eq. (1) the observed experimental value of 
the dielectric constant for quinoline at 18°C, 
\=1.82 m; namely, ¢’ obs. =8.74, we obtain as the 
molecular radius a=2.3X10-* cm. This value 
seems to be small, since the benzene molecule 
has a radius of about 2.3 X 10-* cm and quinoline 
has a benzene and a pyridine ring. We can see, 
therefore, that for such a degree of association 
as is in quinoline, the simple Debye equations 
do not check the experimental data. Quinoline 
shows a small anomalous absorption corre- 
sponding to a conductivity of about 4x10-° 
ohm-! at about 10°C. 


Nitrobenzene 


Nitrobenzene was investigated in two different 
samples. Sample 1 was a C.P. Baker Company 
product, sample 2 was obtained by double dis- 
tillation of Baker’s technical nitrobenzene. The 
experimental values for both samples checked 
well within the limits of the experimental error 
which was about 1 percent. Fig. 8 shows the 
curves for dielectric constants and absorptions 
at different temperatures. Table V, based on a 


TABLE V. Nitrobenzene.* 
a=2.4X10-§ cm; A\=1.82 m. 


(poise) €0 € 


°C e‘obs. €’cale. € €’cale. 
60 0.0109 2.36 28.61 28.58 26.46 0.71 v4 
50 0.0125 2.37 30.03 29.6 26.67 0.82 8.3 
40 0.0144 2.38 31.56 30.8 26.37 1.07 11.1 
30 0.0168 2.39 33.35 32.3 25.50 1.40 13.5 
20 0.0201 2.41 35.46 34.2 23.27 1.96 15.95 


* Debye’s Eqs. (1) and (2) were used in obtaining ¢’catc. 
and e’caic.. (See Table VI.) 
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Fic. 8. Variation with temperature of dielectric constant 
e ( ) and absorption k (----) for nitrobenzene at 
A\=1.82 m. 


value of a@=2.4X10-* cm for the molecular 
radius of nitrobenzene according to Johnston 
and Williams,’ shows great discrepancies between 
experimental and theoretical values. For this 
highly associated fluid, Eqs. (1), (2) and (3) are 
not sufficient. Taking the experimental values 
for the dielectric constant of nitrobenzene, we 
can compute again the molecular radius from Eq. 
(1). This is found to be a= 1.53 X 10-8 cm at 20°C. 
The value computed this way is too small, in 
agreement with what was found above for 
quinoline and by Mizushima for glycerine. 


DISCUSSION OF RESULTS 


The foregoing shows that as long as the molar | 


polarizations at normal concentration and at 


infinite dilution are not very different as in © 


the case of alcohols and chlorobenzene, Eqs. | 
(1), (2) and (3) yield, at least qualitatively, — 


values for the real and the imaginary parts of 


dielectric constants in agreement with experi-. 


ment. This conclusion checks not only with the 
results of Mizushima, but also with the more 
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recent* results of Goldhammer,® Malsch,’? and 
Girard and Abadie* who investigated alcohols 
and alcohol solutions. The latter, however, were 
tested at different wave-lengths from those used 
in the present investigation. For fluids where 
molar polarization at infinite dilution P,, differs 
very much from that at normal concentration 
P,, as in the case of quinoline, where this ratio 
P.,/P.=1.7 (Rolinski’), and nitrobenzene, where 
P../Pn=3.6, and possibly glycerine, Eqs. (1) and 
(2) need a correction. At first, it seemed that 
such a correction for associations could be easily 
made in accord with Debye’s theory of associ- 
ation,” where the molecules are assumed to 
produce associated molecular groups, having a 
smaller or a larger electric moment. This theory 
was developed by Wolfke'! and the theoretical 
molar polarization curves as a function of the 
concentration for different alcohols agreed well 
with the experimental ones. For chlorobenzene, 
quinoline, and nitrobenzene the computation of 
the degree of association was in good accordance 


_ with the experimental data. If we want to apply, 


however, this theory to measurements at a high 
frequency, we must introduce in the equations 
for the molar polarization, P, a frequency factor 
1/(1+iwt) Debye (p. 90). 

The equation will then be: 


e—1M 4r 


—=— Na, 


e+2 p 3 
1 
+— No ( : ) (4) 
3 3KT 


P\.) =molar polarization at a given frequency; 
M=molecular weight; 
p=density; 
N=total number of molecules per mole; 
N.=number of unassociated molecules; 
«=molecular polarizability taking account of the dis- 
tortion effect; 
u=electric moment; 
w= 27 X frequency. 


(o) = 


where: 


*The present investigation was made in 1931. Since 
that time, several papers appeared dealing with a similar 
subject. 

*Goldhammer, Phys. Zeits. 33, 361 (1932). 

7 Malsch, Ann. d. Physik 12, 865 (1932); Phys. Zeiis. 33, 
19 (1932), 

“ — and Abadie, Comptes Rendus 195, 119, 217 
2). 


Deducing the expression for « we can find, 
however, that this substitution does not change 
the final expressions (1) and (2), as found by 
Debye. The above theory of association cannot, 
therefore, explain the discrepancies encountered 
at a higher frequency. 

Malsch"? recently proposed a correction of the 
Debye equations, based on the introduction of an 
additional internal force, as already used by 
Lorentz.!* This internal force is due to the inter- 
action of the molecules and, therefore, can be 
considered the result of association. Such an 
assumption brings the problem back to a cor- 
rection of the Clausius Mossotti relation, which 
can be then deduced from: 


where 


F=total force acting on a molecule per unit charge; 
E=field strength; 

P=polarization =electric moment per unit volume; 
v=factor for the introduction of an additional force. 


We have then according to Malsch:!? 
(42/3)NaT 


e—1M 
e+2 p obs. 1—v(p/M)(41/3)Na 


where a=ao+p?/3KT. Having the experimental 
values for the molar polarization at normal con- 
centration and at infinite dilution, we can find 
v as a function of the temperature. The cor- 
rection factor in Debye’s Eqs. (1) and (2) would 
then be: 


é+2+(e—1)p 
2= 


instead of (€:+2)/(e+2). Finally the Debye 
equations will have the form: 


= eo + (e:—€0)/[1+ (zt)? ] 
— egzwt/[1+ ]. 


(5) 
(6) 


If we now apply this new equation to nitro- 


® Rolinski, Phys. Zeits. 29, 658 (1928). 

10 Debye, Handbuch d. Radiologie, Vol. 6. 

1 Wolfke, Phys. Zeits. 29, 713 (1928); Comptes Rendus 
de la Soc. Pol. de Phys. 

12 Malsch, Phys. Zeits. 33, 383 (1932). 

13 Lorentz, Wied. Ann. 9, 641 (1880). 
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TABLE VI. Nitrobenzene.* 
a=2.4X10-§ cm; A\=1.82 m. 


obs. €’cale, 


28.58 28.26 
29.60 29.47 
30.8 30.58 
32.3 31.99 
34.2 33.11 


* Debye-Malsch Eqs. (5) and (6) were used in computing 
€'cale. and (See Table W.) 


benzene, we see from Table VI that the agree- 
ment between experimental and_ theoretical 
values is much better for the dielectric constant. 
For the absorption, the discrepancy is still very 
great, although it is better than that in Table V. 
In Tables V and VI, the molecular radius was 
assumed to be 2.4X10-* cm which is a rather 
small value, since x-ray studies indicate that it 
should be about 3.0 10-§ cm—4.0X 10-8 cm in 
which case the correction used above would still 
be insufficient. If this correction were applied to 
the alcohols, then all values would be over- 
corrected unless we would assume higher radii. 
For instance, for ethyl alcohol, with a radius of 
2.0 10-* cm and the correction of Malsch, we 
again find a good agreement between theoretical 
and experimental data. u-Propyl alcohol, how- 
ever, would seem to be overcorrected unless a 
higher radius were taken. The corrections for 
isobutyl alcohol, amy] alcohol, and chlorobenzene 
are negligible. Applying the correction to quino- 
line, we find a much bigger radius, that of 2.6 
X 10-8 cm instead of 2.3 10-8 cm as found with- 
out correction. The bigger radius certainly seems 
more reasonable. From the above, however, we 
can see that without a knowledge of the exact 
molecular radius, no quantitative comparison as 
to the validity of the introduced corrections can 
be made, since even a small change in the radius 
will change considerably the computed results. 
As it was pointed out by Oncley and Williams,“ 
the radii can be several times larger, at low 
temperatures, than kinetic theory values because 
of association effects. This would bring the associ- 
ation phenomena back to the Debye association 


4 Oncley and Williams, Phys. Rev. 43, 341 (1933). 


theory” mentioned above. However, only com- 
paratively small parts of the alcohol molecules 
would have a double moment, as can be seen 
from the comparison of polarization at infinite 
dilution and of that at normal concentration. 
Therefore the correction of the average relaxation 
time would be small. In a triple molecule with 
a zero resultant moment there would be no influ- 
ence on the relaxation time. 

Since there are, at present, insufficient data 
available for molecular radii at different tem- 
peratures for associated fluids, we have to use 
the kinetic theory values as a first approxima- 
tion, or preferably x-ray data, so far as they are 
known. 

Another way of explaining the discrepancy 
between the theoretical and experimental data 
would be that of Oncley and Williams.'* On the 
basis of many experiments with solutions of 
various polar substances in mineral oils of vary- 
ing viscosities, they found that there was no 
dispersion even at 1.727X10’ cycles. These 
results are consistent with the experiments of 
Goldhammer‘® and Weigle and Liithi.'® 

Oncley and Williams give a probable explana- 
tion of the results by the failure of Stokes’ law. 
The reason for such a result may be either 
because the application of a law in a form of 
Stokes’ law is not justified, or because the use of 
macroscopic viscosity is not permissible. The 
experiments of Oncley and Williams on floricin 
and rosin solutions seem to indicate that Stokes’ 
law is true in its general form. Therefore, we 
must assume that for the calculation of the 
relaxation time of an associated fluid, we have to 
use a corrected macroscopic viscosity. 

Such a correction has been used by Gold- 
hammer,® who found a much better agreement 
between experimental and theoretical data for a 
solution of n-butyl alcohol in paraffin oil at 
\=0.72 m if the viscosity were decreased 15 
percent. 

The same procedure applied to nitrobenzene 
would, in our case, give a higher value for the 
molecular radius of nitrobenzene than 2.410 
cm which is rather small. 


% Weigle and Liithi, Comptes Rendus 49, 130 (1932). 
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€’’obs, 
60 0.71 2.9 
50 0.82 3.96 
40 1.07 5.0 
7 30 1.40 6.3 
20 1.96 8.4 
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